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1 Introduction 

A major success of low temperature physics was achieved with the intro- 
duction by Landau of the notion of quasiparticles. According to his hypoth- 
esis, the properties of many body interacting systems at low temperatures 
are determined by the spectrum of some low energy long living excitations 
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(quasiparticles). Another milestone of many body theory is the Mean Field 
Approximation (MFA), which permitted achieving considerable progress in 
the theory of phase transitions. Phenomena which cannot be described by 
the quasiparticle method or by MFA arc usually called fluctuations. The 
BCS theory of superconductivity is a bright example of the use of both the 
quasiparticle description and MFA. The success of the BCS theory for tra- 
ditional superconductors was determined by the fact that fluctuations give 
small corrections with respect to the MFA results. 

During the first half of the century after the discovery of superconduc- 
tivity the problem of fluctuation smearing of the superconducting transition 
was not even considered. In bulk samples of traditional superconductors the 
critical temperature T c sharply divides the superconducting and the normal 
phases. It is worth mentioning that such behavior of the physical character- 
istics of superconductors is in perfect agreement both with the Ginzburg- 
Landau (GL) phenomenological theory (1950) (!]] and the BCS microscopic 



The characteristics of high temperature and organic superconductors, low 
dimensional and amorphous superconducting systems studied today, strongly 
differ from those of the traditional superconductors discussed in textbooks. 
The transitions turn out to be much more smeared out. The appearance of 
superconducting fluctuations above critical temperature leads to precursor 
effects of the superconducting phase occurring while the system is still in 
the normal phase, sometimes far from T c . The conductivity, the heat capac- 
ity, the diamagnetic susceptibility, the sound attenuation, etc. may increase 
considerably in the vicinity of the transition temperature. 

The first numerical estimation of the fluctuation contribution to the heat 
capacity of a superconductor in the vicinity of T c was done by Ginzburg in 
1960 [g. In that paper he showed that superconducting fluctuations increase 
the heat capacity even above T c . In this way fluctuations change the temper- 
ature dependence of the specific heat in the vicinity of critical temperature, 
where, in accordance with the phenomenological GL theory of second order 
phase transitions, a jump should take place. The range of temperatures where 
the fluctuation correction to the heat capacity of a bulk clean conventional 
superconductor is relevant was estimated by Ginzburg [] to be 



where Ep is the Fermi energy. The correction occurs in a temperature range 
ST many orders of magnitude smaller than that accessible in real experiments. 

1 The expression for the width of the strong fluctuation region in terms of the 
Landau phenomenological theory of phase transitions was obtained by Levanyuk 
|Q. So in the modern theory of phase transitions the relative temperature width 
of fluctuation region is called the Ginzburg-Levanyuk parameter Gi(D), where D 
is the effective space dimensionality. 





(1) 
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In the 1950s and 60s the formulation of the microscopic theory of su- 
perconductivity, the theory of type-II superconductors and the search for 
high-T c superconductivity attracted the attention of researchers to dirty sys- 
tems, superconducting films and filaments. In 1968, in the papers of L. G. 
Aslamazov and A. I. Larkin ||, K.Maki || and a little later in the paper of 
R.S.Thompson |7| the fundaments of the microscopic theory of fluctuations 
in the normal phase of a superconductor in the vicinity of the critical tem- 
perature were formulated. This microscopic approach confirmed Ginzburg's 
evaluation Q for the width of the fluctuation region in a bulk clean supercon- 
ductor. Moreover, was found that the fluctuation effects increase drastically 
in thin dirty superconducting films and whiskers. In the cited papers was 
demonstrated that fluctuations affect not only the thermodynamical proper- 
ties of superconductor but its dynamics too. Simultaneously the fluctuation 
smearing of the resistive transition in bismuth amorphous films was experi- 
mentally found by Glover 0, and it was perfectly fitted by the microscopic 
theory. 

In the BCS theory Q only the Cooper pairs forming a Bose-condensate 
are considered. Fluctuation theory deals with the Cooper pairs out of the 
condensate. In some phenomena these fluctuation Cooper pairs behave simi- 
larly to quasiparticles but with one important difference. While for the well 
defined quasiparticle the energy has to be much larger than its inverse life 
time, for the fluctuation Cooper pairs the "binding energy" Eq turns out to 
be of the same order. The Cooper pair life time tql is determined by its 
decay into two free electrons. Evidently at the transition temperature the 
Cooper pairs start to condense and tql = oo. So it is natural to suppose 
from dimensional analysis that tql ~ h/ks(T — T c ). The microscopic theory 
confirms this hypothesis and gives the exact coefficient: 

= Wh (2) 
TGL 8k B (T-T c y [ ' 

Another important difference of the fluctuation Cooper pairs from quasi- 
particles lies in their large size £(T). This size is determined by the distance 
on which the electrons forming the fluctuation Cooper pair move away during 
the pair life-time tql- In the case of an impure superconductor the electron 
motion is diffusive with the diffusion coefficient T> ~ v 2 f t (t is the electron 
scattering time^j), and £,d{T) — \JVtgl ~ vf\/ttgl- In the case of a clean 
superconductor, where fesTr 3> fr, impurity scattering does not affect any 
more the electron correlations. In this case the time of electron ballistic mo- 
tion turns out to be less than the electron-impurity scattering time r and is 
determined by the uncertainty principle: r^ai ~ h/ksT. Then this time has 

2 Strictly speaking r in the most part of future results should be understood as 
the electron transport scattering time Ttr- Nevertheless, as it is well known, in 
the case of isotropic scattering these values coincide, so for sake of simplicity we 
will use hereafter the symbol r. 
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to be used in this case for determination of the effective size instead of t: 
£ C (T) ~ VFy/hTaL/ksT. In both cases the coherence length grows with the 
approach to the critical temperature as (T — T c )~ 1//2 , and we will write it 
down in the unique way (£ = £ c ,d): 

f T — T 

The microscopic theory in the case of an isotropic Fermi surface gives for £ 
the precise expression: 

„9 KpT 2 f , .1 H . ,,1. ft ,',El 

^) = "ir {^2 + " ^ ~ 4^2V* ( 2 } ) ' (4) 

where ip(x) is the digamma function and D = 3, 2, 1 is the space dimension- 
ality. In the clean (c) and dirty (d) limits: 

Cc = 0.133^ JI = 0.7^ JI, (5) 




Here / = uj?t is the electron mean free path and £o = hvp/irA^) is the 
conventional BCS definition of the coherence length of a clean superconductor 
at zero temperature. One can see that (||) and @ coincide with the above 
estimations [J 

Finally it is necessary to recognize that fluctuation Cooper pairs can be 
really treated as classical objects, but these objects instead of Boltzmann 
particles appear as classical fields in the sense of Ray leigh- Jeans. This means 
that in the general Bose-Einstein distribution function only small energies 
E(p) are involved and the exponent can be expanded: 



3 



Let us stress some small numerical difference between our Exp. (j4j) and the usual 
definition of the coherence length. We are dealing near the critical temperature, 
so the definition (^) is natural and permits us to avoid many numerical coeffi- 
cients in further calculations. The cited coherence length £o = hvF /nA(0) — 
0.18hvF/kBT c , as is evident, was introduced for zero temperature and an 
isotropic 3-D superconductor. 

It is convenient to determine the coherence length also from the formula for 
the upper critical field: H c2 {T) = A(T)$ /2n£ 2 (T). A(T C ) = 1, while its value at 
T = depends on the impurities concentration. For the dirty case the appropriate 
value was found by K.Maki Q Ad(0) = 0.69, for the clean case by L.Gor'kov |§ 
A 2 C D (0) = 0.59, A 3 C D (0) = 0.72. 
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This is why the more appropriate tool to study fluctuation phenomena is not 
the Boltzmann transport equation but the GL equation for classical fields. 
Nevertheless at the qualitative level the treatment of fluctuation Cooper pairs 
as particles with the density np) = J n(p) , 2 often turns out to be useful 

& 

Below will be demonstrated both in the framework of the phenomeno- 
logical Ginzburg-Landau theory and the microscopic BCS theory that in the 
vicinity of the transition 

E(p) = ak B (T T c ) + JE^ = ^ [hye (T) + P 2 ] ■ (8) 

Far from transition temperature the dependence n(p) turns out to be more 
sophisticated than (^), nevertheless one can always write it in the form 

n( P ) = ^e(T)f{^). (9) 

In classical field theory the notions of the particle distribution function 
n(p) (proportional to E^ 1 (p) in our case) and Cooper pair mass m* are 
poorly determined. At the same time the characteristic value of the Copper 
pair center of mass momentum can be defined and it turns out to be of 
the order of pa ~ h/£(T). So for the combination m*E(po) one can write 
m*E(j>o) ~ Pq ~ h 2 /^ 2 (T). In fact the particles density enters into many 
physical values in the combination n/m*. As the consequence of the above 
observation it can be expressed in terms of the coherence length: 

^ = r ~ k -fe- D (n do) 

m m E(po) a n 

Pq here estimates the result of momentum integration. 

For example we can evaluate the fluctuation Cooper pairs contribution to 
conductivity by using the Drude formula 

a = ^ ^d D ^e- D (T)(2efr GL {e) ~ (11) 

m ft 



4 This particle density is denned in the (D)- dimensional space. This means that it 
determines the normal volume density of pairs in 3D case, the density per square 
unit in 2D case and the number of pairs per unit length in ID. The real three 
dimensional density n can be defined too: n = d D_3 n( D ), where d is the thickness 
of the film or wire. 
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Analogously a qualitative understanding of the increase in the diamagnetic 
susceptibility above the critical temperature may be obtained from the well- 
known Langevin expression for the atomic susceptibility^]: 

X = A—* ( R2 ) => -^^d D -^- D (T) ~ -e D ' 2 - 2 . (12) 
c m c n 

Besides these examples of the direct influence of fluctuations on supercon- 
ducting properties, indirect manifestations by means of quantum interference 
in the pairing process and of renormalization of the density of one-electron 
states in the normal phase of superconductor take place. These effects, be- 
ing much more sophisticated, have a purely quantum nature, and in contrast 
to the direct Cooper pair contributions require microscopic consideration. 
This is why, in developing phenomenological methods through the first five 
Sections of this review, we will deal with the direct fluctuation pair contri- 
butions only. The sixth Section is devoted to the microscopic justification of 
the time dependent Ginzburg-Landau equation. The description of the micro- 
scopic theory of fluctuations, including indirect fluctuation effects, discussion 
of their manifestations in various physical properties of superconductors will 
be given in the Sections 7-8. 

The first seven Sections are written in detail, so they can serve as a 
textbook. On the contrary, in the last (eighth) Section a wide panorama 
of fluctuation effects in different physical properties of superconductors is 
presented. Thus this Section has more of a handbook character, and the 
intermediate calculations often are omitted. 

Finally we would like to mention that the number of articles devoted to 
superconducting fluctuations published in the last 33 years is of the order of 
ten thousands, so our bibliography list does not pretend nor at the complete- 
ness nor at the establishment of the rigorous priorities. 



2 Ginzburg-Landau formalism. Thermodynamics 
2.1 Fluctuation contribution to heat capacity 

GL functional. The complete description of the thermodynamic proper- 
ties of a system can be done through the exact calculation of the partition 
function^]: 

I / ■Q\ I 

(13) 

5 This formula is valid for the dimensionalities D — 2, 3, when fluctuation Cooper 
pair has the possibility to "rotate" in the applied magnetic field and the average 
square of the rotation radius is < R 2 >~ tf(T). "Size" effects, important for low 
dimensional samples, will be discussed later on. 

6 Hereafter K = Ub = c = 1. 
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As discussed in the Introduction, in the vicinity of the superconducting 
transition, side by side with the fermionic electron states, fluctuation Cooper 
pairs of a bosonic nature appear in the system. As already mentioned, they 
can be described by means of classical bosonic fields ^(r) which can be treated 
as " Cooper pair wave functions ". So the calculation of the trace in ([l3]) can be 
separated into a summation over the " fast " electron degrees of freedom and 
a further functional integration carried out over all possible configurations of 
the "slow" Cooper pairs wave functions: 



D&(r)Z[V(r)], (14) 



where 



^(r)]=exp(-^M) (is) 

is the system partition function in a fixed bosonic field ^(r), already summed 
over the electronic degrees of freedom. Here it is assumed that the classical 
field dependent part of the Hamiltonian can be chosen in the spirit of the GL 
approach^ 

F[¥(r)]=F N + J rfV- |«|^(r)| 2 + ||^(r)| 4 + _l_| W( r)| 2 | .. (16) 

Let us discuss the coefficients of this functional. In accordance with the 
Landau hypothesis, the coefficient a goes to zero at the transition point and 
depends linearly on T — T c . Then a = aT c e; all the coefficients a, b and m 
are supposed to be positive and temperature independent. Concerning the 
magnitude of the coefficients it is necessary to make the following comment. 
One of these coefficients can always be chosen arbitrary: this option is re- 
lated to the arbitrariness of the Cooper pair wave function normalization. 
Nevertheless the product of two of them is fixed by dimensional analysis: 
ma r~j £~ 2 (T). Another combination of the coefficients, independent of the 
wave function normalization and temperature, is a 2 /b. One can see that it 
has the dimensionality of the density of states. Since these coefficients were 
obtained by a summation over the electronic degrees of freedom, the only 
reasonable candidate for this value is the one electron density of states v. 
The microscopic theory gives the precise coefficients for the above relations: 

AmaT c = £~ 2 ; a 2 /b = ^-v, (17) 

where C,(x) is the Riemann zeta function, £(3) = 1.202. One can notice that 
the arbitrariness in the normalization of the order parameter amplitude leads 



For simplicity the magnetic field is assumed to be zero. 



Fluctuation Phenomena in Superconductors 9 

to the unambiguity in the choice of the Cooper mass introduced in ( |l6| ) as 
2m. Indeed, this value enters in ( p"^ ) in the product with a so one of these 
parameters has to be fixed. In the case of a clean D-dimensional supercon- 
ductor it is natural to suppose that the Copper pair mass is equal to two free 
electron masses what results in 

2Dtt 2 T c . , 

a ^=7CME- F (18) 

As the first step in the Landau theory of phase transitions & is supposed 
to be independent of position. This assumption in the limit of sufficiently 
large volume V of the system permits a calculation of the functional integral 
in (|l4|) by the method of steepest descent. Its saddle point determines the 
equilibrium value of the order parameter 

l^ 2 = {of>o 6,e< °- ( 19 ) 



Choosing a in accordance with (|18|) one finds that this value coincides with 
the superfluid density n s of the microscopic theory || . 

The fluctuation part of the free energy related to the transition is deter- 
mined by the minimum of the functional (Jlq) : 



p = (mUu = m = \l N ~ rt^ 6 < ■ (20) 




From the second derivative of ( pp| ) one can find an expression for the jump 
of the specific heat capacity at the phase transition point: 



"-^-to-nwj-nw)- (2i) 

d' 2 F\ a\ 8,' 

= -r-Tc = —7—s vl r . 



VT C \ de 2 J b 7C(3) 

Let us mention that the jump of the heat capacity was obtained because 
of the system volume was taken to infinity first, and after this the reduced 
temperature e was set equal to zero. 



Zero dimensionality: the exact solution. In a system of finite volume 
fluctuations smear out the jump in heat capacity. For a small superconducting 
sample with the characteristic size d -C S,(T) the space independent mode 
&o = ^\/V defines the main contribution to the free energy: 
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Fig. 1. Temperature dependence of the heat capacity of superconducting grains in 
the region of the critical temperature 



Z (0) = 



(a|^ | 2 + & |-^o| 4 ) 



TT 3 VT 

exp(x 2 )(l-erf(x))l =a ^. (22) 



2b 

By evaluating the second derivative of this exact result jllj] one can find the 
temperature dependence of the superconducting granular heat capacity (see 
Fig. |l|). It is evident that the smearing of the jump takes place in the region 
of temperatures in the vicinity of transition where x ~ 1, i.e. 



where we have supposed that the granule is clean; T c0 and £ are the critical 
temperature and the zero temperature coherence length (see the footnote 3 
in Introduction) of the appropriate bulk material. From this formula one can 
see that even for granule with the size d ~ £o the smearing of the transition 
still is very narrow. 

Far above the critical region, where GiiQ) <C e <C 1, one can use the 
asymptotic expression for the erf(x) function and find 

F (0) =-TlnZ (0) =-Tln— . (23) 

ae 
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Calculation of the second derivative gives an expression for the fluctuation 
part of the heat capacity in this region: 

SC (0) = 1/e 2 . (24) 

The experimental study of the heat capacity of small Sn particles in the 
vicinity of transition was done in p2[ 



Arbitrary dimensionality: case T > T c . It is possible to estimate the 
fluctuation contribution to the heat capacity for a specimen of an arbitrary 
effective dimensionality on the basis of the following observation. The volume 
of the specimen may be divided into regions of size £(T), which are weakly 
correlated with each other. Then the whole free energy can be estimated as 
the free energy of one such zero-dimensional specimen (p3|), multiplied by 
their number N {D) = V£~ D (T) : 



F [D) =~TVC D (T)ln— . (25) 
ae 



This formula gives the correct temperature dependence for the free energy 
for even dimensionalities. A more accurate treatment removes the In e depen- 
dence from it in the case of the odd dimensions. 

Let us begin with the calculation of the fluctuation contribution to the 
heat capacity in the normal phase of a superconductor. We restrict ourselves 
to the region of temperatures beyond the immediate vicinity of transition, 
where this correction is still small. In this region one can omit the fourth 
order term in ^(r) with respect to the quadratic one and write down the GL 
functional, expanding the order parameter in a Fourier series: 

F[m =F N + J2l*+ ^Hl^kl 2 = aT c £ (e + £ 2 k 2 ) |^ k | 2 . (26) 

k k 

Here $k = J !f - (r)e _4kr (iV" and the summation is carried out over the 
vectors of the reciprocal space. Now we see that the free energy functional 
appears as a sum of energies of the independent modes k. The functional 
integral for the partition function ([15]) can be separated to a product of 
Gaussian type integrals over these modes: 



(27) 



Carrying out these integrals, one gets the fluctuation contribution to the free 
energy: 

F(6>0) = -Tln^ = -r^ln * k2 (28) 

k a(e + 4maTc ) 



12 Anatoly Larkin and Andrei Varlamov 



The appropriate correction to the specific heat capacity of a superconductor 
at temperatures above the critical temperature may thus be calculated. We 
are interested in the most singular term in e _1 , so the differentiation over the 
temperature can be again replaced by that over e : 

1 ( d 2 F\ 1 ^ 1 

^-^y = 4 ( C+ ^ - (29) 

The result of the summation over k strongly depends on the linear sizes 
of the sample, i.e. on its effective dimensionality. As it is clear from (p9|). 
the scale with which one has to compare these sizes is determined by the 
value (4maT c e)"5 which, as was already mentioned above, coincides with 
the effective size of Cooper pair £(T). Thus, if all dimensions of the sample 
considerably exceed £(T), one can integrate over (2ir)~ 3 L x L y L z dk x dk y dk z 
instead of summing over n x ,n y ,n z . In the case of arbitrary dimensionality 
the fluctuation correction to the heat capacity turns out to be 

sc + = -JVT^t^ = «°-^^- (30) 

where Vd = V, S, L, 1 for D = 3,2, 1, 0. For the coefficients "do it is convenient 
to write an expression valid for an arbitrary dimensionality D, including 
fractional ones. For a space of fractional dimensionality we just mention that 
the momentum integration in spherical coordinates is carried out according 
to the rule: / d D k/ (2n) D =fJ, D J k°- 1 dk, where 

MO = 2 D n D /ir(D/2 + 1) (31) 



and r(x) is a gamma-function. The coefficient in (|30|) can be expressed in 
terms of the gamma-function too: 

_ r(2-£>/2) 

yielding &i = 1/4, i? 2 = 1/4tt and i9 3 = 1/8tt. 

In the case of small particles with characteristic sizes d < £(e) the ap- 
propriate fluctuation contribution to the free energy and the specific heat 
capacity coincides with the asymptotics of the exact results ( ^3|) and (^4|). 
From the formula given above it is easy to see that the role of fluctuations in- 
creases when the effective dimensionality of the sample or the electron mean 
free path decrease. 
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Arbitrary dimensionality: case T < T c . The general expressions ( |l4"| ) 
and (|l^) allow one to find the fluctuation contribution to heat capacity below 
T c . For this purpose let us restrict ourselves to the region of temperatures 
not very close to T c from below, where fluctuations are sufficiently weak. In 
this case the order parameter can be written as the sum of the equilibrium 
& (see (pj|)) and fluctuation tp(r) parts: 

W{t)=9 + iI3{t). (33) 

Keeping in ([l6]) the terms up to the second order in ip( r ) an d up to the fourth 
order in & , one can find 



Z[9] = exp( ±-L ) J] / rfRe^kdWk x (34) 

k ^ 

f 1 k 2 ~ k 2 

x exp <^ -~[(3b& 2 + a + — )Re 2 ^k + (b^ 2 + a+ — Jlm 2 ^" 
[_ J 4to 4m 

Carrying out the integral over the real and imaginary parts of the order 
parameter one can find an expression for the fluctuation part of the free 
energy: 




(35) 



Let us discuss this result. It is valid both above and below T c . The two terms 
in it correspond to the contributions of the modulus and phase fluctuations 
of the order parameter. Above T c \P = and these contributions are equal: 
phase and modulus fluctuations in the absence of \P represent just two equiv- 
alent degrees of freedom of the scalar complex order parameter. Below T c , the 
symmetry of the system decreases (see (|l9|)). The order parameter modulus 
fluctuations remain of the same diffusive type as above T c , while the char- 
acter of the phase fluctuations, in accordance with the Goldstone theorem, 
changes dramatically. 

Substitution of to ([35]) results in the disappearance of the temper- 
ature dependence of the phase fluctuation contribution and, calculating the 
second derivative, one sees that only the fluctuations of the order parame- 
ter modulus contribute to the heat capacity. As a result the heat capacity, 
calculated below T c , turns out to be proportional to that found above: 

5C- = 2%- 2 5C+. 



Hence, in the framework of the theory proposed we found that the heat 
capacity of the superconductor tends to infinity at the transition temperature. 
Strictly speaking, the restrictions of the above approach do not permit us to 
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discuss seriously this divergence at the critical point itself. The calculations in 
principle are valid only in that region of temperatures where the fluctuation 
correction is small. We will discuss in the next Section the quantitative criteria 
for the applicability of this perturbation theory. 



2.2 Ginzburg-Levanyuk criterion 

The fluctuation corrections to the heat capacity obtained above allow us to 
answer quantitatively the question: where are the limits of applicability of 
the GL theory? 

This theory is valid not too near to the transition temperature, where 
the fluctuation correction is still small in comparison with the heat capacity 
jump. Let us define as the Ginzburg-Levanyuk number Gi^r>) P fl th e value 
of the reduced temperature at which the fluctuation correction (|30[) equals 
the value of AC Mf\ 



G 



a 



Vn d R—i 

V -£-§ D b{±m)-T c - 



(36) 



Substituting into this formula the microscopic values of the GL theory pa- 
rameters (n7h one can find 



Gi{D) = 



8tt 2 



Vd 
V 



1 



(37) 



Since vnT c ~ vdVf/£,c ~ Pp 1 £ c 1 ~ 0,1 D £c 1 one can convert this formula 
to the form 



G 



1(D) 



7((3)$ D (Vd\ Zc.a n 11 
V 



8tt 2 



where a is the interatomic distance, ft is worth mentioning that in bulk 
conventional superconductors, due to the large value of the coherence length 
(£ c ~ 10 -6 -j- 10 _4 cm), which drastically exceeds the interatomic distance 
(a ~ 10~ 8 cm), the fluctuation correction to the heat capacity is extremely 
small. However, the fluctuation effect increases for small effective sample di- 
mensionality and small electron mean free path. For instance, the fluctuation 
heat capacity of a superconducting granular system is readily accessible for 
experimental study. 

8 One can see that some arbitrariness occurs in this definition. For instance the Gi 
number could be defined as the reduced temperature at which the AL correction 
to conductivity is equal to the normal value of conductivity (as it was done 
in |l3|,|l4j]). Such definition results in the change of the numerical factor in Gi 
number: 

Gi(2,cr) = l-44Gl(2,h. c .). 
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Using the microscopic expression for the coherence length the Ginzburg- 
Levanyuk number (^) can be evaluated for different cases of clean (c) and 
dirty (d) superconductors of various dimensionalities and geometries (film, 
wire, whisker and granule are supposed to have 31? electronic spectrum): 



Gi (3) 






G«( ) 


80(|^) 4 ,(c) 


(&) ■ W 


0.5, (c) 

1.3( P 2,5)" 2/3 (T c t)" 1/3 , (d) wire 
2.3 (p%S) ~ 2/3 , (c) whisker 


\/ 7 C( 3 ) i _ 

2ir y/vT c V 

i ■? "iTsO-x la 



Table 1. 



One can see that for the three dimensional clean case the result coincides 
with the original Ginzburg evaluation and demonstrates the negligibility of 
the superconducting fluctuation effects in clean bulk materials. 



2.3 Scaling and renormalization group 

In the above study of the fluctuation contribution to heat capacity we have 
restricted ourselves to the temperature range out of the direct vicinity of the 
critical temperature: |e| > Gi^y As we have seen the fluctuations in this 
region turn out to be weak and neglecting their interaction was justified. In 
this Section we will discuss the fluctuations in the immediate vicinity of the 
critical temperature ( |e| < Gi/jj\) where this interaction turns out to be of 
great importance. 

We will start with the scaling hypothesis, i.e. with the belief that in the 
immediate vicinity of the transition the only relevant length scale is £(T). The 
temperature dependencies of all other physical quantities can be expressed 
through £(T). This means, for instance, that the formula for the fluctuation 
part of the free energy (p5| ) with the logarithm omitted is still valid in the 
region of critical fluctuations^] 

F (D) ~ -r%), (38) 

the coherence length is a power function of the reduced temperature: £(e) ~ 
e~' y . The corresponding formula for the fluctuation heat capacity can be 
rewritten as 

sc ~_4r?~^- 2 . (39) 

As was demonstrated in the Introduction, the GL functional approach, 
where the temperature dependence of £(T) is determined only by the dif- 
fusion of the electrons forming Cooper pair, £(e) ~ e -1 / 2 and SC ~ e -1 / 2 . 

9 The logarithm in ( pa ) is essential for the case D = 2. This case will be discussed 
later. 
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Fig. 2. Examples of diagrams for the fluctuation contribution to b. 

These results are valid for the GL region ( |e| > Gi) only, where the inter- 
action between fluctuations can be neglected. In the immediate vicinity of 
the transition (so-called critical region), where |e| < Gi, the interaction of 
fluctuations becomes essential. Here fluctuation Cooper pairs affect the co- 
herence length themselves, changing the temperature dependencies of £(e) 
and (5C(e). In order to find the heat capacity temperature dependence in the 
critical region one would have to calculate the functional integral with the 
fourth order term, accounting for the fluctuation interaction, as was done 
for OD case. For the 3-D case up to now it is only known how to calculate 
a Gaussian type functional integral. This was done above when, for the GL 
region, we omitted the fourth order term in the free energy functional (|l6|). 

The first evident step in order to include in consideration of the critical 
region would be to develop a perturbation series in b Any term in this 
series has the form of a Gaussian integral and can be represented by a dia- 
gram, where the solid lines correspond to the correlators (&(r)iP*(r )\. The 
"interactions " b are represented by the points where four correlator lines 
intersect (see Fig. ||). 

This series can be written as 



For e > Gi it is enough to keep only the first two terms to reproduce the per- 
turbational result obtained above. For e < Gi all terms have to be summed. 
It turns out that the coefficients c„ can be calculated for the space dimen- 
sionality D — > 4 only. In this case the complex diagrams from Fig. || (like the 
diagram similar to an envelope) are small by the parameter e = 4 — D and 

10 Let us mention that this series is an asymptotic one, i.e. it does not converge 
even for small b. One can easily see this for small negative b, when the integral 
for the partition function evidently diverges. This is also confirmed by the exact 
OD solution (E|). 



4--P 
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in order to calculate c n it is sufficient to sum the relatively simple "parquet" 
type diagrammatic series. Such a summation results in the substitution of the 
"bare" vertex b by some effective interaction b which diminishes and tends to 
zero when the temperature approaches the transition point. Such a method 
was originally worked out in quantum field theory [ p^|p^ , p7| | . For the problem 
of a phase transition such a summation was first accomplished in p8| . 

Instead of a direct summation of the diagrams it is more convenient and 
physically obvious to use the method of the renormalization group. In the case 
of quantum field theory this method was known long ago |19| , ^of . For phase 
transition theory it was proposed in [|l8|]2lf| but the most simple and evident 
formulation was presented by Wilson |22| . The idea of the renormalization 
group method consists in separating the functional integration over "fast" 
{^\k\>A) an d "slow" (V'lfc^yi) fluctuation modes. If the cut off A is large enough 
the fast mode contribution is small and the integration over them is Gaussian. 
After the first integration over fast modes the functional obtained depends 
on the slow ones only. They can, in their turn, be divided on slow (\k\ < Ai) 
and fast (A\ < \k\ < A), and the procedure can be repeated. Moving step by 
step ahead in this way one can calculate the complete partition function. 

As an example of the first step of renormalization, the partition function 
calculation below T c can be recalled. There we separated the order parameter 
into the space-independent part & ("slow" mode) and the fluctuation part 
ip(r) ("fast" mode) which was believed to be small in magnitude. Being in 
the GL region it was enough to average over the fast variables just once, 
while in the critical region the renormalization procedure requires subsequent 
approximations . 

The cornerstone of the method consists in the fact that in the critical 
region at any subsequent step the free energy functional has the same form. 
For D close to 4 this form coincides with the initial free energy GL functional 
but with the coefficients a A and b A depending on A. We will perform these 
calculations by the method of mathematical induction. Let us suppose that 
after the (n — 1)— st step the free energy functional has the form: 



^[^U-J = F N,A n -i + f dv[a An _ 1 \W An _ 1 \ 2 + (40) 



Writing i ?a„_ 1 in the form ^a„_! = &A„ + Vm„ and choosing A n close enough 
to A n -i it is possible to make ipA n so small, that one can restrict the func- 
tional to the quadratic terms in ipA only and perform the Gaussian inte- 
gration in complete analogy with (|34|) . The important property of the spaces 
with dimensionalities close to 4 is the possibility to choose A n <C A n —\ and 
still to have ipA <C &A n - In this case 'P Ari can be taken as coordinate inde- 



pendent, and one can use the result directly following from (34) 
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t:T c 



-*rV (41) 



7rT r 

4m 



Expanding the last term in (41) in a series in \&a„ one can get for F^P An ] the 
same expression ( |40| ) with the substitution of /l„_i — > vl n . From ( f4l| ) follows 
that 



yl„<|k|<A„_ 1 

&A„ 



a A „ = + 2T £ (42) 



^A = b A , — 5T A " : ; . 

A„<|k|<A n _ 1 



Passing to a continuous variable vl n — * A one can rewrite these recursion 
equations as the set of differential equations: 

g^=-I>^ln , , (43) 



D — 1 



~SA-- bl ^ D (a(A) + A± } 2 



SAT ~ - 2T VD- 

' " 4 W . (44) 



These renormalization group equations are evidently valid for small enough 
A only, where the transition from discrete to continuous variables is justified. 
This means at least 

A 2 /Am^T c0 Gi {D) . (45) 

in order to move away from the first approximation. 

Let us recall that in the framework of the Landau theory of phase tran- 
sitions the coefficient a(T c o) = at the transition point and this can be 
considered as the MFA definition of the critical temperature T c $. The same 
statement for the function a — a(A) in the framework of the renormalization 
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group method can be written as the a(T c o,A ~ = 0. With the decrease 
of A the effect of critical fluctuations more and more is taken into account and 
the renormalized value of the critical temperature decreases, being defined by 
the equation: a(T c (A),A) = 0. Finally, after the application of the complete 
renormalization procedure, one can define the real critical temperature T c , 
shifted down with respect to T c q due to the effect of fluctuations, from the 
equation: 

a(T c ,A = 0) = 0. (46) 

It is easy to find this shift in the first approximation. Indeed, let us integrate 
the first equation of (Q) over A in limits [0, The main contribution to 
the integral will be determined by the region where a(A) <C Being far 
from the critical point one can assume that the coefficient b — const and then 



fair 1 ) /•!/? 

a(C 1 ) = aST c = da = -8mTu D b A 

Ja(0) JO 



D - 3 dA. 



For the 3D case this gives the shift of the critical temperature 8T C due to 
fluctuations [^] 

<5T C (3) 2mb b 1 7C(3) 



T c nail 2irT c a 2 £ 3 167r 3 z/T c £ 3 



-jGi^y (48) 



Let us come back to study of properties of the system of equations (| 
One can find its partial solution at T — T c in the form: 



a ^= , m[ ^-D)] A2 (49) 

b (T c ,A) = 1 4 ~ D ^ A 4 - D . 

y ° ; \&m 2 Tn D [5 + (4 - D)] 2 

These power solutions are correct in the domain of validity of the system ( |44"|) 
itself, i.e. for small enough A defined by the condition ((45[). Nevertheless, in 
a space of dimensionality close to 4 (D = 4 — e,e <C 1) it is possible to 
extend their validity up to the GL region and to observe their crossover to 

11 In the same way we can analyze the shift of the critical temperature in 2D case 
too and obtain: 



As we will show below both 3D and 2D results for ST C coinside with those 
obtained by the analysis of the effect of fluctuations on superconducting density 
in the perturbation approach. 
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the GL results: a(T c ) = 0; b(T c ) = bo = const. Indeed, in this case, due to 
proportionality of a(A) to e — ► 0, one can omit it in the denominators of the 
system (Q) and to write down the solution for b(A) in the form 

b~\T c , A) = + ^^TM^-* - (50) 

We have chosen the constant of integration as &q 1 — Tfi]j£ 4 ~ D in order 

to match the renormalization group and GL solutions at the value of A — 
A ~ 

J1 max s 

Now let us pass to study of the function a(T, A) for the same interesting 
case of space dimensionality D — > 4 for temperatures slightly different ( but 
still close enough) from T c , where one can write 

a{T, A) = a(T c , A) + a(T c , A)T c e. 

The first term on the right hand side is determined by Eq. (|49|) . In order to 
determine a(T c ,A) let us expand the first equation in (H) in terms of e 



da(T c ,A) _ orr b(T c ,A)A D - 1 
For 



— - 2Tn D i ' A2 , 2 a(T c , A). (51) 

OA n (T A\-lAL\ 



A z /4m>a(T c ,A)T c e (52) 

we can again use the solution ( |50| ) for b(T c , A) and omit a(T c , A) in the 
denominator o f (|5l| ). The constant of integration, appearing in the process 
of solution of (|5l|), is chosen in accordance with the condition that for A = 
A max ~ we match a(T c , A) = a(T c o, £ _1 ) = ao with the GL theory: 



a(T c , A) = a Q 



80m 2 b D _ 4 ,_ D 

1 + Ja^d) T ^ d{a ~t ] 



-2/5 



(53) 



The condition (p2h can be written as A > £ _1 (T), where £(T) is the 
generalized coherence length, determined by the equation: 

r 2 (T)=4ma (T c , £ _1 (T)) T c e. 

Such a definition is valid at any temperature. For example, far enough from 
the critical point, in the GL region, a (T C ,£ -1 (T)) = ao and one reproduces 
the result (|^). Vice versa, in the critical region the main contribution on the 
right hand side of the Eq. (|5^) results from the second term containing A D ~ 4 
so, putting £ -1 (T) = A, one can rewrite the self-consistent equation for £(T) 
and get 

f(T) = (4m)^>/ 2 L_ e -», (54) 

20o i ^DV-l c c<o 
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where 2u = [1 — (4 — I?)/5] _1 . As was already mentioned, strictly speaking 
this result was carried out for e = 4 — D -C 1, so it is confident up to the 
first in e expansion only: v = 1/2 + e/10. Nevertheless extending it to e = 1 
(D = 3) one can obtain 1/3 = 3/5. 

Let us pass to the calculation of the critical exponent of the heat capacity 
in the immediate vicinity of the transition. For this purpose one can calculate 
the second derivative of equation (E3h with respect to e: 



dC(A) _ 2 A D ~W(A) 



T , , JW1 — ; A2 , (55) 



The heat capacity renormalized by fluctuations has the value C(A = 0) which 
is the result of integration over all fluctuation degrees of freedom. Carrying 
out the integration of ([35|) over all A < one can divide the domain of 
integration on the right hand side in two: A < £ _1 (T) and £ _1 (T) < A < 
In the calculation of the integral over the region £ > A > £ _1 (T) the 
inequality a(A)T c e<§; 4— holds, and the function a(A) can be omitted in the 
denominator. In the numerator of (^5|) one can use for a(A) the solution ([53]). 
In the region A < £ _1 (^) one nas to use the partial solution (Ej?) for a(A) 
and can find that the contribution of this domain has the same singularity 
as that from the region A > £ _1 (T), but with a coefficient proportional to 
(4 — D) 2 — e 2 , hence negligible in our approximation. The result is: 



C(A = 0) =al[{AmTf^ D ]f 



4 (4 - D) 

5 b 



4/5 



4—D 



4-D 



(T). (56) 



Substituting the expression for £(T) one can finally find 



C = 2^al[^ D -^±-,]h-^ (57) 



m 2 T 2 
'bW~D) 

confirming the validity of the scaling hypothesis and the relation (|39|). The 
critical exponent in (15^) is 



' 4 - D > -./10. 



10[1- (4-D)/5] 



One can see that generally speaking the critical exponents ^ and a appear 
in the form of series in powers of e. More cumbersome calculations permit 
finding the next approximations for them in e = 4 — D. Nevertheless it is 
worth mentioning that even the first approximation, giving ^3 = 3/5 and 
o?3 = 1/10 for e = 1, is already weakly affected by the following steps of the 
expansion in powers of e p3[ . 

One can notice that the exercise performed in this Section has more aca- 
demic than practical character. Indeed, the results obtained turn out appli- 
cable to the analysis of the critical region of a 3D superconductor only if 
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Gi is so small that the theoretical predictions are hardly experimentally ob- 
servable. Nevertheless we demonstrated the RG method which helps to see 
the complete picture of the fluctuations manifestation in the vicinity of the 
critical temperature. 



2.4 Fluctuation diamagnetism 

Qualitative preliminaries. In this Section we discuss the effect of fluctua- 
tions on the magnetization and the susceptibility of a superconductor above 
the transition temperature. Being the precursor effect for the Meissner dia- 
magnetism, the fluctuation induced magnetic susceptibility has to be a small 
correction with respect to the diamagnetism of a superconductor but it can 
be comparable to or even exceed the value of the normal metal diamagnetic 
or paramagnetic susceptibility and can be easily measured experimentally. 
As was already mentioned in the Introduction the temperature dependence 
of the fluctuation induced diamagnetic susceptibility can be qualitatively an- 
alyzed on the basis of the Langevin formula, but some precautions in the case 
of low dimensional samples have to be made. 

As regards the 3D case we would like just to mention here that Exp.(|l2|), 
presented in terms of £(T), has a wider region of applicability than the GL 
one. Namely, the scaling arguments are valid for diamagnetic susceptibility 
too and one can write the general relation 

X(3) ~ -e 2 T£(T) ~ -xpe- 1 ' 2 { ' /T/2-* < n . , (58) 

which is valid in the region of critical fluctuations in the immediate vicin- 
ity of the transition temperature too. Here, in order to define the scale of 
fluctuation effects, we have introduced the Pauli paramagnetic susceptibility 
Xp — e 2 VF I '47T 2 . Moreover, the Langevin formula permits us to extend the 
estimation of the fluctuation diamagnetic effect to the other side beyond the 
GL region: to high temperatures T ^> T c . The coherence length far from the 
transition becomes a slow function of temperature. In a clean superconductor, 
far from T c , £(T) ~ vf/T, so one can write 

X(3c) (T » T c ) ~ -e 2 T£(T) ~ - Xp (59) 

and see that the fluctuation diamagnetism turns out to be of the order of the 
Pauli paramagnetism even far from the transition. More precise microscopic 
calculations of X(3)C^ ^ ^c) l ea d to the appearance of \n 2 (T/T c ) in the 
denominator of fl59). 

In the 2D case Exp.(p^) is applicable for the estimation of X(2) m the case 
when the magnetic field is applied perpendicular to the plane, permitting 2D 
rotations of fluctuation Cooper pairs in it: 

X(2c) (T) ~ e 2 — < i? 2 >^ e 2 T£ 2 (T) ~ - X p . (60) 
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This result is valid for a wide range of temperatures and can exceed the Pauli 
paramagnetism by factor ^f- even far from the critical point (we consider the 
clean case here). 

For a thin film (d <C £,(T)) perpendicular to the magnetic field the fluc- 
tuation Cooper pairs behave like effective 2D rotators, and the formula ( p"2| ) 
still can be used, though one has to take into account that the susceptibility 
in this case is calculated per unit square of the film. So for the realistic case 
(from the experimental point of view) of the dirty film, one has just use in 
( |60| ) the expression (|J) for the coherence length: 

e 2 T 2 f l\ T r 



XM~^?W~-X P [ 1 ) ¥ =t o . (61) 

Let us discuss now the important case of a layered superconductor (for 
example, a high temperature superconductor). It is usually supposed that 
the electrons move freely in conducting planes separated by a distance s. 
Their motion in the perpendicular direction has a tunneling character, with 
effective energy J. The related velocity and coherence length can be estimated 
as v z = dE(p)/dp± ~ J/p± ~ sJ and £z,(c) ~ sJ/T for clean case. In dirty 
case the anisotropy can be taken into account in the spirit of formula (^|) 
yielding ~ y/DjT ~ S J^/7jT. 

We start from the case of a weak magnetic field applied perpendicular 
to layers. The effective area of a rotating fluctuation pair is £e(e)fi/(e)- The 



density of Cooper pairs in the conducting layers (10) has to be modified for 



the anisotropic case. Its isotropic 3D value is proportional to l/£(e), that 
now has to be read as ~ 1/ y/£x{e)£ y (e)- The anisotropy of the electron 
motion leads to a concentration of fluctuation Cooper pairs in the conduct- 
ing layers and hence, to an effective increase of the Cooper pairs density of 
yj £x ( e )£» ( e ) /£* ( e ) times its isotropic value. This increase is saturated when 
£ 2 (e) reaches the interlayer distance s, so finally the anisotropy factor appears 
in the form \/£,x(£)£y(t) / max{s, £ z (e)} and the square root in its numerator 
is removed in the Langevin formula (fL2), rewritten for this case 



X(la yer ,±)(e,H^ 0) ~ -e 2 T (62) 

max{s,^ z (e)| 

The existence of a crossover between the 2D and 3D temperature regimes 
in this formula is evident: as the temperature tends to T c the diamagnetic 
susceptibility temperature dependence changes from 1/e to lf-y/e . This hap- 
pens when the reduced temperature reaches its crossover value e cr = r 
(C^( £ cr) ~ s). The anisotropy parameter 

4^(0) J 2 f *f, Tt«1 

plays an important role in the theory of layered superconductors |^[ 

12 We use here a definition of r following from microscopic theory (see Section 6). 
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It is interesting to note that this intrinsic crossover, related to the spec- 
trum anisotropy, has an opposite character to the geometric crossover which 
happens in thick enough films when £(T) reaches d . In the latter case the 
characteristic 3D 1 / y/e ~ dependence taking place far enough from T c (where 
£(T) <C d), is changed to the 2D 1/e law (see (|6l"|)) in the immediate vicin- 
ity of transition ( where £(T) ^> d) ^4|. It is worth mentioning that in a 
strongly anisotropic layered superconductor the fluctuation-induced suscep- 
tibility may considerably exceed the normal metal dia- and paramagnetic 
effects even relatively far from T c p5| , pc| . 

Let us consider a magnetic field applied along the layers. First it is nec- 
essary to mention that the fluctuation diamagnctic effect disappears in the 
limit J ~ £ 2 — > 0. Indeed, for the formation of a circulating current it is 
necessary to tunnel twice, so 

X(layer,\\) ~ ~e 2 T mjJx{ ^ (r)} ~ ~Xp( — ) ^^{^ J /T} ' 

In the general case of an anisotropic superconductor, choosing the z axis 
along the direction of magnetic field H, the following extrapolation of the 
results obtained may be written 

x ~-e 2 T ^(<0£„(e) ^ 

max{a, £ x (e)} max{b, £ y (e)} max{s, 6(e)} 

This general formula is useful for the analysis of the fluctuation diamagnetism 
of anisotropic superconductors or samples of some specific shape: granular, 
quasi-l-D, quasi-2£>, and 313. It is also applicable to the case of a thin film 
(d <C £z( e )) placed perpendicular to the magnetic field: it is enough to replace 
£z( £ ) by d in (|64|). Nevertheless the formula (|(34| ) cannot be applied to the 
cases of thin films in parallel fields, wires and granules. In those cases the 
Langevin formula ( |l2|) can still be used with the replacement of (-R 2 ) — * d 2 



X(D) XP 



— ) e- D d D - X ~ e D ' 2 -\ 



v F 



The magnetic field dependence of the fluctuation part of free energy in 
these cases is reduced only to account for the quadratic shift of the critical 
temperature versus magnetic field. 

For 3D systems or in the case of a film in a perpendicular magnetic field 
the critical temperature depends on H linearly, while the magnetic field de- 
pendent part of the free energy for H <C H* 2 (—e) (the line H* 2 (—e) is mirror- 
symmetric to the H C 2(e) with respect to y-axis passing through T — T c ) is 
proportional to H 2 . This is why the magnetic susceptibility is determined by 



Eq.(64) for weak enough magnetic fields H <C ^o/[£x(e)£y( e )] = H c2 {e) <C 



H C 2 (0) only. In the vicinity of T c these fields are small enough. 
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Zero-dimensional diamagnetic susceptibility. For quantitative analysis 
of the fluctuation diamagnetism we start by writing down the GL functional 
for the free energy (see Exp.(jl6|)) in the presence of the magnetic field 



F[#(r)]=F n + J dV^(r)\ 2 + ^(r)\ 4 + ^-\(-t\7-2eA)W(r)\ 



2 



where A is vector potential. As long as fluctuation effects are comparatively 
small, the average magnetic field in the metal B may be assumed to be equal 
to the external field H. Thus we omit the last two terms in ( |65| ) (see later 
on). 

The fluctuation contribution to the diamagnetic susceptibility in the sim- 
plest case of a "zero-dimensional" superconductor (spherical superconducting 
granule of diameter d <C £(e)) was considered by V.Shmidt [jllj. In this case 
the order parameter does not depend on the space variables and the free 
energy can be calculated exactly for all temperatures including the critical 
region in the same way as was done for the case of the heat capacity in the 
absence of a magnetic field. Formally the effect of a magnetic field in this 
case is reduced to the renormalization of the coefficient a, or, in other words, 
to the suppression of the critical temperature. This is why one can use the 



same formula (22) for the partition function with the critical temperature T c 
shifted by magnetic field asp} 

T c (ff)=T c (0)(l-^<A 2 >). (66) 

Here = § is the magnetic flux quantum and < > means the averaging 

over the sample volume. 

Such a trivial dependence of the properties of OD samples on magnetic 
field immediately allows one to understand its effect on the heat capacity 
of a granular sample. Indeed, with the growth of the field the temperature 
dependence of the heat capacity presented in Fig. [ijjust moves in the direction 
of lower temperatures. 

In the GL region Gi(o) < e one can write the asymptotic expression ( p3[ ) 
for the free energy: 

F {0) (e,H) = -Tin- ^ — -. 

<*{e+%$- < A 2 >) 

In the case of a spherical particle the relation < A 2 >= j^H 2 d 2 can be used 
in full analogy with the calculation of the moment of inertia of a solid sphere. 

13 Let us stress the difference between the H 2 shift of the critical temperature for 
a zero-dimensional granule and the linear shift in the case of bulk material. 
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In this way an expression for the OD fluctuation magnetization valid for all 
fields H -C H C 2(0) can be found: 

2£!£ ,2 



■V/> 



dF (0) (e,H) 

L(n\\t, U ) — — — 1 ttzt\ 



M (o) (C # ) = ^77^ = ~T , . & (67) 



One can see that the fluctuation magnetization turns out to be negative 
and linear up to some crossover field, which can be called the temperature 
dependent upper critical field of the granule i? C 2(o)( e ) ~ df^ry = §#c2(0)-\/e 
at which it reaches a minimum. At higher fields H c2 ^(e) < H <C H C 2(0) the 
fluctuation magnetization of the OD granule decreases as 1/7?. In the weak 
field region H <C -ff C 2(o) ( e ) the diamagnetic susceptibility is: 

which coincides with our previous estimate in its temperature dependence 
but the numerical factor found is very large. Let us underline that the tem- 
perature dependence of the OD fluctuation diamagnetic susceptibility turns 
out to be less singular than the OD heat capacity correction: e™ 1 instead of 
e- 2 . 

The expression for the fluctuation part of free energy ( p8| ) is also applicable 
to the cases of a wire or a film placed in a parallel field: as was already 
mentioned above all its dependence on magnetic field is manifested by the 
shift of the critical temperature (|66|). In the case of the wire in a parallel 
field one has to choose the gauge of the vector-potential A=|Hxr yielding 
(A 2 ) ( , tuire ,. . = (the calculation of this average is analogous to that of 

the moment of inertia of a solid sphere). For a wire in a perpendicular field, or 
a film in a parallel field, the gauge has to be chosen in the form A =(0, Hx, 0) 
(to avoid the appearance of currents perpendicular to surface). One can find 
(A 2 ), . . , = ^ri- for a wire and (A 2 ), ... ... = for a film. 

\ / (tmre,-L) lb \ / ( f il m.W ) 12 

Calculating the second derivative of Eq.(|2g|) with the appropriate mag- 
netic field dependencies of the critical temperature one can find the following 
expressions for the diamagnetic susceptibility: 



-4=, wire in parallel field 



X(r>)(e) = -27T— xp < wire in perpendicular field . (68) 



■gg In j , film in parallel field 



GL treatment of fluctuation magnetization. Let us analyze quantita- 
tively, on the basis of the GL functional, the temperature and field depen- 
dencies of the fluctuation magnetization. We will carry on the discussion for 
a layered superconductor. As was already mentioned this system has a great 
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practical importance because of its direct applicability to high temperature 
superconductors, where the fluctuation effects are very noticeable. Moreover, 
the general results obtained will allow us to analyze 31? and 2D situations 
as limiting cases. The effects of a magnetic field are more pronounced for 
perpendicular orientation, so let us consider first this case. 

The generalization of the GL functional for a layered superconductor 
(Lawrence-Doniach (LD) functional p7|) in a perpendicular magnetic field 
can be written as 

T LD [#]=Y,J d2 r (a\# l \ 2 + ~m i + ±\(V l] -2ieA ]] )\* l \ 2 

+J'|*i+i-*i| 2 ) , (69) 

where \Pi is the order parameter of the l—th superconducting layer and the 
phcnomenological constant J is proportional to the Josephson coupling be- 
tween adjacent planes. The gauge with A z = is chosen in (|6|). In the 
immediate vicinity of T c the LD functional is reduced to the GL one with 
the effective mass M = (AJs 2 )^ 1 along c-direction, where s is the inter-layer 
spacing. One can relate the value of J to the coherence length along the 
c-direction: J = 2aT c (l/s 2 . Since we are dealing with the GL region the 
fourth order term in (|6^) can be omitted. 

As it is well known the Landau representation is the most appropriate for 
problems related with the motion of a charged particle in a uniform magnetic 
field. The fluctuation Cooper pair wave function <fi n k z (r) can be written as 
the product of a plane wave propagating along the magnetic field direction 
and a Landau state wave function. Let us expand the order parameter ^(r) 
on the basis of these eigenfunctions: 

Wi{t) = ^n,fc>»fe,(r)exp(i^O, (70) 

n,k z 

where n is the quantum number related with the degenerate Landau state 
and k z is the momentum component along the direction of the magnetic field. 
Substituting this expansion into (|6^) one can find the LD free energy as a 
functional of the ^ n ,fe s coefficients: 

V- f H ( 1\ , „ ssl , ,9 



Tld [9 n , km ] = 2^ [ aT - € + [ n + 2 ) +J t 1 - cos ( fc , s ))| l -( 71 ) 

In complete analogy with the case of an isotropic spectrum the functional 
integral over the order parameter configurations x P n _k z in the partition func- 
tion can be reduced to a product of ordinary Gaussian integrals, and the 
fluctuation part of the free energy in a magnetic field takes the form: 

F[6 ' H) = _ ^ T S ln -T c e+ ^-An + % + J(l-cos(k z s)y i72) 
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n)= I /(»)~24 /(^ + 1/2)-/ (-1/2) 



Here the summation over the degenerate states of each Landau level was 
performed (S is the sample cross-section) and results in appearance of the 
number of particle states (2nHS/<Po) with the definite quantum numbers n 
and k z . The summation over n has to be performed through all occupied 
states, i.e. the upper limit of the sum is N ~ Itxi^qEf j H. 

In the limit of weak fields one can carry out the summation over the 
Landau states by means of the Euler-Maclaurin's transformation 

N fN+l/2 i 

£/( ' ' 

n=0 J' 1 /* 

and obtain 

, . irSTH 2 r /s Nsdk z f 1 1 ,„ s 

Here M is the total number of layers. After the momentum integration one 
gets: 

F (e ,H)=F M ) 1 



m 



with the anisotropy parameter defined ai 

r = 2 J = 4g(0) 
aT s 2 

The magnetic susceptibility in a weak field turns out |2q,P9| to be 



(74) 



e 2 T £ 2 
3tts + r ) 



These results confirm the qualitative estimation (|62|) additionally providing 
the exact value of the numerical coefficient and the temperature dependence 
in the crossover region. In the limit r 3> e Exp.(|75|) transforms into the 
diamagnetic susceptibility of the 31? anisotropic superconductor pOI] . 



14 Let us stress the difference between J and J in the two definitions (jS3|) and 
( jri] ) of the anisotropy parameter r. The first one was introduced as the electron 
tunneling matrix element, while the second one enters in the LD functional as 
the characteristic Josephson energy for the order parameter. Later on, in the 
framework of the microscopic theory, it will be demostrated that, in accordance 
with our qualitative definition, r ~ J 2 , while J turns out to be proportional to J 2 
too. In the dirty case it depends on the relaxation time of the electron scattering 
on impurities: J ~ a J 2 max{r, 1/T}. Hence both definitions (|63|), appearing 
in the qualitative consideration, and (jT^), following from the LD model, are 
consistent. 
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For a film of thickness d the integral over k z in Exp.(|73|) has to be replaced 
by a summation over the discrete k z and when £ z (T) 3> d only the term with 
k z — has to be taken into account: 

e 2 T £ 2 

X(fiim,±) = -^—- (76) 

Note that these formulas predict a nontrivial increase of diamagnetic suscep- 
tibility for clean metals p9| . The usual statement that fluctuations are most 
important in dirty superconductors with a short electronic mean free path 
does not hold in the particular case of susceptibility because here £ turns out 
to be in the numerator of the fluctuation correction. 

Now we will demonstrate that, besides the crossovers in its tempera- 
ture dependence, the fluctuation induced magnetization is a nonlinear func- 
tion of magnetic field too, and these nonlincarities, different for various 
dimensionalities, take place at relatively low fields. This, strong in com- 
parison with the expected scale of H C 2(0), manifestation of the nonlinear 
regime in fluctuation magnetization and hence, field dependent fluctuation 
susceptibility, was the subject of the intensive debates in early seventies 
[|Fj3^,||,||||||,||,|3j3|,|| (see also the old but excellent review of W.J. 
Skocpol and M. Tinkham gfj) and after the discovery of HTS pl ^^ Pl 



( see also very recent detailed essay of T.Mishonov and E.Penev M5| with 
references there). We will mainly follow here the paper of Buzdin et al. |Q, 
dealing with the fluctuation magnetization of a layered superconductor, which 
permits observing in a unique way all variety of the crossover phenomena in 
temperature and magnetic field. 

Let us go back to the general expression (|7^) and evaluate it without 
taking the magnetic field to small. The difficulty in dealing with it consists 
in the divergence of the sum over Landau levels n. This divergence can be 
regularized (see po] , p7[ ), but let us observe that in order to calculate the 
magnetization we must know the magnetic field dependent part of the free 
energy only. So a very convenient method to bypass the divergence problem 
jl8) is to calculate the difference F(H) — ^(0), turning the sum over Landau 
states in F(e, 0) into an integral and then, in its turn, turning this integral 
into a sum of integrals over the unit length intervals x S [n — 1/2, n + 1/2]. 
Then 



F(e,0 = - lim — T / — ^ V / dx x 

xln ttT 
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and by introducing the dimensionless variable [^] 

h = H c2 (0) = 2maT c /e = <P /2<l y , (77) 

one can write 

F(e,H)-F(e,0)= (78) 

TV , r j V- f 1/2 , , (2n+l + 2x)h + r/2(l-cosz) + e 
= —7; — T^-n I dzy I axm — — — . 

1™e xy J-* ^ 7-1/2 (2n + I)/i + r/2(I-cosz) + e 

Performing the integrations over z and x in (|7^) and differentiating with 
respect to h we finally obtain a very convenient general expression for the 
fluctuation magnetization in a layered superconductor: 



mi m T V^J 1 viRn + l) 



i- + ^+1/2 1 

V(i«n + l/2) V(^n + 1/2) 2 - P 2 J 



with (p(x) = x + -\/a; 2 — p 2 , i?„ = n + e/2/i + p and p = r/2h. The sum in ([79]) 
converges as 1/n 2 and it provides a volume magnetization expression that 
can be compared with experiment. 

Let us comment on the different crossovers in the M(e, H) field depen- 
dence analyzing the general formula (f79j) . Let us fix the temperature e <C r. In 
this case the c-axis coherence length exceeds the interlayer distance (£ 2 ^> s 
) and in the absence of a magnetic field the fluctuation Cooper pairs motion 
has a 3D character. Supposing the magnetic field to be not too high (h <C r) 
we may perform an expansion in ( [79] ) inl / p and obtain 

T\/2h 1/2 
M (3) (e < r,h< r) = —— — x 



*o6 

1/2 3 



e 

2li 



1/2 



1 2 1 2/i 

For weak fields (h <C e) the magnetization grows linearly with magnetic 
field, justifying our preliminary qualitative results. Nevertheless, this linear 

15 Let us remind that the exact definition of H C 2(0) contains the numerical coeffi- 
cient A(0) (see footnote 3). 
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growth is changed to the nonlinear 3-D high field regime M ~ \[H already 
in the region of a relatively small fields H C 2(e) < H (e < h). The further 
increase of magnetic field leads to the next crossover in the magnetization field 
dependence at h ~ r . However the Exp.(|S0|) was obtained in the assumption 
< r as an expansion over l/p, and it does not work any more. Handling 
with the Hurvitz zeta-function the summation in (|T9|) for 3-D case can be 
carried out for an arbitrary field J34[ : 



M (3 )(e<r, h) = 

\c(- 



T 



$ s 
1 1 e 

2' 2 + 2h 



1/2 



(81) 



-c 



1 1 

2' 2 



e 

2h 



e 

6h 



One can see from this formula that for large fields the magnetization saturates 
at the value M Xj lUiJl : 



T 



M(h > r) -> A/oo = = -0.346 , 

1 '' $ s 



(82) 



In 2 T 
~2~l^~s 

that is a typical for 2D superconductors. Therefore at h *~ r we have a 
3-D — ► 2-D crossover in M(H) behavior in spite of the fact that all sizes 
of fluctuation Cooper pair exceed considerably the lattice parameters. The 
effective " bidimensionalization" of the fluctuations is related to the effect of 
a strong magnetic field which "freezes out" the rotations along its direction. 
Let us stress that this crossover occurs in the region of already strongly non- 
linear dependence of M{H) and therefore for a rather strong magnetic field 
from the experimental point of view in HTS. 

Fixing the temperature e 3> r in the formula ([79]) one can find the general 
formula for 2D fluctuation regime [H : 



^>>^) = ^{inr(i + 4 



-ln(27r) 



e 

2h 



1> 



e 

2h 



(83) 



where r (x) is the Euler gamma- function and "0 (x) = dhi T (x) / dx is already 



cited in Introduction digamma- function. Using (83) one can directly pass 
from the linear regime in a weak magnetic field corresponding to (|6^) to the 
saturation of magnetization (82) in strong fields. 

Near the line of the upper critical field (/i C 2(e) = — e) the contribution of 
the term with n — in the sum ( |7^ ) becomes most important and for the 
magnetization the expression 
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Fig. 3. Schematic representation of the different regimes for fluctuation magneti- 
zation in the (H,T) diagram. The line H* 2 (T) is mirror-symmetric to the H C 2(T) 
line with respect to a t/-axis passing through T = T c . This line defines the crossover 
between linear and non-linear behavior of the fluctuation magnetization above T c . 



can be obtained. It contains the already familiar for us "OD" regime (r <C 
h — h C 2 <C 1), where the magnetization decreases as —M(h) ~ h _ i h (com- 
pare with (|67[)), while for h — h c2 « r the regime becomes "ID" and the 
magnetization decreases slower, as —M(h) ~ ■ 

Such an analogy is observed in the next orders in Gi too. In the Ref. 
the analogy was demonstrated for the example of the first eleven terms for 
the 2D case and nine for the 3D case. Summation of the series of high order 
fluctuation contributions to the heat capacity by the Pade-Borel method re- 
sulted in its temperature dependence similar to the 01? and II? cases without 
magnetic field. Nevertheless a considerable difference has not be forgotten: 
in the OD and ID cases no phase transition takes place while in the 2D and 
31? cases in a magnetic field a phase transition of first order to the Abrikosov 
vortex lattice state occurs. 

In conclusion, the fluctuation magnetization of a layered superconductor 
in the vicinity of the transition temperature turns out to be a complicated 
function of temperature and magnetic field, and it evidently cannot be fac- 
torized on these variables. The fit of the experimental data is very sensitive 
to the anisotropy parameter r and allows determination of the latter with 
a rather high precision f|l],[52| . In Fig. || the successful application of the 
described approach to fit the experimental data on YBa2Cu3C"7 is shown 

m. 



3 Ginzburg-Landau theory of fluctuations in transport 
phenomena 

The appearance of fluctuating Cooper pairs above T c leads to the opening 
of a "new channel" for charge transfer. In the Introduction the fluctuation 
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Fig. 4. Fluctuation magnetization of a YBaC0123 normalized on y~H as the 
function of temperature in accordance with the described theory shows the crossing 
of the isofield curves at T = T c (0) = 92.3K. The best fit obtained for anisotropy 
parameter 1—0.09. In the inset the magnetization curves as the function of magnetic 
field are reported. 



Cooper pairs were treated as carriers with charge 2e while their lifetime tql 
was chosen to play the role of the scattering time in the Drude formula. 
Such a qualitative consideration results in the Aslamazov-Larkin (AL) pair 
contribution to conductivity ( |ll|) (so-called paraconductivity^). 

Below we will present the generalization of the phenomenological GL func- 
tional approach to transport phenomena. Dealing with the fluctuation order 
parameter, it is possible to describe correctly the paraconductivity type fluc- 
tuation contributions to the normal resistance and magnetoconductivity, Hall 
effect, thermoelectric power and thermal conductivity at the edge of tran- 
sition. Unfortunately the indirect fluctuation contributions are beyond the 
possibilities of the description by Time-Dependent GL approach and they 
will be calculated in the framework of the microscopic theory (see Sections 
6-8). 



3.1 Time dependent GL equation 

In previous Sections we have demonstrated how the GL functional formalism 
allows one to account for fluctuation corrections to thermodynamical quan- 
tities. Let us discuss the effect of fluctuations on the transport properties of 
a superconductor above the critical temperature. 

In order to find the value of paraconductivity, some time-dependent gen- 
eralization of the GL equations is required. Indeed, the conductivity char- 
acterizes the response of the system to the applied electric field. It can be 
defined as E = —dA/dt but, in contrast to the previous Section, A has to 

16 This term may have different origins. First of all, evidently, paraconductivity is 
analogous to paramagnetism and means excess conductivity. Another possible 
origin is an incorrect onomatopoeic translation from the Russian "paroprovodi- 
most' " that means pair conductivity 
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be regarded as being time dependent. The general non-stationary BCS equa- 
tions are very complicated, even in the limit of slow time and space variations 
of the field and the order parameter. For our purposes it will be sufficient, 
following p6| , |67p^ , |6^ , [70 71,73], to write a model equation in the vicinity of 



T c , which in general correctly reflects the qualitative aspects of the order 
parameter dynamics and in some cases is exact. 

Let us keep in mind the GL functional formalism introduced above. If a 
deviation from equilibrium is assumed, then it is no more possible to derive 
the GL equations starting from the condition that the variational derivative 
of the free energy is zero. At the same time, in the absence of equilibrium 
\P begins to depend on time. For small deviations from the equilibrium it 
is natural to assume that in the process of order parameter relaxation its 
time derivative d'P/dt is proportional to the variational derivative of the free 
energy SJ-'/SW*, which is equal to zero at the equilibrium. But this is not 
all: side by side with the normal relaxation of the order parameter the effect 
of thermodynamical fluctuations on it has to be taken into account. This 
can be done by the introduction the Langevin forces £(r, t) in the right-hand 
side of equation describing the order parameter dynamics. Finally, gauge 
invariance requires that d^/dt should be included in the equation in the 
combination d\P/dt+ 2ze<ptf', where ip is the scalar potential of the electric 
field. By concluding all these speculations one can write the model time- 
dependent GL equation (TDGL) in the form 

-7GL(-g- t +2ieipU=— + C(r,t). (84) 



with the GL functional T determined by (|16|),(|65|),(|69|). The dimensionless 
coefficient ^gl i n the left-hand-side of the equation can be related to pair 
life-time tgl (0): Igl — cxT c €Tgl — irot/8 by the substitution in ( pi[ ) of the 
first term of ([!(]) onlyPj 

Neglecting the fourth order term in the GL functional, equation ( pi] ) can 
be rewritten in operator form as 

[L- 1 - 2ie lG Mr, t)]#(r, t) = C(r, t) (85) 

with the TDGL operator L and Hamiltonian TL defined as 



L 



8 J.G 



H = aT c 



£-?(V-2ieA) 2 . (86) 



17 It will be shown below that taking into account the electron-hole asymmetry 
leads to the appearance of an imaginary part of ~/gl proportional to the deriva- 
tive d\n(pv 2 T) /8E\e f ~ 0(1/Ep). This is important for such phenomena as 
fluctuation Hall effect or fluctuation thermopower and, having in mind writing 
the most general formula, we will suppose jgl = Tra/S+ilmjcL, where necessary. 
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We have introduced here the formal operator of the coherence length £ to 
have the possibility to deal with an arbitrary type of spectrum. For example, 
in the most interesting case for our applications to layered superconductors, 
the action of this operator is defined by the Exp.(|69|). 

In the absence of an electric field one can write the formal solution of 
equation (|85|) as 



<f (0 >(r,t) = L£(r,t). (87) 

The Langevin forces introduced above must satisfy the fluctuation-dissipation 
theorem, which means that the correlator < Wp^* (t')Wp\t) > at coinciding 
moments of time has to be the same as < Itf'pl 2 >, obtained by averaging 
over fluctuations in thermal equilibrium. This requirement is fulfilled if the 
Langevin forces £(r, i) and C*(r, t) are correlated by the Gaussian white-noise 
law 

<TM)C(r',<') >=2TRe lGL S(r-v')S(t-t). (88) 

To show it let us restrict ourselves for sake of simplicity to the case of A = 
and calculate the correlator 

< &*(r,t)&(r',t) >=< C*M)L*£C(r',i) >= (89) 

. m^/ ^<«~'> £ f i-ta IW» 0). 

L(p, D) can be found by making the Fourier transform in 

L{p,n) = {-iiglV + Ep)-\ (90) 
resulting in 

t°° dfi 1 

{W*(r,tp(r',t)) p = 2TB^ GL = (|^ p | 2 ) , (91) 

where 

£ P -«T c (e + cV). (92) 
is the Cooper pair energy spectrum. 



3.2 Paraconductivity 

Let us try to clear up the reason why the simple Drude formula works so 
well for complex phenomenon like paraconductivity. For this purpose let us 
try to derive the Boltzmann master equation for the fluctuation Cooper pair 
distribution function 

«p(t)= / <^(r,t)^(r',i))exp(- ip (r-r')Kr-r'). (93) 
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Let us recall that in the state of thermal equilibrium n p ^ =< |^ p | 2 >= T/e p . 

In order to determine the electric field dependence of n p let us write its 
time derivative using ( |S4| ) 



dn p (t) 
dt 



< 



dt 



V*(r',t) > + < <P(r',t) 



dW*(r,t) 
dt 



d(r 



-ip(r— r ) 



+ J_ Re U^ t ) S -^(r',t)) + —Re(ar,t)W*(r',t)) 

TGL \ °* / IGL 



where is determined by (pq). Expressing the scalar potential by the electric 



field E one can transform the first term of the last integral into — 2eE^ £ . 
The term with the variational derivative can be evaluated by means of (^6|) 
and expressed in the form — ^^-e p n p . More cumbersome is the evaluation of 
the last term, containing the Langevin force. To the first approximation it is 
possible to use here the order parameter )p( \r,t) (see (jg7|)) unperturbed by 
the electric field as the convolution L((r, t). In this way, using (|8|) and ( |B9"| ) 
we calculate the last average in (p3) 



> 



2-|Wr) - <p(r')] < &(r,m*(r',t) > + (94) 
i 



— — — [ d(r - r')e-' lp< - r - r '^Re(e(r,t)L((r,t) 

IGL J \ 



(0) 



c 7) 

P*) D ~ —y-oc 2^-^'"' 7GL 7GL P P 

and obtain the transport equation 



<9t dp 7gl 



(np-n(°))=-^(n p -<)). (95) 



In the absence of a magnetic field e p was determined by (£^) and the mo- 
mentum dependent lifetime, corresponding to the Ginzburg-Landau one, can 
be introduced: 

, r GL (e) 

rp=7Gi/£p= mw 

Let us stress the appearance of the coefficient 2 on the right hand side of the 
equation (|9^). This means that the real Cooper pair lifetime, characterizing 
its density decay, is r p /2 . 

The effect of a weak electric field on the fluctuation Cooper pair distribu- 
tion function in the linear approximation is determined by 

eE7GL dn ( p 0) eTj GL de p 



dp e p op 
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Substituting this formula into the expression for the electric current side by 
side with the Cooper pair velocity v p = de p /dp one can find 



p 



(2ev a n v )=a afi E fi , (97) 
where the paraconductivity tensor components are: 

a/3 T 2 T V V p^ 

(D) = 4 e aT l^—„ 



^ = ~^V^. (98) 



p P 

In the case of isotropic spectrum 



ad 71 1 



" V 

= < 



e 2 



p 



'P 



32^ 

me 

Tre 2 ^ 1 
I 165 e 3 / 2 



3D case, 

2D film, thickness :d -C £, 

ID wire, cross — section :5 -C £ 2 . 



One can compare this result with that carried out in the Introduction 
from qualitative consideration, based on the Drude formula. Those simple 
speculations reflect correctly the physics of the phenomenon but were carried 
out with the assumption of the momentum independence of the relaxation 
time t p , taken as To = gpycrj \ ■ As we have just seen, in reality r p decreases 
rapidly with increase of the momentum, the excess "2" appeared in ( p5| ) 
because of the wave nature of fluctuation Cooper pairs; accounting for this 
circumstance results in the precise coefficients of (p9|), different from (p"T|). 

An especially simple form the paraconductivity results in the 2D case, 
where, calculated per unit square, it depends on the reduced temperature 
only: 

^ = I^i' (1 ° 0) 

The coefficient in this formula turns out to be a universal constant and is given 
by the value h/e 2 = A.lkH. For electronic spectra of other dimensionalities 
this universality is lost, and the paraconductivity comes to depend on the 
electron mean free path. 

Let us compare <tq (T) with the normal electron Drude part a n = n e e 2 r /m 
by writing the total conductivity 

n 2nn loe 

One sees that at e cr — 0A/(ppl) ~ Gin c ) the fluctuation correction reaches 
the value of the normal conductivity. Let us recall that the same order of 
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magnitude for the 2D Ginzburg-Levanyuk number was obtained above from 
the heat capacity study. We will discuss the region of applicability of (101) 
below in Section 8.6. 

It is worth mentioning that the results derived here for paraconductivity 
are valid with the assumption of weak fluctuations: for the temperature range 
e < Gi(Dj they are not anymore applicable. Nevertheless, one can see that for 
not very dirty films, with p 2 F ld 3> 1, a wide region of temperatures Gir^d) ^ 
e <C 1 exists where the temperature dependence of conductivity is determined 
by fluctuations and in this region the localization effects are negligible. 

The transport equation ( 1951) w as originally derived many years ago by 
L.G.Aslamazov and A. I. Larkin [|73|. Recently T.Mishonov et al. |7J] rederived 
Eq.(95) and solved it for n p in the case of an arbitrary electric field. 



3.3 General expression for paraconductivity 

Unfortunately the applicability of the master equation derived is restricted 
to weak magnetic fields (H H C 2{e)). For stronger fields H C 2{e) < H <C 
H c2 (0) the simple evaluation of averages in (^4|) turns out to be incorrect, 
the density matrix has to be introduced and the master equation loses its 
attractive simplicity. At the same time, as we already know, namely these 
fields, quantizing the fluctuation Cooper pair motion, present special interest. 
This is why in order to include in the scheme the magnetic field and frequency 
dependencies of the paraconductivity, we return to the analysis of the general 
TDGL equation (|84|) without the objective to reduce it to a Boltzmann type 
transport equation. 

Let us solve it in the case when the applied electric field can be consid- 
ered as a perturbation. The method will much resemble an exercise from a 
course of quantum mechanics. To carry out the necessary generality side by 
side with a formal simplicity of expressions we will introduce some kind of 
subscript {i} which includes the complete set of quantum numbers and time. 
By a repeated subscript a summation over a discrete and integration over 
continuous variables (time in particular) will be supposed. 

We will look for the response of the order parameter to a weak electric 
field applied in the form 

!^,*) (102) 

where So? is determined by (]§7|). Substituting this expression into (|8^) and 
restricting our consideration to linear terms in the electric field we can write 

(L-% k} W$ } = 2ie lGLna} ^ 

with the solution in the form 

&{{} = 2ie lGLL{ik}tp{ki}L{i m yC{ m }- (103) 
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Let us substitute the order parameter (102) in the quantum mechanical ex- 
pression for current 



(o) 



(104) 



where V{ ik y is the velocity operator which can be expressed by means of the 
commutator of r with Hamiltonian (|S6|): 



V{ifc} = i{H,r} {ik} . 



(105) 



The second term of (104) can be written by means of a transposed velocity 
operator (which is Hermitian) as the complex conjugated value of the first 
one: 



which results in 



(106) 
(107) 



Let us average now ( |107[ ) over the Langevin forces moving the operator L*^ kj ^ 
from the beginning to the end of the trace and using (|88|) . One finds 



8e 2 Im{7 GL L| fei jV {i;} L {iTO} 93 { „ m} L { „ p }<'{fe}C{p}}- 



- 16Te 2 Re(7Gi)Im{7 GL v {il} L {im} ip {mn} L {np} L* {pi} } 



(108) 



Now we choose the representation where the £{/ m } operator is diagonal (it 
is evidently given by the eigenfunctions of the Hamiltonian (pfl) ): 



£{m} - i^lGL ' 



(109) 



where £{ m } are the appropriate energy eigenvalues. Then we assume that 
the electric field is coordinate independent but is a monochromatic periodic 
function of time: 



ip(r,t) = -E^r 13 exp(-iojt). 



(110) 



In doing the Fourier transform in (108) one has to remember that the time 
dependence of the matrix elements <P{ mn } results in a shift of the frequency 
variable of integration fl — > i? — uj in both L-operators placed after <fs mn \ 
or, what is the same, to a shift of the argument of the previous L{/ m } for uj : 



jZ = 16Te 2 Re( 7GL ) 



(111) 



_K{ 7G L^ } L {/} (^ + w )[-^ i} ; 
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where »/(w) = [/(w) + /*(-w)]/2. 

Let us express the matrix element r^j by means of ^{u} using the com- 



mutation relation (105). One can see that in the representation chosen 

tf,.,=i V{H} (112) 



and, carrying out the frequency integration in (Jlllj), finally write for the 
fluctuation conductivity tensor (j^ = a a/3 (uj)E^) : 



a a P{e,H,Lu) = 8e 2 TRe{ lGL ) ^ 

{M}=0 



(113) 



5? 



1GL 



r {ti} W {H} 



£{i}(lGLe {i} + l* GL £ {i} ~ AlGL?u){e { i} - e {i} ) 



This is the most general expression which describes the d.c, galvanomagnetic 
and high frequency paraconductivity contribution. In the case when we are 
interested in diagonal effects only, where it is enough to accept ^gl as real 
ilGL = Re7c£ = na/8) omitting its small imaginary part, the last expression 
can be simplified by means of symmetrization of the summation variables): 



oo 



*{e,H,uj) = -ae 2 T ^ 
{i,i}=a 



V" V 



£{i}£{i}(£{i} + £{1} - ilGLu) 



(114) 



Let us demonstrate the calculation of the d.c. paraconductivity in the 
simplest case of a metal with an isotropic spectrum. In this case we choose a 
plane wave representation. By using e p defined by ( |92"| ) one has 



de 

9 1pp'} = v pV' v p = "ST = 2 «^ 2 p. 



dp 



(115) 



We do not need to keep the imaginary part of Jgl, which is necessary to 
calculate particle-hole asym metric effects only. Then the fluctuation conduc- 
tivity calculated from (114) coincides exactly with (|9^). 



3.4 Fluctuation conductivity of layered superconductor 



Let us return to the discussion of our general formula ( |114|) for the fluctuation 
conductivity tensor. A magnetic field directed along the c-axis still permits 
separation of variables even in the case of a layered superconductor. The 
Hamiltonian in this case can be written as in (IT1 



H = aT c (e - & y (V xy - 2ieA xy 



cos(A: z s)^ 



(116) 
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It it is convenient to work in the Landau representation, where the summation 
over {i} is reduced to one over the ladder of the Landau levels i = 0, 1,2.. 
(each is degenerate with a density 2eH per unit square) and integration over 
the c-axis momentum in the limits of the Brillouine zone. The eigenvalues of 



the Hamiltonian (116) can be written in the form 



£„ = aT c [e + - cos{k z s)) + h(2n + 1)] = e kz + aT c h{2n + 1), (117) 

where h = 2 maT was already defined by Eq. ([77|) . For the velocity operators 
one can write 

v x ' y = _HV-2zeA)^; v z = ~^T C sin(fc zS ). (118) 

In-plane conductivity. Let us start from the calculation of the in-plane 
components. The calculation of the velocity operator matrix elements requires 
some special consideration. First of all let us stress that the required matrix 
elements have to be calculated for the eigenstates of a quantum oscillator 
whose motion is equivalent to the motion of a charged particle in a magnetic 
field. The commutation relation for the oscillator's velocity components is 
well known (see f7j|): 

[^]=i^ = —h. (119) 
2m z to 

In order to calculate the necessary matrix elements let us present the velocity 
operator components in the form of boson-type creation and annihilation 
operators a + , a with commutation relation [a,a + ] = 1 : 



-x.y _ otT c h I a++a 
V 2m \ ia + — ia 

One can check that the correct commutation relation ( |119| ) is fulfilled. Taking 
into account that 

< l\a\n >=< 77|a + |Z >= y/n5 nj i+i 

it is seen that the only non-zero matrix elements of the velocity operator are 



< l\v x ' v \n >= ■ ' aTcfl ( ^ l - n+1 + V^n,l+l 



2m \iy/l5i,n+i — Wn&n,l+\ 

Using these relations the necessary product of matrix elements can be calcu- 
lated: 

< l\v x \n >< n\v x \l >= ^(lS l>n+1 + nS n , l+1 ). (120) 

2777 
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Its substitution into (114) and accounting of the degeneracy of the Landau 
levels 2eH = AmaT c h gives for the diagonal component of the in-plane con- 
ductivity tensor 



ira 2 T 2 e 2 ^ -y^ ^. (^z,n+i + nS n j + i) 

4m !lo inn £ l £ n{ £ l + £n ~ ijGL^) 

{nj— {/}— 



(121) 



7re^ (aT c )° h z I - — -] 



Expanding the denominator into simple fractions we reduce the problem to 
the calculation of the c-axis momentum integral, which can be carried out in 
the general case by use of the identity: 



dx 



1 



1 



2n cos x 



1 



(122) 



valid for any complex parameter z ^ 1 with the proper choice of the square 
root branch. Using it we write the general expression for the in-plane com- 
ponent of the fluctuation conductivity tensor 



p 2 h ^ I 
o™( E ,h,cj) = —T,(n+l) 



1 



16s 
1 



h-iuj y/[ e + h(2n + 1)] [r + e + h(2n + 1 
1 



h + iw y/[ e + h(2n + 3)] [r + e + h{2n + 3)] 

2h 1 
h 2 + uj 2 y/[ e + h(2n + 2) - iu] [r + e + h(2n + 2) - iu] 



(123) 



where u> — jfjp. 



Out-of plane conductivity. The situation with the out-of plane component 
of paraconductivity turns out to be even simpler because of the diagonal 



structure of the v 



{in} 



- s ^T c sia(k z s) x S in x 8{k z — k z >). Taking into 



account that the Landau state degeneracy 2eH = h^ 2 we write 



1 oo 

a zz (e,H,uj) = -Tzae 2 T ^ 3t 



{i,l}=0 

ire 2 (aT c ) 3 ( sr 

16 



'{il} V {h} 



£ {i} £ {i}{ £ {i} + £ {i} -iaT c w) 



dk, 



sin 2 (A; z s) 



_2i 2-k e 2 l (k z )[e n {k z ) - iaT c uj] 
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The following transformations are similar to the calculation of the in-plane 
component: we expand the integrand into simple fractions and perform the 
fc 2 — integration by means of the identity (122). The final expression can be 
written as 



9 / \ 2 oo 

7re / sr 



32s V £x 



5? 



1 



n=0 



OX 



(124) 



1 



X + iujJ [ y/(e + h(2n + 1) + A)(e + h(2n + 1) + A + r) 
1 1 , 



y/(e + h(2n + 1) - iu)(e + h(2n +1) -iuj + r) 



\x=o 



Analysis of the general expressions. In principle the expressions derived 
above give an exact solution for the a.c.(w <C T) paraconductivity tensor of a 
layered superconductor in a perpendicular magnetic field H <C H C 2 (h <C 1) 
in the vicinity of the critical temperature (e -C 1). The interplay of the 
parameters r,e,ui h entering into (123)-(124), as we have seen in the example 
of fluctuation magnetization, yields a variety of crossover phenomena. 

1. The simplest and most important results which can be derived are the 
components of the d.c. paraconductivity (ui — 0) of layered superconductor 
in the absence of magnetic field. Keeping oj = and setting h — > one can 
change the summations over Landau levels into integration and find 



'(£,/!-> 0,LU = 0) 



! (e,/i^ 0,w = 0) 



1 



16* v/R 



•r/2 



32^ V [e(e + r)]V2 



2. The Aslamazov-Larkin contribution to the magnetoconductivity can be 
studied by putting u) — and keeping magnetic field as ar bitrary. We will not 
go into details and just report the results (following [ 1 54 1 with some revision 
of the coefficient in 3D case) 



h < e 



e < /i< r 
(3D) 



max{e, r} <C h 
(2D) 



4s VSgr 



£11 

8s h 



: (e,h = Q) 



e 2 s r z {t+r/2) ■ 



3.24e z s 



7C(3)e^ r 2 
2' J Z 2 W 



Table 2. 
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Here it is worth making an important comment. The proportionality of the 
fluctuation magnetoconductivity to h 2 is valid when using the parametriza- 
tion e = (T — T c q)/T c q only. As it well known, a weak field shifts the criti- 
cal temperature linearly, which often makes it attractive to analyze the ex- 
perimental data by choosing as the reduced temperature parameter e/, = 
(T — T C (H))/T C (H). In this parametrization one can get a term in the mag- 
netoconductivity linear in h , which previously was cancelled out by the 
magnetic field renormalization of the critical temperature. So it is important 
to recognize that the effect of a weak magnetic field on the fluctuation con- 
ductivity cannot be reduced to a simple replacement of T c q by T C (H) in the 
appropriate formula without the field. Vice versa, this effect is exactly com- 
pensated by the change in the functional dependence of the paraconductivity 
in magnetic field, and finally it contains the negative quadratic contribution 
only. 

3. Letting the magnetic field go to zero and considering nonzero frequency 
of the electromagnetic field one can find general expressions for the compo- 
nents of the a.c paraconductivity tensor. They are cumbersome enough and 
in the complete form can be found, for instance, in fr(| . We recall here the 
simplified asymptotics for the Recr in the 2D regime only: 



16s e 



2e uj 

— arctan 

uj e 



^) In 

uj 



(125) 















1 + 


(r 











e 

2^ 



ST 
£,xy 



ln[l- 



(126) 



The general formulas ( |123[ )-(124) allow one to study the different crossovers 
in the a.c. conductivity of layered superconductor in the presence of magnetic 
field of various intensity. We leave this exercise for the reader having some 
practical interest in the problem. 



3.5 Magnetic field angular dependence of paraconductivity 

We have seen above that in the case of a geometry with a magnetic field 
directed along the c-axis many sophisticated fluctuation features of layered 
superconductors can be studied in the most general form. Nevertheless even 
the attempt to explore the d.c. conductivity in a longitudinal magnetic field 
(directed in ab plane) |77| or, moreover, with the field directed at some arbi- 
trary angle 9 with the c-axis leads to the appearance of the a vector potential 
component in the argument of cos(k z s) and the problem requires a nontrivial 
calculation of the matrix elements over the Mathieu functions. 

We already learned that at temperatures very near to the critical one (e <C 
r) the 3D fluctuation regime takes place. Here the size of the Cooper pairs 
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along the c-axis is so large that the peculiarities of the layered structure do 
not play any more role. This means that only small values of k z are important 
in the ^-integrations, where the cos(fc z s) in ( |7l"| ) can be expanded and the LD 
functional is reduced to its traditional GL form with an anisotropic effective 
mass tensor: 



= j d 3 rl a|<Pf 



B, 



m 



i 



3 

y 



1 d 

i dx,. 



2eA u V 



+ 



8^ 



H B 

47T 



(127) 



We will demonstrate below that in this case a scaling approach provides a 
direct access to the most general results by rescaling the anisotropic problem 
to the corresponding isotropic one on the initial level of the GL approach 

s 

Let us suppose that the external field H is chosen to lie in the y — z 
plane and makes angle 6 with the z— axis. For sake of simplicity and because 
the oxide superconductors are within high accuracy uniaxial materials, we 
choose ra x = m y — m* , while m^ 1 = 2as 2 r (compare with (|S9|)). The 
effective anisotropy parameter 7^ = m* /m z = 2as 2 rm* < 1 is introduced. In 



(127) the anisotropy enters only in the gauge-invariant gradient term, so the 
simple rescaling of the coordinate axes: x = x,y = y, z = j a z together with 
the scaling of the vector potential: A = (A x , A y , Az/ja) will render this term 
isotropic. The magnetic field evidently is rescaled to B = (B x /j a , B y /j a , B z ) 
and the last three terms in (127), describing the magnetic-field energy, are 
transformed to 



+ {-SL + Bl)- 

7a 



In short, we have removed the anisotropy from the gradient term but rein- 
troduced it into the magnetic energy term. In general it is not possible to 
make both terms isotropic in the Gibbs energy simultaneously. However, de- 
pending on the physical question addressed, we can neglect fluctuations in 
the magnetic field, as was mostly done above. 

Let us demonstrate how the method works for the example of the d.c. 
fluctuation conductivity tensor which was calculated above for a magnetic 
field directed along the c-axis. We restrict our consideration to the 3D region 
(e <C r). One can write the scaling relations between the electric field and 
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current components before and after the scaling transformation by means of 
a conductivity tensor and the anisotropy parameter: 



jx.y = jx,y jz ~ GV Z ~ J a j z (128) 

dw 1 ~ 

F — F F ~ — — ~ F, 

dz 7q 

Now let us rewrite the relations between the current and electric field vectors 
before and after the scale transformation 



3 a — VapEp 

ja = VupEp. 



Comparing them with (128) and introducing the operator of the direct scaling 
transformation T a p 



T, 




a/3 



one can write j a — T atl jp, E a = (T 1 ) Q/1 i? AI and express the conductivity 
tensor as 

c Q /3 = Ta^a^pTpp. 

Now let us work in the already isotropic coordinate frame. We suppose 
that initially the magnetic field was directed along the c-axis and now we 
rotate it in the X-Z plane by the angle 9 with respect to the initial direction. 
The conductivity tensor will be transformed by the usual matrix law: 

a a p(d) = R T a ^pp{0)R pP = Rl^T-^p^mT-^sRsp 

and 

cr a /3{9) = T ai a 1 s{9)Ts ! s = T ai R^ fJi (T~ 1 ) t _ ti: a <iri (0)(T~ 1 ) v sRsKT K f3, 
where 

(cos 9 0— sin 9 v 
sin# cos^ 

Finally the fluctuation conductivity tensor cr a p(9) in the initial tetragonal 
system with the magnetic field directed at the angle 9 with respect to the 
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c-axis can be expressed by means of the effective transformation operator 

M af3 : 



a af3 (9) = M^(9)a, v (0,H)M^{9), 
with 

M aP {9) = {T-^sRskTkp = 

The angle 9 can be expressed via the renormalized magnitude of the magnetic 
field H = VW+WI ■ 

~ H z cos 9 ~ 7 a sin 9 

" — ; sm.9 — 




H J cos 2 9 + 7 2 sin 2 9 J cos 2 9 + 7 2 sin 2 9 

In the case of the paraconductivity of a layered superconductor with the 
magnetic field applied at an arbitrary angle 9 the answer can be written in 
the general form by means of the three diagonal components of conductivity 
CTjj(Q, H) in the perpendicular field H: 



Vaf3{9) 



1 

cos 2 9 + 7 2 sin 2 9 

/ a xx cos 9+ 
+lt a zz sin 2 I 






o yy (cos 2 9 + 7 2 sin 2 i 



V°zz7a sin # cos 
In the simplest case of a longitudinal field 9 = 90° : 



cr 2 2 7 a sin 2 9 cos 0- \ 
-a xx sin 9 


cr. 2 COS 2 9 



' 7 2 a Z2 (0, 7 ai?) 







a y y 










(129) 



4 Fluctuations near superconductor-insulator 
transition 



4.1 Quantum phase transition 

It is usually supposed that the temperature of the superconducting transi- 
tion does not depend on the concentration of non-magnetic impurities (An- 
derson's theorem J79lp0]]). Nevertheless when the degree of disorder is very 
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high Anderson localization takes place, and it would be difficult to expect 
that under conditions of strong electron localization superconductivity can 
exist, even if there is inter-electron attraction. This means that at T = the 
phase transition takes place with a change of the disorder strength or carrier 
concentration. Such a transition is called a quantum phase transition since 
at zero temperature the classical fluctuations are absent. Indeed, one can see 
from ([7]) that in the limit T — > the thermal fluctuation Cooper pairs vanish. 

In the metallic phase of a disordered system the conductivity is mostly de- 
termined by the weakly decaying fermionic excitations, their dynamics yield- 
ing the familiar Drude formula (the method which accounts for the fermionic 
excitations will be referred to as the Fermi approach later on). Inside the 
critical region the charge transfer due to fluctuation Cooper pairs turns out 
to be more important. In some approximation, the pairs may be considered 
as Bosc particles. Therefore the approach dealing with the fluctuation pairs 
will be called below the Bose approach. 

Let us suppose that at temperature T — the superconducting state 
occurs in a weakly disordered system. In principle two scenarios of the devel- 
opment of the situation are possible with an increase of the disorder strength: 
the system at some critical disorder strength can go from the superconduct- 
ing state to the metallic state or to the insulating state. The first scenario 
is natural and takes place in the following cases: if the effective constant of 
the inter-electron interaction changes its sign with the growth of the disor- 
der; if the effective concentration of magnetic impurities increases together 
with the disorder growth; if the pairing symmetry of superconducting state 
is nontrivial it can be destroyed even by the weak disorder level. We will 
study here the second scenario where the superconductor becomes an insu- 
lator with disorder increase. This means that at some disorder degree range, 
higher than the localization edge when the normal phase does not exist any 
more at finite temperatures, superconductivity can still survive. From the 
first glance this statement seems strange: what does superconductivity mean 
if the electrons are already localized? And if it really can take place beyond 
the metallic phase, at what value of disorder strength and in which way does 
the superconductivity finally disappears? 

One has to have in mind that localization is a quantum phenomenon in its 
nature and with the approach to the localization edge the coherence length 
of localization t grows. From the insulator side of the transition vicinity 
this means the existence of large scale regions where delocalized electrons 
exist. If the energy level spacing in such regions does not exceed the value 
of superconducting gap Cooper pairs still can be formed by the delocalized 
electrons of this region. 

The problem can be reformulated in other, already familiar, way: how does 
the critical temperature of the superconducting transition decrease with the 
increase of the disorder strength? In the previous Sections we have already 
tried to solve it by discussing the critical temperature fluctuation shift. We 
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have seen that the fluctuation shift of the critical temperature is proportional 
to y/Gi(z) for a 3D superconductor and to Gi( 2 ) ln{l/Gi( 2 )} for 2D. This 
means that the critical temperature is not changed noticeably as long as the 
Ginzburg-Levanyuk number remains small. So one can expect the complete 
suppression of superconductivity when Gi ~ 1 only. For further consideration 
it is convenient to separate the 3D and 2D cases because the physical pictures 
of the superconductor-insulator transition for them are quite different. 



4.2 3D superconductors 

As one can see from Table 1 in the 3D case the Ginzburg-Levanyuk number 
remains small at ppl ~ 1 : Gi^d) ~ (e^) ^ 1- Nevertheless, approaching 
the edge of localization, the width of the fluctuation region increases [^IJ. In 
the framework of the self-consistent theory of localization |82 such growth 



of the width of the fluctuation region was found in paper j83| . 

Instead of the cited self-consistent theory let us make some more general 
assumptions concerning the character of the metal-insulator (M-I) transition 
in the absence of superconductivity [jl3j . We suppose that in the case of very 
strong disorder and not very strong Coulomb interaction the M-I transition 
is of second order. The role of "temperature" for this transition is played 
by the " disorder strength" which is characterized by the dimensionless value 
g = ^t^- With its decrease the conductivity of the metallic phase decreases 
and at some critical value g c tends to zero as 

a = e 2 p F (g-g c )". (130) 

This is the critical point of the Anderson (M-I) transition. We assume that 
the thermodynamic density of states remains constant at the transition point. 

The electron motion in metallic phase far enough from the M-I transition 
has a diffusion character and the conductivity can be related to the diffusion 
coefficient T> = ppl/im by the Einstein relation: a = ve 2 T>. One can say that 
diffusion like "excitations" with the spectrum oj(q) = iDq 2 propagate in the 
system. At the point of the M-I transition normal diffusion terminates and 
conductivity, together with T>, turns zero. In accordance with scaling ideas, 
the diffusion coefficient can be assumed here to be a power function of q: 
T)(q) r~j q z ~ 2 , with the dynamical critical exponent z > 2. The anomalous 
diffusion excitation spectrum in this case would take the form lj ~ q z . 

In the insulating phase (g < g c ) some local, anomalous diffusion, confined 
to regions of the scale £, is still possible. It cannot provide charge transfer 
through out all the system, so V{q = 0) = 0, but for small distances (q > ^~ 1 ) 
anomalous diffusion takes place. Analogously, in the metallic phase (g > g c ) 
the diffusion coefficient in the vicinity of the transition has an anomalous 
dependence on q for q > t and weakly depends on it for q < £ . So one 
can conclude that the diffusion coefficient for q > l^ 1 from both sides of the 
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transition has the same g-dependence as for all q in the transition point. It 
can be written in the form 



PFt 



x, x 3> 1 

ip(x) = { 1, x< 1, 5 > g c , (131) 
0, x<^l,g <g c 



where the dimensionless localization length £, characterizing the spatial scale 
near the transition, grows with the approach to the transition point like 

e(9) = —(9-9c)-^- (132) 

PF 

The critical exponent in this formula is found from the Einstein relation in 
the vicinity of the M-I transition. 

At finite temperatures, instead of the critical point <? c ,a crossover from 
metallic to insulating behavior of a(g) takes place. The width of the crossover 
region is g—g c , where g is determined from the relation T>£~ 2 (g) ~ Ep[pF£(g)]^- 



T (we have used the second asymptotic of (131)). In this region the diffusion 
coefficient is 

Pf\T 



V(T) ~ Tt ~ — -=r (133) 



and it depends weakly on the g — g c . Beyond this region the picture of the 
transition remains the same as at T = 0. 

Let us consider now what happens to superconductivity in the vicinity of 
the localization transition. In the mean field approximation (BCS) the ther- 
modynamic properties of a superconductor do not depend on the character 
of the diffusion of excitations. This should be contrasted with the fluctua- 
tion theory, where such a dependence clearly exists. We will show that the 
type of superconducting transition depends on the dynamical exponent z. If 
z > 3, the transition to superconductivity occurs on the metallic side of the 
localization transition (we will refer to such a transition as S-N transition). 
If z < 3, the transition to superconductivity occurs from the insulating state 
directly (S-I transition). 

Let us study how the superconducting fluctuations affect the transition 
under discussion. In spirit of the GL approach fluctuation phenomena in 
the vicinity of the transition can be described in the framework of the GL 
functional (|2^). The coherence length in the metallic region, far enough 
from the Anderson transition, was reported in Introduction to be equal 
£ 2 = £J = QA2V/T . In the vicinity of the M-I transition we still believe 
in the diffusive character of the electron motion resulting in the pair forma- 
tion. The only difference from the previous consideration is the anomalous 
character of the quasiparticle diffusion. So in order to describe the supercon- 
ducting fluctuations simultaneously near superconducting (in temperature) 
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and Anderson (in g) transitions let us use the GL functional (f26J) with the 
fc-dependent diffusion coefficient ( |131| ). 

The value of Gi can be estimated from the expression for the fluctuation 
contribution to heat capacity ( pfj| ) taken at e ~ Gi, where the fluctuation 
correction reaches the value of the heat capacity jump: 

^ (134) 

vj {TGi + V(q)q 2 ) 2 

with T ~ T c . Let us approach the M-I transition from the metallic side. If 



we are far enough from transition, Gi is small and the integral in (134) is 
determined by the region of small momenta T>(q)q 2 < TGi : 

Two scenarios are possible: Gi becomes of the order of 1 in the metallic 
phase, or it remains small up to the crossover region, where finally reaches 
its saturation value. In the first case we can use the second asymptotic of 



( |131D for D(q) and find: 

Gi ~ m^) 3 *" 6 . (136) 

tip 

l 

One can see that Gi becomes of the order of 1 at pf£m ~ {Ep /T) 3z ~ e . 
Comparing this value with pp£(g) ~ (Ep/T) 1 ^ at the limit of the crossover 
region we see that for z > 3 the first scenario is realized. Concluding the first 
scenario discussion we see that the superconducting critical temperature goes 
to zero at £ = £m, still in the metallic phase, so at T = a superconductor- 
normal phase (S-N) type quantum phase transition takes place. 

The second scenario takes place for z < 3 when Gi remains small even at 
the edge of crossover region, reaching there the value 

, 2(3-*) 



In the crossover region the diffusion coefficient, and hence Gi, almost do 
not vary. This is why the temperature of superconducting transition remains 
almost frozen with further increase of disorder driving the system through the 
Anderson transition. The abrupt growth of Gi and decrease of T c take place 
when the system finally goes from the crossover to the insulating region. In 



the insulator phase the diffusion coefficient D(q < I x ) = and from (134) 
one can find for Gi : 



(138) 
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Comparing this result with the Table 1 it is easy to see that it coincides 
with the Ginzburg-Levanyuk number for a zero-dimensional granule of size 
I -C £(T). Hence we see that in the second scenario the Ginzburg number 

reaches 1 and, respectively, T c — > at pf^i ~ ( ^ ) , which is far enough 
from the M-I transition point. This is why in this case one can speak about 
the realization at T — of a superconductor-insulator (S-I) type quantum 
phase transition. The scale £i determines the size of the "conducting" do- 
mains in the insulating phase, where the level spacing reaches the order of 
the superconducting gap. It is evident that in the domain of scale £ < tj 
superconductivity cannot be realized. 

In the vicinity of a quantum phase transition one can expect the ap- 
pearance of non monotonic dependencies of the resistance on temperature 
and magnetic field. Indeed, starting from the zero resistance superconduct- 
ing phase and increasing temperature from T = 0, the system passes through 
the localization region, where the resistance is high, to high temperatures 
where some hopping charge transfer will decrease the resistance again. The 
analogous speculations are applicable to the magnetic field effect: first the 
magnetic field "kills" superconductivity and increases the resistance, then it 
destroys localization and decreases it. 

The phase diagram in the (T, g) plane has the form sketched in Fig. ^|. 
For gi = g c - (T c / E F ) 3{z ~ 2)/>c , an S-I transition takes place at T = 0. In- 
creasing the temperature from T = in the region < g < gi we remain in 
the insulating phase with exponential dependence of resistance on tempera- 
ture. For gi < g < <?_ at low temperatures < T < T c (g) the system stays 
in the superconducting state which goes to the insulating phase at higher 
temperatures. In the vicinity of the Anderson transition (gL < 5 < <7+) the 
superconducting state goes with growth of the temperature to some crossover 
metal-insulator state which is characterized by a power decrease of the resis- 
tivity with the increase of temperature. Finally at g c < g the superconducting 
phase becomes of the BCS type and at T = T c it goes to a metallic phase. 

The phase diagram in the magnetic field - disorder plane is similar to 
that in the (T, g) plane with the only difference that at T = there is no 
crossover region, instead a phase transition takes place. 



4.3 2D superconductors 

Preliminaries. As was demonstrated in Section 4, according to the con- 
ventional theory of paraconductivity, the sheet conductivity in the vicinity 
of the superconducting transition is given by a sum of the electron residual 
conductivity ge 2 (Fermi par t) a nd the conductivity of the Cooper pair fluc- 
tuations (Bose part) (see (101)). This expression is valid in the Ginzburg 
Landau region when the second term is a small correction to the first one. 
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The width of the critical region can be determined from the requirement of 
equality of both contributions in (101)^]: 



1 

165 



l.ZGi 



(2d)- 



(139) 



In accordance with general scaling ideas one can believe that inside the fluc- 
tuation region the conductivity should obey the form: 



<r(T) = ge 2 f 



Gi 



(2d) 



(140) 



18 It is worth mentioning that this definition of the Ginzburg Levanyuk number 
(Gi(2d) — sTjrl) a S rees with that denned from the heat capacity fluctuations 

(Gi(2d) = 
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Concerning the scaling function f(x), we know its asymptotes in the mean 
field region [x ^> 1) and just above the BKT transition (56lM: 



f{x) - { exp [-b(x - x B kt)- 1/2 )] , I - Lkt = -4 ■ (U1) 

The BKT transition temperature T BKT is determined by 
AmT 

n s2 {T c ) = (142) 



and one can find its value by comparing the superfluid density n s from (142) 
with that found in the BCS scheme: 

T BKT = _ AGi y ^ 14 gj 

Here we assumed that the Ginzburg parameter is small, so that the BKT 
transition temperature does not deviate much from the mean-field BCS tran- 
sition temperature T c q. 



Boson mechanism of the T c suppression. The classical and quantum 
fluctuations reduce n s and therefore, suppress T C BKT . At some g — g c ~ 1 , the 
superfluid density n s , and simultaneously T BKT , go to zero. In the vicinity 
of this critical concentration of impurities T BKT -C T c o. Thus a wide new 
window of intermediate temperatures T BKT «r« T c o opens up. In this 
window, according the dynamical quantum scaling conjecture ^4|, one finds 

° = ^(j^kt)- ( 144 ) 



At T - T C BKT < T C BKT the Berezinski-Kosterlitz-Thouless law (|l40[)-([l4l[) 
should hold, so f(x) = f(x) and is exponentially small. In the intermediate 
region T C BKT < T < T BCS the duality hypothesis gives ip(x) = n/2. Let us 
derive this relation. 

We will start from the assumption that in the region T BKT <C T T c q the 
conductivity is a universal function of temperature which does not depend on 
the pair interaction type. Being in the framework of the classical approach, 
let us suppose that in a weak electric field pairs move with the velocity 
v = F/rj, where F = 2eE is the force acting on the pairs. The current 
density j = 2env =crE, (here n is the pair density), so one can relate the 
conductivity with the effective viscosity r\: a=4e 2 n/r]. 

Let us recall that we are dealing with a quantum fluid, so another, su- 
perconducting, view on the problem of its motion near the quantum phase 
transition exists. One can say that with the increase of Gi the role of quan- 
tum fluctuations grows too and fluctuation vortices carrying the magnetic 
flux quantum = f/e are generated. With electric current flow in the sys- 
tem the Lorentz (Magnus) force acts on a vortex: F = j<pQ. The electric field 
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is equal to the rate of magnetic flux transfer, i.e. to the density of the vor- 
tex current: E =<l > on v v v ='Pon v F/r] v , where n v is the density and r\ v is the 
viscosity of the vortex liquid. As a result E =<P^n v j / rj v = j/cr. So one can 
conclude that for vortices the velocity is proportional to the voltage, and the 
force is proportional to the current. For Cooper pairs (bosons) the situation 
is exactly the opposite. 

The duality hypothesis consists in the assumption that at the critical 
point the pair and the vortex liquid density flows are equal: n v v v = nv. Com- 
paring these quantities, expressed in terms of the conductivity from the above 
relations, one can find a universal value for the conductivity at the critical 
point 



One can restrict oneself to a less strong duality hypothesis, supposing the 
product nrj = CT S with a universal S exponent both for the pair and the 
vortex liquids, while the constant C for them is different. In this case, based 
on duality, is possible to demonstrate that S — and the conductivity is 
temperature independent up to T c q but its value is not universal any more 
and can vary from one sample to another. 

To conclude, let us emphasize that in the framework of the boson scenario 
of superconductivity suppression, the BCS critical temperature is changed in- 
significantly, while the "real" superconducting transition temperature T C BKT — 



Fermion mechanism of T c suppression. Apart from the above fluc- 
tuation (boson) mechanism of the suppression of the critical temperature in 
the 2D case, there exists another, fermionic mechanism. The suppressed elec- 
tron diffusion results in a poor dynamical screening of the Coulomb repulsion 
which, in turn, leads to the renormalization of the inter-electron interaction 
in the Cooper channel, and hence to the dependence of the critical temper- 
ature on the value of the high-temperature sheet resistivity of the film. As 
long as the correction to the non-rcnormalized BCS transition temperature 
T c o is still small, one finds f^,pq,B7) : 




At small enough T c q this mechanism of critical temperature suppression turns 
out to be the principal one. The suppression of T c down to zero in this case 
may happen in principle even at </ > 1. A renormalization group analysis 
gives |Q the corresponding critical value of conductance 



0. 





(146) 




(147) 
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Here we should recall that the typical experimental ]8!|] values of g c are in 
the region g c ~ 1—2, and do not differ dramatically from the predictions of the 
boson duality assumption g c = ~ . If one attempts to explain the suppression 
of T r within the fermion mechanism, one should assume that In tJ— > 5. 
Then, according to Eq. (147), g c > 2/n and the boson mechanism is not 
important. On the contrary, if In -^-^ < 4, then Eq. ( |146| ) gives a small 
correction for T c even for g c — 2/tt and the fermion mechanism becomes 
unimportant. The smallness of the critical temperature T c compared to the 
Fermi energy is the cornerstone of the BCS theory of superconductivity and it 
is apparently satisfied even in high-T c materials. Nevertheless it is necessary 
to use the theoretically large logarithmic parameter with care, if one needs 
In Tj^p to be as large as 4. 



5 Microscopic derivation of the Time-Dependent 
Ginzburg-Landau equation 



5.1 Preliminaries 



We have seen above how the phenomenological approach based on the GL 
functional allows one to describe fluctuation Cooper pairs (Bose particles) 
near the superconducting transition and to account for their contribution to 
different thermodynamical and transport characteristics of the system. Now 
we pass to the discussion of the microscopic description of fluctuation phe- 
nomena in superconductors. The development of the microscopic approach is 
necessary for the following reasons: 

1. This description permits microscopic determination of the values of the 
phenomenological parameters of the GL theory. 

2. This method is more powerful than the phenomenological GL approach 
and permits treatment of fluctuation effects quantitatively even far from the 
transition point and for magnetic fields strong as H C 2, taking into account 
the contributions of dynamical and short wavelength fluctuations. 

3. The electron energy relaxation times in metals are relatively large [r e 3> 
h/T) which causes the electron low frequency dynamics to be sensitive to 
the nearness to the superconducting transition. This is why the temperature 
dependence of fluctuation corrections can be determined generally speaking 
not only by the Cooper pair motion but also by changes in the single-electron 
properties. 

4. There are some fluctuation phenomena in which the direct Cooper pair 
contribution is considerably suppressed or even absent altogether. Among 
them we can mention the nuclear magnetic relaxation rate, tunnel conduc- 
tivity, c-axis transport in strongly anisotropic layered metals, thermoelectric 
power and heat conductivity where the fluctuation pairing manifests itself by 
means of the indirect influence on the properties of the single-particle states 
of electron system. 
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Formally in the above consideration averaging over the superconducting 
order parameter has been accomplished by means of a functional integration 
over all its possible bosonic field configurations. In this description we have 
dealt with the fluctuation Cooper pair related effects only and the method of 
the functional integration turned out to be simple and effective for their de- 
scription. In the following Sections we will develop the diagrammatic method 
of Matsubara temperature Green functions which is more adequate for the 
description of the properties of a Fermi system of interacting electrons. 



5.2 The Cooper channel of electron-electron interaction 

Let us start the microscopic description of fluctuation phenomena in a su- 
perconductor from the electron Hamiltonian. We will choose it in the simple 
BCS form[3 

H=J2 E (p)$Up,* + 9 V'p+q,a^p,-aV'-P',- CT 'V'p'+q, CT '-(148) 

The momentum conservation law side by side with singlet pairing are already 
taken into account in the interaction term. Here E(jp) is the quasiparticle 
spectrum of the normal metal; g is the negative constant of electron-electron 
attraction which is supposed to be momentum independent and different 
from zero in a narrow domain of momentum space in the vicinity of the 
Fermi surface where 

Pf < P , P < Pf H • 

vf vp 

4>p a and V'p.o- are the creation and annihilation field operators in the Heisen- 
berg representation, so the first term is just the kinetic energy of the non- 
interacting Fermi gas. The interaction term is chosen in the traditional form 
characteristic for the electron-phonon mechanism of superconductivity^]. 
For the description of the properties of an interacting electron system with 



the Hamiltonian (148) we will use the formalism of the Matsubara tempera- 



We suppose that reader is familiar with the BCS formulation of the theory of 
superconductivity (see for example, @). 

Fluctuations in the framework of more realistic Eliashberg [Q model of super- 
conductivity were studied by B.Narozhny He demonstrated that the strong 
coupling does not change drastically the results of the weak coupling approxima- 
tion. The critical exponents turn out to be exactly the same as in the framework 
of the GL theory, which provides an adequate description of paraconductivity in 
strong coupling superconductors. The robustness of the critical exponents and 
their dependence in GL region on the space dimensionallity only was stressed 
in J&l in relation to the discussion of the paraconductivity at the edge of the 
superconductor-insulatortransition. 
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ture diagrammatic technique. The state of a non-interacting quasiparticle is 
described by its Green function 

g(p ' £ " ) = ^?(p)' (149) 

where e n = (2n + l)nT is a fermion Matsubara frequency and £(p) = 
E(p)—Ep is the quasiparticle energy measured from the Fermi level. 

As it is well known the effective electron-electron attraction leads to a 
reconstruction of the ground state of the electron system which formally 
manifests itself by the appearance at the critical temperature of a pole in the 
two particle Green function 

where T T is the time ordering operator and 4D vector notations are used |)2| . 
As it well known the two particle Green function can be expressed in terms 
of the vertex part |9^] . In the case under consideration it is the vertex part of 
the electron-electron interaction in the Cooper channel L(q, f?k), which will 
be called below the fluctuation propagator. The Dyson equation for L(q, 14), 
accounting for the e-e attraction in the ladder approximation, is represented 
graphically in Fig.^J . It can be written down analytically as 

L- 1 (q,n k ) = g- 1 -n(q,[2 k ), (150) 

where the polarization operator i7(q, 14) is defined as a loop of two single- 
particle Green functions: 

77(q,/2 fe )=T^ J ^ G (p + q, e „ +fc )G(-p,e-„). (151) 

Let us emphasize, that the two quantities introduced above, C(p,p',q) 
and L(q), are closely connected with each other. The former being integrated 
over momenta p and p' becomes an average of the product of two order 
parameters: 

dpdp'£(p,p',q) = ^ (A q A* q ), (152) 

where A q is the superconducting gap proportional to the condensate wave 
function ty. Thus, this quantity represents the coefficient in the linear term 
in the GL equation. In terms of the polarization operator introduced above 
it can be written as 



/ 



dpdp' C (p,p' , q) oc 



i- g n 
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Fig. 6. The Dyson equation for the fluctuation propagator (wavy line) in the lad- 
der approximation. Solid lines represent one-electron Green functions, bold points 
correspond to the model electron-electron interaction. 



Comparing this equation with Eq. ( |150| ) for the fluctuation propagator we see 
that the corresponding expressions are very similar. After analytical contin- 
uation to the real frequencies the fluctuation propagator L(q,ifi) coincides 
with the quantity defined by Eq.(152) (up to a constant). 

One can calculate the propagator (15C) using the one-electron Green func- 
tions of the normal metal ( |149| ) . For sake of convenience of future calculations 
let us define the correlator of two onc-clcctron Green functions 



V(q,ei,s 2 )= / — —3-G(p + q,£i)G(-p,£ 2 ) 
J (27r) 

/ 1 

= 2ltV<9{ — £\£~i) 



(153) 



ex -e 2 | +iZ\£(q,p)| £(p)=i 



f / F.S. 



where 0(— £\£2) is Heavyside step function, v is the one-electron density of 
states, <>f.s.= J -ajr means the averaging over the Fermi surface, 



^(q,p) 



= K(q + P) - £{-p)]\e(p)=E F ~(v p q) ( 



(p)=o- 



(154) 



The last approximation is valid not too far from the Fermi surface, i.e. when 



It is impossible to carry out the angular averaging in (153) for a general 
anisotropic spectrum. Nevertheless in the following calculations of fluctuation 
effects in the vicinity of critical temperature only small momenta v p q <C 
T will be involved in the integrations, so we can restrict our consideration 
here to this region, where one can expand the integrand in powers of v p q. 
Indeed, the presence of 0(— £\£2) leaves the difference of the two fermionic 
frequencies in (|153|) to be of the order of the temperature which permits 
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this expansion. The first term in v p q will evidently be averaged out, so with 
quadratic accuracy one can find: 

7>(q, s u e 2 ) = ^v^—^ (I - 2 ^ j . (155) 

Now one can calculate the polarization operator 



J7( q , n k ) =Tj2n^ 



(156) 



E- 

n>0 



i /9 4- HM 

1 / Z + 4ttT 



- 2 



((y P q) 2 )F.s. 

(4ttT) 2 



X 

E 



'o(„+l/2 + J^ 



The calculation of the sums in (156) can be carried out in terms of the 
logarithmic derivatives of the P-function ^ n \x). It worth mentioning that 
the first sum is well known in the BCS theory, one can recognize in it the 
so-called "Cooper logarithm"; its logarithmic divergence at the upper limit 
(ip(x ^> 1) w In a;) is cut off by the Debye energy (N max = f^) and one gets: 



^77(q, O k ) = i> Q 



\oa 



4ttT 
((v P q) 2 



up 
2^f 



4> 



F.S. 



(4ttT) 2 



1 

2 



K4j 

4ttT 



(157) 



4vrT/ 



The critical temperature in the BCS theory is determined as the temperature 
T c at which the pole of L(0,0,T c ) occurs 



L _1 (q =0, O k = 0,T C 



77(0, 0, T c ) =0, 



T 2 7B / 1 
i c = cjd exp 

7T \ 1/0 



(158) 



where 75 = 1.78 is the Euler constant. Introducing the reduced temperature 
e = In (7^) one can write the propagator as 



L~\q, O k ) 



1 



Ok 



e + tb(- + 

^ v 2 4ttT' 



<(v P q) 2 



"F.5. 



(4ttT) : 



1> 







a- 


471-17 



We found (15E) for bosonic imaginary Matsubara frequencies iO k = 2niTk. 
These frequencies are necessary for the calculation of fluctuation contribu- 
tions to any thermodynamical characteristics of the system. 

In the vicinity of the transition point one can restrict oneself in sum- 
mations of the expressions with £(q, Ok) ovcr Matsubara frequencies to the 
so-called static approximation, taking into account the term with O k = 



.(159) 
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only, which turns out to be the most singular term in t < 1 . This ap- 
proximation physically means that the product of Heisenberg field operators 
ip PtCr il)-p^- a appears here like a classical field which in the phcnomcnolog- 
ical approach describes the Cooper pair wave function and in the vicinity of 
critical temperature is proportional to the fluctuation order parameter. Hav- 
ing in mind namely this GL region of temperatures we restricted ourselves 
above by the assumption of small momenta v p q <C T. In these conditions 
the static propagator reduces to 

L(q ' 0) = ^7rW- (160) 

With an accuracy of a numerical factor and the total sign this correlator 
coincides with the expression for (l^ql 2 )- By this expression we also have 
finally obtained the microscopic value of the coherence length £ for a clean 
superconductor with an isotropic D-dimensional Fermi surface which was 
often mentioned previously (compare with (|^)) 

2 7C(3)vg. 

In order to describe the fluctuation contributions to transport phenomena 



one has to start from the analytical continuation of the propagator (159) 
from the discrete set of Qk > to the whole upper half-plane of imaginary 
frequencies. The analytical properties of (x) — functions (which have poles 
at x = 0, — 1, —2...) permit one to obtain the retarded propagator L R (q, —if}) 
by simple substitution ifl^ — > fl . For small Q <C T the ^—functions can be 
expanded in —ifljAnT and the propagator acquires the simple pole form : 

1 1 1 

L*(q, Q) = , 2 = — — 8Tc2 2V (162) 

^-Jf^ + e + £ 2 q 2 TXViQ- (r G ^ + ^-£ 2 q 2 ) 

This expression provides us with the microscopic value of the GL relaxation 
time tgl — 8(t-t ) ' w idely used above in the phenomenological theory. 



Moreover, comparison of the microscopical ly d erived (162) with the phe 



nomenological expressions (|36|), ( |90| ) and (109) shows that aT c — v and 
Igl = ttv/8T c . 

In evaluating L(q, Qk) we neglected the effect of fluctuations on the one- 
electron Green functions. This is correct when fluctuations are small, i.e. not 
too near to the transition temperature. The exact criterion of this approxi- 
mation will be discussed in the following. 



5.3 Superconductor with impurities 



Account for impurities. In order to study fluctuations in real systems 
like superconducting alloys or high temperature superconductors one has to 
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Fig. 7. The equation for the vertex part X(q, uii, W2) in the ladder approximation. 
Solid lines correspond to bare one-electron Green functions and dashed lines to the 
impurity random potential correlators. 



perform an impurity average in the graphical equation for the fluctuation 
propagator (see Fig. [|). This procedure can be done in the framework of the 
Abrikosov-Gorkov approach M] , which we shortly recall below. 

Let us start from the equation for the electron Green function in the 
potential of impurities U(r): 



(e-U{v)-H) G e {t,t')=5{v-v') 



(163) 



If we solve this equation using the perturbation theory for the impurity poten- 
tial and average the solution, then the average product of two Green func- 
tions, can be presented as series, each term of which is associated with a 
graph drawn according the rules of diagrammatic technique (see Fig. (?]). In 
this technique solid lines correspond to bare Green functions and dashed lines 
to random potential correlators. We assume that the impurity system random 
potential U(r) is distributed according to the Gauss S— correlated law. Then 
all the correlators can be represented as the products of pair correlators 



(U(r))=0, {U{r)U(r'))={U 2 )8{r~r') 1 



(164) 



where the angle brackets denote averaging over the impurity configuration. 
Equation 164 corresponds to the Born approximation for the electron interac- 



tion with short range impurities, and (U 2 ") — Ci rnp (J V(r)dr) where Ci mp 
is the impurity concentration and V(r) is the potential of the single impurity. 
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In conductors (far enough from the metal-insulator transition) the mean 
free path is much greater than the electron wavelength / 3> A = 2-k /pp (which 
in practice means the mean free path up to tens of interatomic distances). 
As is well known j9^] for the electron spectra with dimensionality D > 1 
the angular integration in momentum space reduces considerably the con- 
tribution of the diagrams with intersecting impurity lines what permits to 
omit them to the leading approximation in (ppl) ■ In this approxim ation 



the one-electron Green function keeps the same form as the bare one (149) 
with the only substitution 

£ n =>• £ n = £« + ^-sign(s n ), (165) 
It 

where 1/r = 2-kv (U 2 ) is the frequency of elastic collisions. 

Another effect of the coherent scattering on the same impurity by both 
electrons forming a Cooper pair is the renormalization of the vertex part 
A(q, £i, £2) in the particle-particle channel. Let us demonstrate the details 
of its calculation. The renormalized vertex A(q, £i,£2) is determined by a 
graphical equation of the ladder type (see Fig. j?j ). Here after the averaging 
over the impurity configurations the value (U/ = -^hr * s associated with 
the dashed line. In the momentum representation this, generally speaking, 
integral equation is reduced to the algebraic one 

A-^cbe^ea) = 1 - - Pfaei.eb), (166) 

l-KVT 



where P(q, £1,62) was defined above by ( |153| ) 

Now one has to perform a formal averaging of the general expression (153) 
over the Fermi surface (( ...) f.s.)- Restricting consideration to small momenta 

^(q,p)|| P |=p, « |ei-e 2 |. (167) 
the calculation of A(q, loi,u>2) for the practica lly i mportant case of an ar- 



bitrary spectru m can be done analogously to ( 155 ). Indeed, expanding the 



denominator of (153) one can find 



A(q,wi,w 2 ) = 77-7-7 — r, 757 . (168) 

l £ l - £ Z\ + rl^-cDaP <y(-£l£ 2 J 

It is easy to see that assumed restriction on momenta is not too severe and is 
almost always satisfied in calculations of fluctuation effects at temperatures 
near T c . In this region of temperatures the effective propagator momenta 
are determined by |q| e // ~ [£, GL {T)] _1 = ^ while the Green 

function q-dependence becomes important for much larger momenta q ~ 



min{£ ,/ }, which is equivalent to the limit of the condition (167) 
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The average in ( |168| ) can be calculated for some particular types of spec- 
tra. For example in the cases of 2D and 31? isotropic spectra it is expressed 
in terms of the diffusion coefficient £>(_d) : 

((^(q,p)l|p|=p F ) 2 >i,s.(tf) =r- 1 V (D) q 2 = ^f. (169) 

Another important example is already familiar case of quasi-two-dimensional 
electron motion in a layered metal: 

Z{p)=E( Pll )+Jcos(p z s)-E F , (170) 

where -E(p||) = p|/(2m), p = (p\\,p z ), P\\ = (px,P y ), J is the effective 
nearest-neighbor interlayer hopping energy for quasiparticles. We note that 
J characterizes the width of the band in the c-axis direction taken in the 
strong-coupling approximation and can be identified with the effective energy 
of electron tunneling between planes (see (|6^) and footnote 14) . The Fermi 
surface, defined by the condition £(p) = 0, is a corrugated cylinder (see Fig. 



In this case the average (169) is written in a more sophisticated form: 



((Z}£(q,p)|| p | =PF ) 2 >F.s. = ^(v 2 F q 2 + 4J 2 sin 2 (q^^^r^Vq 2 , (171) 

where we have introduced the definition of the generalized diffusion operator 
T> in order to deal with an arbitrary anisotropic spectrum. 



Propagator. In Section 4, in the process of the microscopic derivation of 
the TDGL equation, the fluctuation propagator was introduced. This object 
is of first importance for the microscopic fluctuation theory and it has to 
be generalized for the case of an impure metal with an anisotropic electron 
spectrum. This is easy to do using the averaging procedure presented in 



the previous Section. Formally it is enough to use in equation (150) the 
polarization operator 77(q, f2k) averaged over impurity positions, which can 
be expressed in terms of P(q, £n+fc, £-n) introduced above: 

77(q, Qk) — M 1 ?' £ n+k, £-„)"P(q, uJ n+k , w_ n ) = (172) 

UJ n 



[P(q,e„ +fc ,e_ 



-l l 



For relatively small q (A^(q,p)\\ E ( p )\=E F < \ £ n+k - £-n| ~ max{T, r 1 }) 
and one can find an expression for the fluctuation propagator, which 

can be useful in studies of fluctuation effects near T c (e <C 1) for the dirty 



and intermediate but not very clean case (TV 1/^/e). Expanding (168) in 
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Fig. 8. The Fermi surface in the form of a corrugated cylinder 



powers of (/\£(q, p)||e( p )|=e f /|2£„ + &k\) 2 it is p ossible write L R (q, Q) in a 
form almost completely coinciding with Exp. (162): 



1 



L*(q,/2)=-~ 



ve-m+e{T T )c? 



(173) 



Let us stress that the phenomenological coefficient 70 l turns out to be equal 
to the same value |^ as in clean case, and hence does not depend on the 
impurity concentration. The only difference in comparison with the clean 
case is in appearance of a dependence of the natural effective coherence length 
on the elastic relaxation time. In the isotropic D-dimensional case it can be 
written as 



%d)(Tt) 



(4maT) 1 = i] (D) 



(174) 



D 



Hi 



1 



4itTt 



4ttTt y v 2 ' 



(we introduced here the parameter 77(D) frequently used in the microscopic 
theory)^. 

21 Let us recall that its square determines the product of the GL parameter a and 
the Cooper pair mass entering in the GL functional. In clean case we supposed 
the letter equal to two free electron masses and defined a in accordance with 
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The generalization of (fl73j) for the case of a layered electronic spectrum 
is evident: 

L R (q,f2) = = 1 T/ — . (175) 

"e-t^+TJpoqfi +r S m 2 (q z s/2) 



One has to remember that the Exp. (173) was derived in the assumption 
of small momenta A£(q,-p)\\ E ( p )\=E F < |?n+fc - £-n| ~ max{T, r" 1 }, so 
the range of its applicability is restricted to the GL region of temperatures 
e = ln(5-) <C 1, where the integrands of diagrammatic expressions have 
singularities at small momenta of the Cooper pair center of mass. 

Finally let us express the Ginzburg-Levanyuk parameter for the important 
2D case in terms of the microscopic parameter 77(2) • In accordance with ( p7|) 
and the definition ( |174| ): 

G H2) (Tt) = ^ = 1 p= r - (176) 
I67H mT c r](2) (Tt) 

One can see that this general definition in the limiting cases of a clean and 
dirty metal results in the same values Gi^c) and Gir^d) &s was reported in 
Table 1. 



5.4 Microscopic theory of fluctuation conductivity of layered 
superconductor 

Qualitative discussion of different fluctuation contributions In Sec- 
tion 4 the direct fluctuation effect on conductivity, related with the charge 
transfer by means of fluctuation Cooper pairs, was discussed in details. Never- 
theless, below in this Section we return to its discussion and will demonstrate 
its calculation by means of the microscopic theory. This will be done in pur- 
pose to prepare the basis for studies of the Aslamazov-Larkin contribution in 
the variety of physical values like magnetoconductivity near the upper crit- 
ical field, conductivity far from transition point, fluctuation conductivity in 
ultra-clean limit, Hall conductivity etc. 

Microscopic approach permits also to calculate the cited above indirect 
fluctuation effects like so called DOS and MT contributions. We will start 
now from their qualitative discussion. 

The important consequence of the presence of fluctuating Cooper pairs 
above T c is the decrease of the one-electron density of states at the Fermi level. 
Indeed, if some electrons are involved in pairing they can not simultaneously 

([L8|). As we just have seen in the case of the impure superconductor £ depends 
on impurity concentration and this dependence, in principle, can be attributed 
both to a or m. For our further purposes it is convenient to leave a in the same 
form ( |l8| ) as in the case of a clean superconductor. The Cooper pair mass in this 
case becomes dependent on the electron mean free pass what physically can be 
attributed to the diffusion motion of the electrons forming the pair. 
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participate in charge transfer and heat capacity as single-particle excitations. 
Nevertheless, the total number of the electronic states can not be changed 
by the Cooper interaction, and only a redistribution of the levels along the 
energy axis is possible |93|,[)4j . In this sense one can speak about the opening 
of a fluctuation pseudo-gap at the Fermi level. The decrease of the one- 
electron density of states at the Fermi level leads to a reduction of the normal 
state conductivity. This, indirect, fluctuation correction to the conductivity 
is called the density of states (DOS) contribution and it appears side by 
side with the paraconductivity (or Aslamazov-Larkin contribution). It has the 
opposite (negative) sign and turns out to be much less singular in (T— T c ) _1 
in comparison with the AL contribution, so that in the vicinity of T c it 
was usually omitted. However, in many cases p9| , p5 96,9^,^,^9), when for 



some special reasons the main, most singular, corrections are suppressed, the 
DOS correction becomes of major importance. Such a situation takes place 
in many cases of actual interest (quasiparticle current in tunnel structures, 
c-axis transport in strongly anisotropic high temperature superconductors, 
NMR relaxation rate, thermoelectric power). 

The correction to the normal state conductivity above the transition tem- 
perature related with the fluctuation DOS renormalization for the dirty su- 
perconductor can be evaluated qualitatively. Indeed, the fact that some elec- 
trons (AA/" e per unit volume) participate in fluctuation Cooper pairing means 
that the effective number of carriers taking part in one-electron charge trans- 
fer diminishes leading to a decrease of conductivity (we deal here with the 
longitudinal component): 

m to 

where n s is the superfluid density coinciding with the Cooper pairs concen- 
tration. The latter can be identified with the average value of the square of 
the order parameter modulus already calculated as the correlator 

(<T (0)^(r)) (T > T c0 ) = £ (|f fe | 2 ) exp (ikr) = Q^) ,(178) 

with r ~ £. For the 2D case, which is of the most interest to us, one finds: 
1 , ^_7C(3)^,_ 1 



Us 



Aira^ 2 s 



:Ko(^) = 4V— In - (179) 



where we have used the explicit expression (|18|) for a and £. As we will see the 
corresponding expression for the fluctuation DOS correction to conductivity 



(177) coincides with the accuracy of 2 with the microscopic expression (191) 
which will be carried out below. 

The third, purely quantum, fluctuation contribution is generated by the 
coherent scattering of the electrons forming a Cooper pair on the same elastic 
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impurities. This is the so called anomalous Maki- Thompson (MT) contribu- 
tion |(|[7j which can be treated as the result of Andreev scattering of the 
electron by fluctuation Cooper pairs. This contribution often turns out to be 
important in conductivity and other transport phenomena. Its temperature 
singularity near T c is similar to that of the paraconductivity, although being 
extremely sensitive to electron phase-breaking processes and to the type of 
orbital symmetry of pairing it can be suppressed. Let us evaluate it. 

The physical origin of the Maki-Thompson correction consists in the fact 
that the Cooper interaction of electrons with the almost opposite momenta 
changes the mean free path (diffusion coefficient) of electrons. As we have 
already seen in the previous Section the amplitude of this interaction increases 
drastically when T — ► T c : 

g 1 T _ 1 

9eff ~ l-gln^p - hT£ ™ T~^T C ~ 7 

What is the reason of this growth? One can say that the electrons scatter 
one at another in resonant way with the virtual Cooper pairs formation. Or 
it is possible to imagine that the electrons undergo the Andreev scattering at 
fluctuation Cooper pairs binding in the Cooper pair themselves. The proba- 
bility of such induced pair irradiation (let us remind that Cooper pairs are 
Bose particles) is proportional to their number in the final state, i.e. n(p) ((?]). 
For small momenta n(p) ~ 1/e. 

One can ask why such interaction does not manifest itself considerably 
far from the transition point? The matter of fact that so intensively interacts 
just small number of electrons with the total momentum q < £ _1 (T). In 
accordance with the Heisenberg principle the minimal distance between such 
electrons is of the order of ~ £,(T). From the other hand such electrons in 
purpose to interact have to approximate one another up to the distance of 
the Fermi length Xp ~ 1/pf- The probability of such event may be estimated 
in the spirit of the self-intersection trajectories contribution evaluation in the 



weak localization theory [100 102] ■ 

In the process of diffusion motion the distance between two electrons 
increases with the time growth in accordance with the Einstein law: R(t) ~ 

1/2 

(Dt) ' " . Hence the scattering probability 

W ~L TkT) VFdt - 

The lower limit of the integral can be estimated from the condition i?(£ m i n ) 
~ £(T) (only such electrons interact in the resonant way). The upper limit 
is determined by the phase breaking time since for larger time intervals 
the phase coherence, necessary for the pair formation, is broken. In result 
the relative correction to conductivity due to such processes is equal to the 
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product of the scattering probability on the effective interaction constant: 
5a MT /a = W geff-ln the 2D case 

6a MT --In ^ 



em 



This result will be confirmed below in the frameworks of the microscopic 
consideration. 



Generalities Let us pass to the microscopic calculation of the fluctuation 
conductivity of the layered superconductor. We begin by discussing the quasi- 
particle normal state energy spectrum. While models with several conducting 
layers per unit cell and with either intralayer or intcrlayer pairing have been 



considered |l03|| , it has been shown [104] that all of these models give rise to 



a Josephson pair potential that is periodic in k z , the wave- vector component 
parallel to the c-axis, with period s, the c-axis repeat distance. While such 
models differ in their superconducting densities of states, they all give rise 
to qualitatively similar fluctuation propagators, which differ only in the pre- 
cise definitions of the parameters and in the precise form of the Josephson 
coupling potential. Ignoring the rather unimportant differences between such 
models in the Gaussian fluctuation regime above T C {H), we therefore con- 
sider the simplest model of a layered superconductor, in which there is one 
layer per unit cell, with intralayer singlet s-wave pairing. These assumptions 



lead to the simple spectrum (170) and hence to a Fermi surface having the 
form of a corrugated cylinder (see Fig.|J). 

Some remarks regarding the normal-state quasiparticle momentum relax- 
ation time are necessary. In the "old" layered superconductors the materials 
were generally assumed to be in the dirty limit (like TaS^ (pyridine) 1/2)- I n 
the high-T c cuprates, however, both single crystals and epitaxial thin films 
are nominally in the "intermediate" regime, with l/£ X y ~ 2 — 5. In addition, 
the situation in the cuprates is complicated by the presence of phonons for 
T ~ T c ~ 100-ftT, the nearly localized magnetic moments on the Cu 2+ sites, 
and by other unspecified inelastic processes. In this Section we assume simple 
elastic intralayer scattering and restrict our consideration to the local limit 
in the fluctuation Cooper pair motion. This means that we consider the case 
of not too clean superconductors, keeping the impurity concentration m and 
reduced temperature such that the resulting mean-free path satisfies the re- 
quirement I < £ X y{T) = %| and the impurity vertex can be taken in the local 



form (168) with ((A£ (q, p)|) 2 )f.s. determined by (171). The phase-breaking 
time Ttp is supposed to be much larger than r. 

The most general relation between the current density j(r,t) and vector- 
potential A(r ,t') is given through the so-called electromagnetic response 
operator Q a p(r, r', t, t') §1: 



j(r,t) = J Q Q(3 (r,r',M')A(r',0 dr' dt, 



70 Anatoly Larkin and Andrei Varlamov 

Assuming space and time homogeneity, one can take the Fourier transform 
of this relation and compare it with the definition of the conductivity tensor 
ja = CapEp. This permits us to express the conductivity tensor in terms of 
the retarded electromagnetic response operator 

<? af 3(u) = - — [Q aP ] R (uj) . (180) 
iuj 

The electromagnetic response operator Q a f){wv), defined on Matsubara 
frequencies lo v = (2v + l)itT, can be presented as the correlator of two ex- 
act one-electron Green functions ]92f averaged over impurities and accounting 
for interactions, in our case the particle-particle interactions in the Cooper 
channel. The appropriate diagrams corresponding to the first order of per- 
turbation theory in the fluctuation amplitude are shown in Fig. ([)]). 

With each electromagnetic field component A a we associate the external 
vertex ev a (p) = e ■ For the longitudinal conductivity tensor elements 
(parallel to the layers, for which a = x,y), the resulting vertex is simply 
ep a /m. For the c-axis conductivity, the vertex is given by 

ev z (p)=e^- = -eJssin(p z s). (181) 

Each solid line in the di agram s represents a one-electron Green function 
averaged over impurities (|149|), a wavy line represents a fluctuation propa- 



gator i(q, f2 k ) (J1T3J) , three-leg vertices were defined by the Exp.( |l68| ). The 
four-leg impurity vertex, appearing in diagrams 3-4, 9-10 of the Fig. p| i s 
called the Cooperon in the weak localization theory (see, for example, Jl05| ) 
It is easy to see that it differs from the above three-leg vertex only by the 
additional factor (2-kvt)~ 1 . We do not renormalize the current vertices: it is 
well known (sec ]9^|) that this renormalization only leads to the substitution 
of the scattering time r by the transport one Tt r . We integrate over the inter- 
nal Cooper pair momentum q and electron momentum p and sum over the 
internal fermionic and bosonic Matsubara frequencies, with momentum and 
energy conservation at each internal vertex (fluctuation propagator endpoint) 
in the analytical expressions for the diagrams presented in Fig. ([)]). 

After these necessary introductory remarks and definitions we pass to the 
microscopic calculation of the different fluctuation contributions. 



Aslamazov-Larkin contribution We first examine the AL paraconduc- 
tivity (diagram 1 of Fig.(^|)). Actually this contribution was already studied 
in the Section 4 in the framework of the TDGL equation but, in order to 
demonstrate how the method works, we will carry out here the appropriate 
calculations in the microscopic approach, as was originally done by Aslama- 
zov and Larkin ||. 
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Fig. 9. Feynman diagrams for the leading-order contributions to the fluctuation 
conductivity. Wavy lines are fluctuation propagators, thin solid lines with arrows 
are impurity-averaged normal-state Green's functions, shaded semicircles are vertex 
corrections arising from impurities, dashed lines with central crosses are additional 
impurity renormalisations and shaded rectangles are impurity ladders. Diagram 1 
represents the Aslamazov-Larkin term; diagrams 2-4 represent the Maki-Thompson 
type contributions; diagrams 5-10 arise from corrections to the normal state density 
of states. 
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The AL contribution to the electromagnetic response operator tensor has 
the form: 

Q£p(u») = 2e 2 T^ J ^!S_i? Q (q,a,^)i(q,a) x (182) 

xB)3(c[,n k ,uj u )L(q,n k +lu v ), 
where the three Green function block is given by 

B a (c[, Q k , u v ) = T^2\(q,e n+U ,n k - e„)A(q, e n , fl k - e„) x (183) 

f d 3 p 

x J -7^T3Ua(p)G(p,en +I /)G(p,e„)G(q-p, fl k - e n ). 

Expanding G(q — p, Q k — e n ) over q one find that the angular integration 
over the Fermi surface kills the first term and leaves nonzero the second 
term of the expansion only. Then the £— integration is performed by means 
of the Cauchy theorem. The further summation over the fermionic frequency 
is cumbersome, so we will show it for the example of the simplest case of 
a dirty superconductor with Tt <C 1. In this case the main sources of the 



e n — dependence in (183) are the A-vertices and that originating from the 
Green functions can be neglected by the parameter Tt -C 1 (indeed, one can 
see that e n ~ T are important in vertices, while in Green functions e n > t _1 
only). The remaining summation in (|183|) is performed in the same way as 



was done in (173) and results in 



B a (q,f2k,to v ) = v-^r (vaqpvp) 



FS 



(184) 



8T 



\-cj u + Dq 2 _ I \n k \+Dq 2 
AttT ' V[ 2 4ttT ' 



Dq 2 



,,i , \n k+v \ +uj u + D q 2 i \\n k+u \ + 

4^ } -^ ( 2 + 4^T 



Now let us return to the general expression for Qaai^i/) an< i transform 
the f2 k — summation into a contour integral, using the identity [106] 



Tj2fm = ^ jf dzcoth^/(-iar), 



where z = iQ k is a variable in the plane of complex frequency and the contour 
C encloses all bosonic Matsubara frequencies over which the summation is 
carried out. In our case the contour C can be chosen as a circle with radius 
tending to infinity (see Fig. fio| ): 
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Fig. 10. The contour of integration in the plane of complex frequencies. 



e f da f z 
Q^W = -j- J j> dz coth — B a (q, -iz + u v , -iz) x (185) 

xi(q, — iz)Bp(q L , —iz + ui v , —iz)L(q, —iz + u) v ). 



One can see that the integrand function in ( 18E ) has breaks of analyticity 
at the lines Imz = and Imz = —%u) v . Indeed, the fluctuation propaga- 
tor L(q, Qk) and Green function blocks i? Q (q, fik, w^) were defined on the 
bosonic Matsubara frequencies only, while now we have to use them as func- 
tions of the continuous variable z. As it is well known from the properties 
of Green functions in the complex plane z , two analytical functions, related 
with L(q, J2fc), can be introduced. The first one, L R (q, — iz) (retarded), is 
analytic in the upper half- plane (Imz > 0), while the second one, L (q, —iz) 
(advanced), has no singularities in the lower half-plane (Imz < 0). As we 
have seen above the same lines separate the domains of the analyticity of the 
Green function blocks, so the functions B RR , B RA , B AA analytic in each do- 
main can be introduced (with the appropriate choices of the | fik+v \ and | fik \ 
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signs in the arguments of the ip— functions, see (184)). This means that by 
cutting the z-plane along the lines Imz = and Imz = — iu> u we can reduce 
the calculation of the contour integral in (185) to the sum of three integrals 
along the contours C\ , C2 , C3 which enclose domains of well defined analyt- 
icity of the integrand function. The integral along the large circle evidently 
vanishes and the contour integral is reduced to four integrals along the cuts 
of the plane in Fig.ffl 



I(q,u u )= f dzcoth — B a (q, -iz)L(q, -iz)Bp(q, -iz)L(q, -iz + lo v ) 

JC1+C2+C3 11 

= / dz coth — L*(q,-i* + w„) [B RR Bf R L R (q,-iz)- 

J — OO ^ 

/OO — ZUJ U 
dz coth ±=L A (q, -iz)x 
-OQ—iu) u ^ 

[B RA B RA L R (d, -iz + u„) - B AA B AA L A (q, -iz + u v )] .. 

Now one can shift the variable in the last integral to z = z' — iuj u , take into 
account that iu> v is the period of coth ^ and get an expression analytic in 

iuj u — ► LU. 

In the v icini ty of T c , due to the pole structure of the fluctuation prop- 
agators in ( 182 ), the leading contribution to the electromagnetic response 
operator Q^a^ arises from them rather than from the frequency depen- 
dence of the vertices B a , so we can neglect the fl^- and w^-dependencies of 
the Green functions blocks and use the expression for B a (q, 0, 0) valid for 
small q ab only: 

fl a (q) = -2i/5W/ vh a , a = X ,y 
Vp [sj smg z s, a — z 

Detailed calculations demonstrate that this result can be generalized to an 
arbitrary impurity concentration just by using the expression (174) for T]( 2 )- 
Finally: 



+L A (q, -iz + iu>)] ImL^q, -iz). 

Being interested here in the d.c. conductivity one can expand the inte- 
grand function in to. It is possible to show that the zeroth order term is 
cancelled by the same type contributions from all other diagrams (this can- 
cellation confirms the absence of anomalous diamagnetism above the critical 
temperature). The remaining integral can be integrated by parts and then 
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carried out taking into account that the contribution most singular in e comes 
from the region z ~ e <C T : 



•s 



(27T) 2 [(r ?(2) q2 + e )( ?7(2) q2 + e + r)]3/2 
e 2 1 _eMl/>/iF, e«r 

16s [efe + rJjVa 16s \ 1/e, e>r ' 1 ; 

where the Lawrence-Doniach anisotropy parameter r p7| was already defined 

by®. 

In the same way one can evaluate the AL contribution to the transverse 
fluctuation conductivity J^Mpi : 



al = 7re 2 sr 2 f d 2 q 1 = 

32 J (2tt) 2 [( r/(2)q 2 +e ) (r?(2)q 2 + e + r) ]3/ 2 ^ > 



e 2 s ( e + r/2 \ e 2 s ( y/r/e, for e <C r 

V fefe + r)l 1 /2 ~ J ~* "64^ 1 fr/2e) 2 , for e > r 



32ry (2) V [e(e + r)]V2 y 64t7 (2) 1 (r/2e)' , for e » 

Note, that contrary to the case of in-plane conductivity, the critical exponent 
for a zz above the Lawrence-Doniach crossover temperature Tld (for which 
c(Tld) = r) is 2 instead of 1, so the crossover occurs from the OD to 3D 
regimes. This is related with the tunneling (so from the band structure point 
of view effectively zero dimensional) character of electron motion along the 
c-axis. 



Contributions from fluctuations of the density of states In original 
paper of Aslamazov and Larkin || the most singular AL contribution to 
conductivity, heat capacity and other properties of a superconductor above 
the critical temperature was considered. The diagrams of the type 5-6 were 
pictured and correctly evaluated as less singular in e. Nevertheless the spe- 
cific form of the AL contribution to the transverse conductivity of a layered 
superconductor, which may be considerably suppressed for small intcrlayer 
transparency, suggested to re-examine the contributions from diagrams 5-10 
of Figj^ which are indeed less divergent in e, but turn out to be of lower or- 
der in the transmittance and of the opposite sign with respect to the AL one 
[ p5p6] |. These, so-called DOS, diagrams describe the changes in the normal 
Drude-type conductivity due to fluctuation renormalization of the normal 
quasiparticles density of states above the transition temperature (see Section 
8.5). In the dirty limit, the calculation of contributions to the long itud inal 
fluctuation conductivity a xx from such diagrams was discussed in |10S,99f. 
Contrary to the case of the AL contribution, the in-plane and out-of-planc 
components of the DOS contribution differ only in the square of the ratio of 
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effective Fermi velocities in the parallel and perpendicular directions. This 
allows us to calculate both components simultaneously. The contribution to 
the fluctuation conductivity due to diagram 5 is 

Q^(w„) = 2e 2 Tj2 J 4^s L (<l, a)T^A 2 (q, £ „,a-£„) x 

O k En 

and diagram 6 gives an identical contribution. Evaluation of the integrations 
over the in-plane momenta p and the summation over the internal frequencies 
e n are straightforward. Treatment of the other internal frequencies fl k is less 
obvious, but in order to obtain the leading singular behavior in the vicinity 
of transition it suffices to set fl k — |Q. After integration over q z , we have 
[110: 

5+6__^_ A f 1 

e 2 Ki A , / 2 



A <*P 111 Ml , / , M/2 . ( 189 ) 



8s ap \e U2 + (e + r) 1 / 2 J ' 
where A xx = A yy = l,A zz = (sJ / 'v F ) 2 , A a ^p = and 



2(vft) 2 



1 2 4ttTt / AitTt * \ 2 



In order to cut off the ultra-violet divergence in q we have introduced here a 
cut off parameter g max = = V^/ 2 m com plete agreement with Section 3. 
Let us stress that in the framework of the phenomenological GL theory we 
attributed this cut off to the breakdown of the GL approach at momenta as 
large as q ~ £ . The microscopic approach developed here permits to see 
how this cut off appears: the divergent shortwave-length contribution arising 
from GL-likc fluctuation propagators is automatically restricted by the q- 
dependencies of the impurity vertices and Green functions, which appear at 
the scale q ~ l^ 1 . 

In a similar manner, the equal contributions from diagrams 7 and 8 sum 

to 

Tre 2 f d 2 q 1 



7+8 ne a f 



'M<*»- ( 2lT ) 2 [(e + ^(2)q 2 )(e+^ + ^(2)q 2 )] 1/2 
e 2 n 2 . , / 2 



8s ^ ln Ui/2 + (e+r )i/ 2 )> ( 19 °) 



{v F rf (I 
K2 =2^TV^ 2 
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Comparing (18E) and ( |190| ), we see that in the clean limit, the main contri- 
butions from the DOS fluctuations arise from diagrams 5 and 6. In the dirty 
limit, diagrams 7 and 8 are also important, having -1/3 the value of diagrams 
5 and 6, for both a xx and a zz . Diagrams 9 and 10 are not singular in e << 1 
at all and can be neglected. The total DOS contribution to the in-plane and 
c-axis conductivity is therefore 



<3° 5 = --<Tr)A ali In ( ^ ■ - ■ ) , (191) 



where 

k(Tt) = Kl + K 2 



< 3 [V(§ + 4^)-^(l)-^'(l)] 

56C(3)/tt 4 « 0.691, Tr<l 
8tt 2 (Tt) 2 / [7C(3)] « 9.384 (TV) 2 , 1 « Tt « l/^fe 



is a function of tT only. As it will be shown below at the upper limit 
Tt ~ l/\/e the DOS contribution reaches the value of the other fluctuation 
contributions and in the limit of Tt — > oo exactly eliminates the Maki- 
Thompson one. 



Maki-Thompson contribution We now consider another quantum cor- 
rection to fluctuation conductivity which is called the Maki-Thompson (MT) 
contribution (diagram 2 of Fig. ||). It was firstly discussed by Maki H in 
a paper which appeared almost simultaneously with the paper of Aslama- 
zov and Larkin [5|. Both these articles gave rise to the microscopic theory 
of fluctuations in superconductor. Maki found that, in spite of the seeming 
weaker singularity of diagram 2 with respect to the AL one (it contains one 
propagator only, while the AL one contains two of them) it can contribute 
to conductivity comparably or even stronger than AL one. 

Since the moment of its discovery the MT contribution became the sub- 
ject of intense controversy. In its original paper Maki found that in 3D case 
this fluctuation correction is four times larger than the AL one. In 2D case 
the result was striking: the MT contribution simply diverged. This paradox 
was, at least at the level of recipe, resolved by Thompson 0: he proposed 
to cut off the infra-red divergence in the Cooper pair center of mass mo- 
mentum integration by the introduction of the finite length l s of inelastic 
scatterings of electrons on paramagnetic impurities. In the further papers of 



Patton 110 , Keller and Korenman [ 1 1 1 1 it was cleared up that the presence 
of paramagnetic impurities or other external phase-breaking sources is not 
necessary: the fluctuation Cooper pairing of two electrons results in a change 
of the quasiparticle phase itself and the corresponding phase-breaking time 
t v appears as a natural cut off parameter of the MT divergence in the strictly 
2D case. The minimal quasi-two-dimensionality of the electron spectrum, as 
we will show below, automatically results in a cut off of the MT divergence. 
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Although the MT contribution to the in-plane conductivity is expected to 
be important in the case of low pair-breaking, experiments on high-temperature 
superconductors have shown that the excess in-plane conductivity can usually 
be explained in terms of the fluctuation paraconductivity alone. Two possible 
explanations can be found for this fact. The first one is that the pair-breaking 
in these materials is not weak. The second is related with the d— wave symme- 



try of pairing which kills the anomalous Maki-Thompson process 112, 113 1. 
We will consider below the case of s-pairing, where the Maki-Thompson pro- 
cess is well pronounced. 

The appearance of the anomalously large MT contribution is nontrivial 
and worth being discussed. We consider the scattering lifetime r and the pair- 
breaking lifetime to be arbitrary, but satisfying t v > r. In accordance with 
diagram 2 of Fig. (||) the analytical expression for the MT contribution to the 
electromagnetic response tensor can be written as 



Q%?(u v ) = 2e 2 Tj2 J ■^L(q,n k )l a p(ti,fi k ,uu), (193) 



d 3 q 

where 

7 Q /3(q, fl k , w u ) — n—v 

)\(q,e n ,!l k _ n ) x (194) 



f d 3 p 

In the vicinity of T c , it is possible to restrict consideration to the static 



limit of the MT diagram, simply by setting fl k = in (193). Although dy- 
namic effects can be important for the longitudinal fluctuation conductivity 
well above Tld, the static approximation is correct very close to T c , as shown 



in |)8 121 ]. The main q— dependence in ( |193| ) arises from the propagator and 



vertices A. This is why we can assume q — in Green functions and to calcu- 
late the electron momentum integral passing, as usual, to a £(p) integration: 

l a p(q,0,u u ) = nu(v a (p)vp(q-p)) FS x (195) 

1 1 1 



(|2 £ „ + „| + iv) (|2 £n | + i3<r 



2 e 



n-\-v | ~r \c-n \ 



In evaluating the sum over the Matsubara frequencies e n in (195) it is useful 
to split it into the two parts. In the first e n belongs to the domains ] — oo, — u> u [ 
and [0, oo[, which finally give two equal contributions. This gives rise to the 
regular part of the MT diagram. The second, anomalous, part of the MT 
diagram arises from the summation over e„ in the domain [— w„,0[. In this 
interval, the further analytic continuation over uj u leads to the appearance of 
an additional diffusive pole: 
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I a p(q,Q,U») = lap fa, "is) + lap' fa, "is) = v ( v a (p)vp fa - p)) Fs 

1 1 1 



Aan), 



(2e n+v + 5g 2 ) (2e„ + 5g 2 ) 



2e n + cj„ + t 



1 



1 



71 — — U I ^C n _|_ 



The limits ol summation in the first sum do not depend on uj v , so it is an 
analytic function of this argument and can be continued to the upper half- 
plane of the complex frequency by the simple substitution oj v — > —ito. Then, 
tending lj — > 0, one can expand the sum over powers of u> and perform the 
summation in terms of digamma-function: 



J tp fa^) = v(v a (p)v (q-p)) FS x 

d 




(196) 



The values of characteristic momenta q <C I 1 are determined by the do- 
main of convergibility of the final integral of the propagator L(q, 0) in (193) 
(analogously to (19C) and one can neglect Dq 2 with respect to t _1 . In result 



J i7 s)J? (^ w °) = v K(p)M<i - P)> FS X 



(197) 



const 



1 2 AtttT 



r (j) 

AttTt 



The appearance of the constant in Q a p(co^) was already discussed in the case 
of the AL contribution and, as was mentioned there, it is cancelled with the 
similar contributions of the other diagrams |9§[ ] and we will not consider it 
any more. 

Now let us pass to the calculation of ijffi fa, oj v ). Expanding the summing 
function in simple fractions one can express the result of summation in terms 
of digamma-functions 



(198) 
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Doing the analytical continuation icj v — > u> — > and taking into account that 
in the further q-integration of I^^ R {q, lj — > 0), due to the singular at small 
q propagator, the important range is Dq 2 <C T, one can find 

T (an)R, "TUT v (^a(p)^(q ~ P))fS nn^ 

Because of the considerable difference in the angular averaging of the different 
tensor components we discuss the MT contribution to the in-plane and out 
of plane conductivities separately. 

Taking into account that (v x (p)v x (c[ — p)) FS = —v 2 F /2 one can find that 
the calculation of the regular part of MT diagram to the in-plane conductivity 
is completely similar to the corresponding DOS contribution and here we list 
the final result |)7| only: 

a MT(re 9 ) = _^~ ln 



where 
k(Tt) 



2s Ve 1/2 + (e + r) 1 ? 2 , 

- 2 [* (i + 4^t) - * (i) - £t# m 

28C(3)/tt 4 ps 0.346, for Tt < 1 
tt 2 / [14C(3)] m 0586, for K Tt < 



(200) 



is a function only of tT. We note that this regular MT term is negative, as 
is the overall DOS contribution. 

For the anomalous part of the in-plane MT contribution we have: 



MT(<m) 

XX 



d 3 q 1 




4s(e - 7 P ) 



?7(2)q 2 + |(1 ~ cosq z s)] 

(201) 



where, in accordance with , the infra-red divergence for the purely 2D case 
(r = 0) is cut off at Dq 2 ~ 1/ rj^|. The dimensionless parameter 

2?7 tt f 1, Tt < 1 

7v ~ 4^ ~" "8TV^ t 7C(3)/ (2^ 3 Tr) , 1«Tt « l/y/l 



22 



The detailed study of the phase-breaking t ime , its energy dependence and the 



effect on the MT contribution was done in [115] 
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is introduced for simplicity. If r 7^ the MT contribution turns out to be 



hnite even with t v — oo. Comparison of the expressions ( |187| ) and ([201]) indi- 
cates that in the weak pair-breaking limit, the MT diagram makes an impor- 
tant contribution to the longitudinal fluctuation conductivity: it is four times 
larger than the AL contribution in the 3D regime, and even logarithmically 
exceeds it in the 2D regime above Tld. For finite pair-breaking, however, the 
MT contribution is greatly reduced in magnitude. 

We now consider the calculation of the MT contribution to the transverse 



conductivity. The explicit expressions for v z (p) and v z (q—p) (see Exp. ( 181 )), 
result in (v x (p)v x (q — p)) FS = \ J 2 s 2 cos q z s. We take the limit Jt << 1 in 
evaluating the remaining integrals, which may then be performed exactly. 



The regular part of the MT contribution to the transverse conductivity is 
e 2 s 2 irrk(TT) f d 3 q cosq z s 



MT(rcg) 



(27T) 



?7(2)q z 



§(1 - cosg z s) 



e 2 srk{TT ) ( (e + r) 1 / 2 - e 1 / 2 
16r?( 2 ) 



„l/2 



This term is smaller in magnitude than is the DOS one, and therefore makes 
a relatively small contribution to the overall fluctuation conductivity. In the 
3D regime below Tld, it is proportional to J 2 , and in the 2D regime above 
Tljj, it is proportional to J 4 . 

For the anomalous part of the MT diagram one can find 



MT(a 



111 ! 

ire J s t 



ire 2 s 



d 3 q 



cos q z s 



(2?r) 3 [1/ Tv + Dg 2 ][e + ry (2) q 2 + §(1 - cosq z s) 



d 2 q 



l v + ?7(2)Q 2 +r/2 



[h v + V(2)<l 2 )(l v + ?7(2)q 2 + r)] 1 / 2 



e + ^(2)q 2 +r/2 



[(e + ^(2)q 2 )(e + ??(2)q 2 +^)] 1/2 
7^ + r + e 



16^) V W + r)Y/ 2 + [ lv {l v +r)]^ 



(202) 
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In examining the limiting cases of ( p02| ), it is useful to consider the cases of 
weak (j v « r, ^=> J 2 tt v » 1/2) and strong (j v >> r,<=> J 2 tt v « 
1/2) pair-breaking separately^. For weak pair- breaking, we have 



MT(an) c a I f~~T~ 



16V W [r/(2e) 7 7v «r«e 

In this case, there is the usual 3D to 2D dimensional crossover in the anoma- 
lous MT contribution at Tld- There is an additional crossover at T v (where 
T c < T v < Tle>), characterized by e(T ¥ ,) = 7 V , below which the anomalous 
MT term saturates. Below Ild, the MT contribution is proportional to J, 
but in the 2D regime above Tld, it is proportional to J 2 . 
For strong pair-breaking 



IT 



MT(an) £ '' S f r/^, £ < T < 7 V 



32^(2) \ r 2 / (4 7¥ ,e) , r<min{7 v ,e} 

In this case, the 3D regime (below Tld) is not singular, and the anomalous 
MT contribution is proportional to J 2 , rather than J for weak pair-breaking. 
In the 2D regime, it is proportional to J 4 for strong pair-breaking, as op- 
posed to J 2 for weak pair-breaking. In addition, the overall magnitude of 
the anomalous MT contribution with strong pair-breaking is greatly reduced 
from that for weak pair-breaking. 

Let us now compare the regular and anomalous MT contributions. Since 
these contributions are opposite in sign, it is important to determine which 
will dominate. For the in-plane resistivity, the situation is straightforward: 
the anomalous part always dominates over the regular and the latter can be 
neglected. The case of c— axis resistivity requires more discussion. Since we 
expect t v > r, strong pair-breaking is likely in the dirty limit. When the 
pair-breaking is weak, the anomalous term is always of lower order in J than 
the regular term, so the regular term can be neglected. This is true for both 
the clean and dirty limits. The most important regime for the regular MT 
term is the dirty limit with strong pair-breaking. In this case, when t v T ~ 1, 
the regular and anomalous terms are comparable in magnitude. In short, it 
is usually a good approximation to neglect the regular term, except in the 
dirty limit with relatively strong pair-breaking and only for the out-of-plane 
conductivity. 

Finally let us mention that the contributions from the two other diagrams 
of the MT type (diagrams 3 and 4 of Fig. ||) in the vicinity of critical temper- 
ature can be omitted: one can check that they have an additional square of 

23 

Phys ically the value J 2 r. characterizes the effective interlayer tunneling rate 
[ [)7 1 14 1 . When l/r^ << J 2 r << 1/r, the quasiparticles scatter many times 



before tunneling to the neighboring layers, and the pairs live long enough for 
them to tunnel coherently. When J 2 r << l/r^>, the pairs decay before both 
paired quasiparticles tunnel. 
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Fig. 11. The normalised excess conductivity for samples of YBCO-123 (triangles), 
BSSCO-2212(squares) and BSSCO-2223 (circles) plotted against e = In T/T c on a 
ln-ln plot as described in [124]. The dotted and solid lines are t he AL theory in 3D 
and 2D respectively. The dashed line is the extended theory of [125 .. 



the Cooper pair center of mass momentum q in the integrand of ^-integration 
with respect to diagram 2 and hence turn out to be less singular in e. 



Discussion Although the in-plane and out-of plane components of the fluc- 
tuation conductivity tensor of a layered superconductor contain the same 
fluctuation contributions, their temperature behavior may be qualitatively 
different. In fact, for a[ l x , the negative contributions are considerably less than 
the positive ones in the entire experimentally accessible temperature range 
above the transition, and it is a positive monotonic function of the temper- 
ature. Moreover, for HTS compounds, where the pair-breaking is strong and 
the MT contribution is in the overdamped regime, it is almost always enough 
to take into account only the paraconductivity to fit experimental data. Some 
examples of the experimental fi ndings for in-plane fluctuation c onductivity 
of HTS materials on can see in \\L 1 6j|l 1 7| , |T]l| , pT9|Jl^ , [l^ , |i^jT23| . 



In Fig. 



11 



the fluctuation part of in-plane conductivity a[ l x is plotted 
as a function of e = \nT/T c on a double logarithmic scale for three HTS 
samples (the solid line represe nts th e 2D AL behavior (1/e), the dotted line 
represents the 3D one: 3.2/^/e) 124]. One can see that paraconductivity of the 
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less anisotropic YBCO compound asymptotically tends to the 3D behavior 
(1/e 1 / 2 ) for e < 0.1, showing the LD crossover at e w 0.07; the curve for 
more anisotropic 2223 phase of BSCCO starts to bend for e < 0.03 while the 
most anisotropic 2212 phase of BSCCO shows a 2D behavior in the whole 
temperature range investigated. All three compounds show a universal 2D 
temperature behavior above the LD crossover up to the limits of the GL 
region. It is interesting that around e ~ 0.24 all the curves bend down 
and follow the same asymptotic 1/e 3 behavior (dashed line). Finally at the 
value e « 0.45 all the curves fall down indicating the end of the observable 
fluctuation regime. 



Reggiani et al. [125] extended the 2D AL theory to the high temperature 
region by taking into account the short wavelength fluctuations. The following 
universal formula for 2D paraconductivity of a clean 2D superconductor as a 



function of the generalized reduced temperature e = In T/T c was obtained^: 

41 = £/M 

with /(e) = e-\ e < 1 and /(e) = e" 3 , e > 1. 

In the case of the out-of-plane conductivity the situation is quite different. 
Both positive contributions (AL and anomalous MT) are suppressed by the 
interlayer transparency, leading to a competition between positive and nega- 
tive terms. This can lead to a maximum in the c-axis fluctuation resistivity 
which occurs in the 2D regime (in the case discussed Jr << l,r« << 1 and 
7 v k > 1) : 



/ ~ 1 _L 



1 

27^ 



This nontrivial effect of fluctuations on the transverse resistance of a layered 
superconductor allows a successful fit to the data observed on optimally and 
overdoped HTS samples (see, for instance, FigJT^) where the growth of the 
resistance still can be treated as the correction. 

The fluctuation mechanism of the growth of the transverse resistance can 
be easily understood in a qualitative manner. Indeed to modify the in-plane 
result ( |100| ) for the case of c-axis paraconductivity one has to take into ac- 
count the hopping character of the electronic motion in this direction. If the 
probability of one-electron interlayer hopping is V\, then the probability of 
coherent hopping for two electrons during the fluctuation Cooper pair lifetime 

tql is the conditional probability of these two events: Vi = V\i^P\TGh)- The 

?2 1 
1 T 5 '' 



transverse paraconductivity may thus be estimated 



in complete accordance with (18S). We see that the temperature singularity 



of <r^ L turns out to be stronger than that in &f L , however for a strongly 

24 In Section 7.2 we will demonstrate how such a dependence (l/ln 3 (T'/T' c )) ap- 
pears by accounting for short wavelength fluctuations for the 2D fluctuation 
susceptibility. 
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Fig. 12. Fit of the temperature dependence of the transverse resistance of an 
unde rdoped BSCCO c-axis oriented film with the results of the fluctuation theory 
126]. The inset shows the details of the fit in the temperature range between T c 
and 110K. 



anisotropic layered superconductor cr^ is considerably suppressed by the 
square of the small probability of inter-plane electron hopping which enters 
in the pre-factor. It is this suppression which leads to the necessity of taking 
into account the DOS contribution to the transverse conductivity. The latter 
is less singular in temperature but, in contrast to the paraconductivity, man- 
ifests itself in the first, not the second, order in the interlayer transparency 
a® os ~ — V\ In -. The DOS fluctuation correction to the one-electron trans- 
verse conductivity is negative and, being proportional to the first order of Pi, 
can completely change the traditional picture of fluctuations just rounding 
the resistivity temperature dependence around transition. The shape of the 
temperature dependence of the transverse resistance mainly is determined 
by competition of the opposite sign contributions: the paraconductivity and 
MT term, which are strongly temperature dependent but are suppressed by 
the square of the barrier transparency and the DOS contribution which has 
a weaker temperature dependence but depends only linearly on the barrier 
transparency. 



6 Manifestation of fluctuations in various properties 

In this Section we will demonstrate the applications of the microscopic theory 
of fluctuations. The limited volume does not permit us to deliver here the 
systematic review of the modern theory and we restrict ourselves only by 
presentation of the several representative recent studies. The details one can 
find in the articles cited in subtitles. 

It is necessary to underline that the comparison of the results of fluc- 
tuation theory with the experimental findings on HTS materials has to be 
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considered sooner in qualitative than quantitative context. Indeed, as is clear 
now, the superconductivity in the most of HTS compounds has the non- 
trivial symmetry. Moreover, as was discussed in the previous Section, these 
compounds are rather clean than dirty. Both these complications can be taken 

,p7[) but this was not done in the majority 
of the cited papers. 



into account (see for example [112 



6.1 T he effects of fluctuations on magnetoconductivity 
1 97j ^27|Jl5^JT55 1 



The experimental investigations of the fluctuation magnetoconductivity are 
of special interest first because this physical value weakly depends on the nor- 
mal state properties of superconductor and second due to its special sensitiv- 
ity to temperature and magnetic field. The role of AL contribution for both 
the in-plane and out-of-plane magnetoconductivities was studied above in the 
framework of the phenomenological approach. The microscopic calculations of 
the other fluctuation corrections to the in-plane magnetoconductivity conduc- 
tivity show that the MT contribution has the same positive sign and temper- 
ature singularity as the AL one. In the case of weak pair-breaking it can even 
considerably exceed the latter. The negative DOS contribution, like in the 
case of the zero-field conductivity, turns out to be considerably less singular 
and many authors (see e.g. Refs. [L36lp^ , p^ , |l39Ml42[[i40| , 



141 



143,144,145|, |146,147|) 



successfully explained the in-plane magnetoresistance data in HTS using the 
AL and MT contributions only p4^ , pT9| , ^50p5ll ] . 

Turning to the out-of-plane magnetoconductivity of a layered supercon- 
ductor one can find a quite different situation. Both the AL and MT contri- 
butions turn out to be here of the second order in the interlayer transparency 
and this circumstance makes the less singular DOS contribution, which re- 
mains however of first order in transparency, to be competitive with the main 
terms [127]. The large number of microscopic characteristics involved in this 



competition, like the Fermi velocity, interlayer transparency, phase-breaking 
and elastic relaxation times, gives rise to the possibility of occurrence of dif- 
ferent scenarios for various compounds. The c-axis magnetoresistance of a set 
of HTS materials shows a very characteristic behavior above T c q. In contrast 
to the ab-plane magnetoresistance which is positive at all temperatures, the 
magnetoresistance along the c-axis has been found in many HTS compounds 
(BSSCO @|3^Q^| , LSSCO [H|, YBCO @ and T1BCCO jl3j ) 
to have a negative sign not too close to T c q and turn positive at lower tem- 
peratures. We will show how this behavior find its explanation within the 
fluctuation theory |97]]. 

We consider here the effect of a magnetic field parallel to the c-axis. In this 
case both quasiparticles and Cooper pairs move along Landau orbits within 
the layers. The c-axis dispersion remains unchanged from the zero-field form. 
In the chosen geometry one can generalize the zero-field results reported in 
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the previous Section to finite field strengths simply by the replacement of the 
two-dimensional integration over q by a summation over the Landau levels 



d 2 q H 



t-^t 2irr\,^ 



(let us remind that rjp) — £ 2 y )- ^° * ne general expressions for all fluctuation 
corrections to the c-axis conductivity in a magnetic field can be simply written 
in the form [ p7| : 

^ L = S?E " JT2 ( 203 ) 

{[e + H^n + l)][r + e + h(2n + 1)]} 3/2 



n .ii"s _ e 2 srnh 



1/h 



zz 



a 



^ {[e + h(2n + l)][r + e + fe(2n + 1)]} 1/2 

MT(re 9 ) = y» /" e + fe(2n + 1) + r/2 \ 

,to V{[ £ + M2n+l)][r + e + / l (2n+l)]} 1/2 J 



°zz 



MT(an) = e * sh ( j v + h{2n + 1) + r /2 

8&,(e " 7*0 r to \j[( 7v + fc(2n + l)][ 7v + fc(2n + 1) + r)]} 1/2 



e + /i(2n + l)+r/2 



{[e + h(2n + l))[r + e + h(2n + 1)]} 1/2 



(206) 



For the in-plane component of the fluctuation conductivity tensor the only 
additional problem appears in the AL diagram, where the matrix elements 
of the harmonic oscillator type, originating from the B» (qi\ ) blocks, have to 
be calculated. How to do this was demonstrated in details in Section 4. The 
other contributions are essentially analogous to their c-axis counterparts: 



(207) 



4s ,to \{[e + h(2n + l)][r + e + h(2n + 1)]} 1/2 

2 

{ [e + h(2n + 2)] [r + e + h(2n + 2)]} 1/2 

! ) 

{ [e + h(2n + 3)] [r + e + h(2n + 3)]} 1/2 / ' 

<t£° s +*£ t ™ = - e2h{K + k) V 1 (208) 

2* ^ {[e + h(2n + l)][r + e + h(2n + 1)]} 1/2 
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and 



MT(an) 



e 2 h 



4s(e 



Tip) S 



{[( lv + h(2n + l)][ 7v + h(2n + 1) + r)]} 



1/2 



{[e + h(2n + l)][r + e + fr(2n + 1)]} 



1/2 



(209) 



These results can in principle be already used for numerical evaluations and 
fitting of the experimental data which was indeed successfully done in a series 
of papers j|l 34| , |l 5 2jjl35|| . 

The detailed comparison of the cited results with the experimental data 



[134 153], especially in strong fields, raised the problem of regularization of 
the DOS contribution. If in the absence of the magnetic field its ultra-violet 
divergence was successfully cut off at q ~ in the case under consideration 
the cut off parameter depends on the magnetic field and makes the fitting 
procedure ambiguous. The solution of this problem was proposed in [ 154 j, 
where the authors cal culat ed th e diff erence Aa zz os = a zz os (h, e)-a z J z os (0, e) 
applying to formulas (204) and ( ^08[ ) the same trick which was already used in 
Section 2 for the regularization of the free energy in magnetic field (Eq. (|78|)). 
The corresponding asymptotics for all out-of-plane fluctuation contributions 
are presented in the following table: 
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Table 3 



The procedure described gives an excellent fitting up to very high fields 
| fl55| which is shown in Fig . [l3| 

Let us start the analysis from the 2D case (r <C e). One can see that 
here the positive DOS contribution to magnetoconductivity turns out to be 
dominant. It grows as H 2 up to H C 2(e) and then crosses to a slow logarithmic 
asymptote. At H ~ H c2 (0) the value of Aa?° s (h ~ 1, e) = -cr£ os (0,e) 
which means the total suppression of the fluctuation correction in such a 
strong field. The regular part of the Maki-Thompson contribution does not 
manifest itself in this case while the AL term can compete with the DOS one 
in the immediate vicinity of T c , where the small anisotropy factor r can be 
compensated by the additional e 3 in the denominator. The anomalous MT 
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Fig. 13. Magnetoconductivity versus temperature at 27 T for an underdoped Bi- 
2212 single crystal. The solid line represents the theoretical c alcu lation. The symbols 
are the experimental magnetoconductivity Aa zz (B i c i /) [155 



contribution can contribute in the case of small pairbreaking only, which is 
opposite to what is expected in HTS. 

In the 3-D case (e <C r) the behavior of the magnetoconductivity is more 
complex. In weak and intermediate fields the main, negative, contribution 
to the magnetoconductivity occurs from the AL and MT terms. At H ~ 
H C 2(e)(h ~ e) the paraconductivity is already considerably suppressed by the 
magnetic field and the h 2 — dependence of the magnetoconductivity changes 
through the y/j- tendency to the high field asymptote — aii L \o, e). In this 
intermediate region of fields (e <C h <C r), side by side with the decrease (~ 
v/1T) °f t ne main AL and MT contributions, the growth of the still relatively 
small DOS term takes place. At the upper limit of this region (h ~ r) its 
positive contribution is of the same order as the AL one and at high fields 
(r -C h <C 1) the DOS contribution determines the slow logarithmic decay of 
the fluctuation correction to the conductivity which is completely suppressed 
only at H ~ H C 2(0). The regular part of the Maki- Thompson contribution is 
not of special importance in the 3D case. It remains comparable with the DOS 
contribution in the dirty case at fields h < r, but decreases rapidly (~ f ) at 
strong fields ( h > r), in the only region where the robust Aa^ os (h, e) ~ In - 
shows up surviving up to h ~ 1 . 

The temperature dependence of the different fluctuation contributions to 
the magnetoconductivity calculated for an underdoped Bi-2212 single crystal 
at the magnetic field 27 T is presented in Fig. [l4]. 

The formulas for the in-plane magnetoconductivity are presented in the 
table: 
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Fig. 14. Decomposition of the calculation of total theoretical magnetoconductivity 
for an underdoped Bi-2212 single crystal at 27 T. The inset shows the regular and 
anomalous parts of the MT contribution which are t oo small to be presented in 
the same scale as the AL and DOS contributions [1531. 
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Table 4 



Analyzing it one can see that in almost all regions the negative AL and 
MT contributions govern the behavior of in-plane magnetoconductivity. Nev- 
ertheless, similar to the c-axis case, the high field behavior is again determined 
by the positive logarithmic A(a x J x 0S + a x % T< " re9 ' > ) contribution, which is the 
only one to survive in strong field. It is important to stress that the sup- 
pression of the DOS contribution by a magnetic field takes place very slowly. 
Such robustness with respect to the magnetic field is of the same physical 
origin as the slow logarithmic dependence of the DOS-type corrections on 
temperature. 

Another important problem which appears in the fitting of the resistive 
transition shape in relatively strong fields with the fluctuation theory is the 
much larger broadening of the transition than predicted by the Abrikosov- 
Gorkov theory [156|. Kim and Gray |Q explained the broadening of the 
c-axis peak with increasing magnetic field in terms of Josephson coupling, de- 
scribing a layered superconductor as a stack of Josephson junctions. In Rcfs. 
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[157 7^| the self-consistent Hartree approach was proposed for the extension 
of fluctuation theory beyond the Gaussian approximation. It results in the 
considerable shift of T C (H) toward low temperatures with a corresponding 
broadening of the transition. The renormalized reduced temperature Eh is 
determined according to the self-consistent equation (7^[ : 

l/h 



£h — £h 



1 



Gi w h 



-! [(e h + hn){e h + h(n + 1) + r)}^ 



(210) 



71=0 



The authors of |15J|, followin g the p roce dure proposed b y Do rsey and Ullah 
\ \H\ , modified the expressions (203)- (209) by account for ( 210 ). In result they 
succeeded to fit quantitatively both in-plane resistivity transition and the 
transverse resistivity peak for BSCCO films strongly broadened by applied 
magnetic field. 



6.2 Fluctuations far from T c or in strong magnetic fields 

As was mentioned above the role of fluctuations is especially pronounced in 
the vicinity of the critical temperature. Nevertheless for some phenomena 
they can be still considerable far from the transition too. In these cases the 
GL theory is certainly unapplicable since the short-wave and dynamical fluc- 
tuation contributions have to be taken into account. It can be done in the 
microscopic approach which we will demonstrate in several examples. 



Fluctuation magnetic susceptibility far from transition [29]. The 

given above qualitative estimations (|58|)-(|64|) for the fluctuation diamagnetic 
susceptibility based on the Langevin formula, demonstrate that even at high 
temperatures T T c it turns to be of the order of xp for clean 3D supercon- 
ductors and exceeds noticeably this value for 2D systems. In order to develop 



the microscopic theory [159,29,161 let us start from the general expression 



for free energy in the one-loop approximation: 

F = T^Tr{ln[l- ff iT(a,r,r')]}, (211) 

where g is the effective interaction constant related with the transition critical 
temperature by ( |158| ). This approximation corresponds to the ladder one (see 
(150)) for the fluctuation propagator. The polarization operator i7(J?fe,r, r') 
is determined by expression (151) but in the case of an applied magnetic field 
the homogeneity of the system is lost and 7T(f2fc,r,r') depends not on the 
space variable difference r— r' but on each separately. Expanding 77(J7i c , r, r') 
one can express the magnetic susceptibility of a layered superconductor in a 
weak magnetic field perpendicular to the layers in terms of the derivatives 
n x = ^TT(q) ||: 

d 2 F 2_ 2m ^ f d 3 q r3j 



x = 'dm = 3 J j^ L n ^ n * n yy ~ n y n ^- ^ 
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The final expressions for the fluctuation diamagnetic susceptibility in the 
clean and dirty cases for wide range of temperatures can be written as: 

v O) m _ XP /0.05(ln- 2 (^/T c ))-ln- 2 (T/T c )), r" 1 « T « uj d 

X f l [ ' ~ 3 1 2VTVI11- 1 (T/T c ), T c « T « r- 1 [ L6) 

(2) 0.05 (E F \ 1 

Let us stress that these results are valid for the fluctuation diamagnetism 
of a normal metal with g > too, if by T c one uses the formal value T c ~ 
^exp(J-). 



Fluctuation magnetoconductivity far from transition [162 . Let us 

discuss the conductivity of the 2D electron system with impurities in a mag- 
netic field at low temperatures. Even in the absence of the field the effects 
of quantum interference of the non-interacting electrons in their scatterings 
on elastic impurities already results in the appearance of a nontrivial tem- 
perature dependence of the resistance. This result contradicts the statement 
of the classical theory of metals requiring the saturation of the resistance at 
its residual value at low temperatures. In a superconductor above the criti- 
cal temperature this, so-called weak localization (WL), effect is amplified by 
the Andreev reflection of electrons on the fluctuation Cooper pair leading 
to appearance of the MT correction to the conductivity. The characteristic 
feature of both the MT and WL corrections is their extreme sensitivity to 
the dephasing time t v and to weak magnetic fields. 

Beyond the GL region (T > T c ) the MT correction is determined by 
the same diagram 2 of Fig. ^| but now the dynamic (Qk 7^ 0) and short- 
wave-length (q <~*j fluctuation modes have to be taken into account. The 



corresponding calculations were performed in |98,162| and the result can be 
written in the form: 

Sa WL+ MT = ^ [« - / 3 ( T )] Y (^ T£ ) . ( 215 ) 

where we introduced the effective Larmour frequency for the diffusion motion 
Ql = 4DeH with the diffusion coefficient T) W\ and the function 

1 1 f T « 1 

Y(x)=]nx + ip{- + -) = i 24' X< ^\ . (216) 

The first term in this formula corresponds to the WL contribution (a = 1 if 
the spin-orbit interaction of the electrons with the impurities is small while 

25 Comparison of the expressions (^), ( |ll7| ) and ( |215[ ) relates the Larmour fre- 
quency with the dimensionless field: h = Ol/2T c introduced in section 2 and the 
diffusion coefficient with the phenomenological GL constants T> = 1/ma. 



Fluctuation Phenomena in Superconductors 93 



in the opposite limiting case a — —1/2), the second describes the MT con- 
tribution to magnetoconductivity. The function /3[ln(T/T c )] was introduced 
in Ref.p62|, At T -> T c f3(x) = l/x and ( |2"l~5| ) reduces to the already studied 
MT correction in the vicinity of critical temperature. For T ^> T c j3(x) = l/x 2 
and the MT contribution gives a logarithmically small correction to the WL 
result. Its zero-field value, being proportional to \aT 2 {T /T c ), decreases with 
the growth of the temperature faster than both the AL contribution (in the 
dirty case 5a al ~ 1/ \n(T/T c )) and the especially slow DOS contribution 
(Sctdos ~ lnln(l/T c r) - lnln(T/T c )) (see Ref. |f||). 

It worth ment io ning that for the regio n of temperatures T 3> T c , analo- 
gous to Exp.( f2~13| )- (214), the result ( |215|) can be applied both to supercon- 
ducting and normal metals {g > 0), if in place of the critical temperature 
the formal value T c ~ Ep exp(-^) is undermined. The interplay of the local- 
ization and fluctuation corrections was extensively studied (see, for example, 



[163 164,165,166 ; 167l]) 



Fluctuations in magnetic fields near H C 2(0) [168]. As one can see from 
(207)- (209), in the vicinity of the upper critical field H C 2(T) the fluctuation 

corrections diverge as 1 for the 2D case and as for the 313 case ^| (it is 
enough to keep just the terms with n = in these formulas). This behavior is 
preserved in strong magnetic fields too, but the coefficients undergo changes. 
A case of special interest is T <C T c (which means H — > H C 2(0)) which repre- 
sents an example of a quantum phase transition [16S|. Microscopic analysis 
of the magnetoconductivity permits us to study the effect of fluctuations in 
magnetic fields of the order of H C 2(0), where the GL functional approach is 
inapplicable. 

We restrict our consideration to the case of a dirty metal (TV <C 1). In 
this limit |w n +^ — £>_ n | ~ r _1 and the Green function correlator ( |153| ) can 
be written in the form 



\si - s 2 \ 



Vq' 



(217) 



Expressing 7T(q, S7k) in terms of V{c\, £i, £2) by means of ( 172 ) and using the 
definition of the critical temperature one can find an explicit formula for the 
fluctuation propagator 




(218) 



The prominent characteristic of this expression is that it is valid even rela- 
tively far from the critical temperature (for temperatures T min{r _1 , lud}) 
and for |q| <C Z~\ <C wj). 

26 eh is the renormalized by the magnetic field reduced temperature = e + h 
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One can rewrite this expression in a magnetic field applied along the c-axis 
in the Landau representation by simply replacing (vq 2 ^ => f^h{n + 1/2) 

|59| : 



T fl 

ln 7; + n2 



4ttT 

72 „;„2 , 



(219) 



f2 L (n + 1/2) + At J 2 sin 2 (g z s/2) \ _ fl 
AttT ) ^ { 2 



In the case of arbitrary temperatures and magnetic fields the expression 
for the AL contribution to the conductivity takes the form: 



Vxx(u v ) = ve 2 T^^ ^2 B n>m (f2 k +uj v , (1 k )L m (f2 k ) x 

£2}- n,m— 



(220) 



(we have restricted our consideration to the 2D case). The expression for 
Bn,m(dk,u u ) can be rewritten as 



2-kv 



m,n+l m,n— 1 J 



with 



A m (ei, £2) 



\ n {Si + Qk — £i)A m (£i, ftk — Ei) 



e(-eiea) 



t|£i-£ 2 | + ^l(to + 1/2) 



(221) 



(222) 



The critical field H c2 {T) is determined by the equation \q z = 0, Qk = 
0) = 0. This is why in the vicinity of H C 2(T) the singular contribution to ( |220p 
originates only from the terms with Lo(0, Ok). The frequency dependencies 
of the functions B ritm (f2k ■!?&) and Li(f2k) are weak although we cannot 
omit them to get nonvanishing answer. It is enough to restrict ourselves to the 
linear approximation in their freq uency dependencies. If the temperature T <C 
T c o the sum over frequencies in (221) can be approximated by an integral. 
Transforming the boson frequency Qk summation to a contour integration 
as was done above and making the analytic continuation in the external 
frequency lo u one can get an explicit expression for the d.c. paraconductivity. 
In the same spirit the contributions of all other diagrams from the Fig. || 
which contribute to fluctuation conductivity in the case under discussion are 
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calculated side by side with the AL one. The final answer can be presented 
in the form: 



5a t0 



2e 2 


( ln^ Tc0 


3tt 2 


+"0 


' T cQ fH 





3 7 eT ( H c2 (T) 



H c2 {T) 



H-H c2 (T) 



( T d 



H c2 (T) 
H c2 (T) 



(223) 



\2 lE T H - H c2 {T) 



1 cO 



H-H c2 (T) 
2 lE T V H c2 (T) 



- 1 



where j E is the Euler constant. Let us consider some limiting cases. If the 
temperature is relatively high T/T c q ^> (H — H c2 (T)) /H c2 (T), we obtain 
the following formula for the fluctuation conductivity: 



, 2 7£ e 2 ... 

0(7 = s — — — 



T 



H c2 (T) 
H-H c2 (T) 



(224) 



If H < H c2 (0), we can introduce T C (H) and rewrite Eq.(|22J) in the usual 
way 



5a 



3e 2 



2 1e tt 2 T-T C {H) 



If H > H c2 (0), in the low-temperature limit T/T c0 < (H 
we have 



5a 



2e 2 



In- 



H c2 {T) 



3tt 2 H-H c2 (T) 



(225) 

H c2 (T))/H c2 (T) 
(226) 



One can see, that even at zero temperature a logarithmic singularity remains 
and the corresponding correction is negative. It results from all three fluctu- 
ation contributions, although the DOS one exceeds the others by numerical 
factor. Let us recall that in the case of the c-axis conductivity of a layered su- 
perconductor, or in granular superconductors abov e T c , the DOS contribution 
exceeds the MT and AL ones parametrically [ 160 ] . 



6.3 The effect of fluctuations on the Hall conductivity [ 174 1 

Let us start with a discussion of the physical meaning of the Hall resistivity 
p xy . In the case of only one type of carriers it depends on their concentration 
n and turns out to be independent of the electron diffusion coefficient: p xy — 
H/ (en) . The fluctuation processes of the MT and DOS types contribute to 
the diffusion coefficient, so their expected contribution to the Hall resistivity 
is zero. For the Hall conductivity in a weak field one can write 

a xy = Pxy a 2 xx = Pxy a x ^ 2 + 2p xy a x n j5a xx = a<$ (l + 2^) (227) 
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so, evidently, the relative fluctuation correction to Hall conductivity is twice 
as large as the fluctuation correction to the diagonal component. This qual- 
itative speculation is confirmed by the direct calculation of the MT type 
diagram 16£]. 

The AL process corresponds to an independent charge transfer which 
cannot be reduced to a renormalization of the diffusion coefficient, ft con- 
tributes weakly to the Hall effect, and this contribution is related to the 
Cooper pa ir particle-hole a sy mmetry. This effect was investigated in a set 
of papers: 16£ , 170, 171 172 71, 173, 174 1. Let us recall that the proper general 
expression describing the paraconductivity contribution to the Hall conduc- 
tivity in the general case of arbitrary magnetic fields and frequencies (in the 
TDGL theory limit s) w as already carried out above in the phenomenological 
approach (see Eq. (|113|)). The microscopic cons ider ation of this value can be 
done in the spirit of the calculation of o^x ( see ( |185| )) and after the analytical 
continuation results in 



„.-.! ._(?M9»V £(» + !) f'-p. 



coth — dz x 
2T 



ImL«(z)— ReL* +1 (z) - Imi* +1 (z)— ReL«(z) 



(228) 



where dimensionless magnetic field h was introduced by ( f77j ) . The phenomeno- 
logical expression (113) can be obtained from this formula by carrying out 
the frequency integration in the same way as was done in the calculation of 



(185) (the essential region of integration is z <C T). 

One can see from ( |228| ) that if ImL^(-z) = -ImL^(z) and ReL^(— z) = 
ReL^(z) the Hall conductivity is equal to zero, or, in terms of the phe- 
nomenological parametrization, the reality of jgl results in a zero Hall ef- 
fect. Physically it is possible to say that this zero is t he d irect consequence 
of electron-hole symmetry. However, from the formula (151) one can see that 
an energy dependence of the density of states or the electron interaction con- 
stant g immediately results in the appearance of an imaginary part of 7gl- 
In the weak interaction approximation 



Im7 GL 



u(0) (d\nT c 



dE 



(229) 



E — Ep 



Usually this value is small in comparison with Re^c l by a ratio of the order 
oiTc/Ep. Taking into account the terms of the order of IrwycL in ( ^28[ ) and 
using the explicit form of the fluctuation propagator for layered supercon- 
ductor (|175|) one can find 



AL c 2rrMlGL h 



2v{Q) 



(230) 
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71/3 dk? 



— 71"/ i 



2tt 



[e + rsin 2 (fc z ,s/2))] 



e + mii%s/2)) 
2 F( 2h } ' 



where 

F(x) = Ax 2 



ip(x) + xip'(x) — 1 — + x) 



(231) 



For H — > the expression for the fluctuation Hall paraconductivity takes the 
form 



Ah 



e 2 T 



{ Im7 G£ 



e + r/2 



e 3 <? r /Ini7 GL \ e + r/2 



6s V f(0) / [e(e + r)] 3/2 6tt S V f(0) / [e (e + r)] 3/2 



r /i.(232) 



One can see that in the 2D case the temperature dependence of the AL 
fluctuation correction to the Hall conductivity 

al „ e 3 ^o ( T, 



° xv " 12tts \E F J e 2 
turns out to be more singular than the MT one. 



6.4 Fluctuations in the ultra-clean case | 175| 



When dealing with the superconductor electrodynamics in the fluctuation 
regime, it is necessary to remember that in the vicinity of the critical tem- 
perature the role of the effective size of a fluctuation Cooper pair is played 
by the GL coherence length £gl(T) = ^o/v^- So, as was already mentioning 
above, the case of a pure enough superconductor with electron mean free 
path £ 3> £o has to be formally subdivided into the clean (£o I <C £gl(T)) 
and ultra-clean (£gl(T) <§C I) limits. The nontrivial cancellation of the con- 



tributions, previously divergent in Tr (see, for example, (192)), will be shown 



in this Section. This results in a reduction of the total fluctuation correction 



in the ultra-clean case to the AL term only. We will base on [175] restricting 
our consideration to the case of a 2D electron system. 

In terms of the parameter Tr, used in the theory of disordered alloys, three 
different domains of the metal purity can be distinguished: Tr <C 1 (dirty 
case), 1 -C Tt <C 1/%/e (clean case) and 1/y/e -C Tt (ultra-clean case of non- 
l ocal electrod ynamics). The latter case was rarely discussed in the literature 



[ 150[l7f ,177j i n spite of the fact that it becomes of primary importance for 



metals of very modest purity, let us say, with Tt ~ 10. Really, in this case 
the condition Tt > 1/%/e , which in terms of the reduced temperature is read 
as 10~ 2 < t « 1, practically covers all the experimentally accessible range of 
temperatures for the fluctuation conductivity measurements. As regards the 
usually considered local clean case (1 <C Tt <§; 1/y/e) for the chosen value 
Tt w 10, it would not have any range of applicability. Indeed, the equivalent 
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condition for the allowed temperature interval is e <C (Tt) -2 , and it almost 
contradicts the 2D thermodynamical Ginzburg-Levanyuk criterion for the 
mean field approximation applicability (Gi(2) = "g^: "C e). Moreover, as we 
will show below, for transport coefficients the higher order corrections become 
comparable with the mean field results much before they are important for 



thermodynamical quantities, namely at e ~ y/Gi^) [ ITS , 179 1 . So in practice 
one can speak about the dirty and the non-local ultra-clean limits only. 

As we saw above the 2D AL contribution turns out to be completely inde- 
pendent of the electron mean free path 1 1| . The anomalous Maki-Thompson 
contribution, being induced by the pairing on the Brownian diffusive trajec- 
tories j7f| , naturally depends on TV, but in an indirect way. It turns out to be 
r— independe nt up to Tt ~ l/y/e, (sec ( |20l| )) and diverges as Trln(Tr) for 



Tt 3> 1/ \[i \ 150| , 176| . The analogous problem takes place in the case of the 



DOS contribution: its standard diagrammatic technique calculations lead to 
a negative correction ( |l92| ) evidently strongly divergent when Tt — > oo. 
In the derivation of all these results the local form of the fluctuation prop- 
agator and Cooperons were used. This is why the direct extension of their 
validity for Tt l/\fe — > oo is incorrect. 



One can notice [ 175 1 that at the upper limit of the clean case, when 
Tt ~ both the DOS and anomalous MT ( |201| ) contributions turn out 

to be of the same order of magnitude but of opposite signs. So one can suspect 
that in the case of a correct procedure of impurity averaging in the ultra-clean 
case the large negative DOS contribution can be cancelled with the positive 
anomalous MT one. In the case of a 2D electron spectrum the Cooperon can 
be calculated exactly for the case of an arbitrary electron mean free path: 

A(q, £l , £2 )=fl- /( 9{ -!%i 5 2 ) (233) 
V TV(ei -e 2 f +v 2 F q 2 l 



One can see that this expression can be reduced to (168) in the case of 
vpq <C \ex — e^l- Let us stress that this result was carried out without any 
expansion over the Cooper pair center of mass momentum q and is valid in 
the 2D case for an arbitrary £q. 

The fluctuation propagator in the 2D case of an arbitrary mean free path 



can be written as 175 



K(q,a-)r^ln| + E{^ (^4) 

n— v ' 



1 I flfc , 1 V , v fV 



' 1 2 4irT 1 4ttTtJ 1 16ir 2 T 2 4ttTt 

Near T c In if- w e and in the local limit, when only small momenta ^« 1 



are involved in the final integrations, the Exp. (234) can be expanded in 



VFq/ max{T, r 1 } and reduces to the appropriate local expression 
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Let us demonstrate the specifics of the non-local calculations for the ex- 
ample of the Maki-Thompson contribution. We restrict our consideration to 
the vicinity of the critical temperature, where the static approximation is 
valid. Using the non-local expressions for the Cooperon and the propagator 
one can find after integration over electronic momentum: 

q(mt) K) = _ 4jTiyv 2 pe 2 T 2 J J^_L(q, 0) x (235) 



where 

M (a,/3,q) 



(2tt)' 



R q (2a)R q (a + f3) - 0{af3)R q (2a)R q (2(3) 
(j3 - a) 2 (R q (2a) - I) (R q {2(3) - i) R q {a + /?) 



and R q (x) — y/x 2 + v^q 2 . The analogous consideration of the DOS diagrams 
2 and 4 which are the leading ones in the clean case [g5| results in similar 
expressions. 

One can see that, after analytical continuation over the external frequency 
uj v — > — and the consequent tending w — > 0, each of the DOS or MT type 
diagrams is written in the form of a Laurent series of the type C_2(TY) 2 + 
C_i(Tt)+Co-|-Ci(Tt) -1 -|-... and is divergent at TV — > oo in accordance with 



(192). Nevertheless the expansion in a Laurent series of the sum of these non- 
local diagrams leads to the exact cancellation of all divergent contributions. 
The leading order of the sum of the MT and DOS contributions in the limit 
of Tt 3> 1 turns out to be proportional to (Tr) _1 only and disappears in the 
ultra-clean limit. So the correct accounting for non-local scattering processes 
in the ultra-clean limit results in a total quantum correction negligible in 
comparison with the AL contribution. Nevertheless, its formal independence 
on imp uriti es concentration (see (|99|)) was re-examined for ultra-clean case 



too in [174] and there it was demonstrated that this statement is valid in a 
rigorous sense only in the case of direct current and absence of a magnetic 
field. Let us recall that the normal Drude conductivity in the ultra-clean case 
takes the form 

G 2 Tll~ I Tfi 

C± M = &xx ± i&xy — z 7J r— , (236) 

1 — iyuo =p uj c )t 

where lo c is the cyclotron frequency. When t — * oo the real part of the 
conductivity vanishes. The analysis of the AL diagram in the ultra-clean 
case demonstrates that each of the Green functions blocks B acquires the 
same denominator. As a result the expression for the fluctuation conductivity 
contains the same Drude like pole but of second order 

AL(l) , . AL(l) 

er ± (CJ) = — rrr (237) 

(1-i(wTw c )t) 2 
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(f^^ is the component of the paraconductivity tensor calculated above in 
the local limit (f)9|-( |232| )). The origin of this pole can be recognized by means 
of the following speculation. The electric held does not interact directly with 
the fluctuation Cooper pairs, but it produces the effect by interaction with the 
quasiparticles forming these pairs only. The characteristic time of the change 
of a quasiparticle state is of the order of t. Consequently the single-particle 
Drude type conductivity in an a.c. field has a first order pole, while in the AL 
paraconductivity it is of second order [174|. In spite of this difference one can 
see that the AL conductivity, like the Drude one, vanishes at tu ^ 0, r — > oo 
because in the absence of impurities the interaction of the electrons does not 
produce any effective force acting on the superconducting fluctuations, while 
the d.c. paraconductivity conserves its usual r— independent form. It is im- 
possible to distinguish the motion of the electron liquid from the condensate 
motion in current experiments without additional scattering. 

The non-local form of the Cooperon and fluctuation propagator have to 
be taken into account not only for the ultra-clean case but in every prob- 
lem where relatively large bosonic momenta are involved: the consideration 
of dynamical and short wavelength fluctuations beyond the vicinity of crit- 
ical temperature, the effect of relatively strong magnetic fields on fluctua- 
tions etc. Recent ly such an approach was developed in a number of studies 
[|!76[|l75|,p0|,^8ll. 



6.5 The effect of fluctuations on the one-electron density of 
states and on tunneling measurements 



Density of states |94|. The appearance of non-equilibrium Cooper pairing 
above T c leads to a redistribution of the one-electron states around the Fermi 
level. A semi-phenomenological study of the fluctuation effects on the density 
of states (DOS) of a dirty superconducting material was first carried out while 
analyzing the tunneling experiments of granular Al in the fluctuation regime 



just above T c [ 182 1 . The second metallic electrode was in the superconducting 
regime and its well developed gap gave a bias voltage around which a struc- 
ture, associated with the superconducting fluctuations of Al, appeared. The 
measured DOS energy dependence has a dip at the Fermi level reaches 
its normal value at some energy Eq(T), show a maximum at an energy value 
equal to several times Eq, finally decreases towards its normal value at higher 
energies. The characteristic energy Eg was found to be of the order of the 
inverse of the GL relaxation time tql introduced above. 

The presence of a depression at E — and of a peak at E ~ (1/ tql) 
in the DOS above T c are precursor effects of the appearance of the super- 
conducting gap in the quasiparticle spectrum at temperatures below T c . The 
microscopic calculation of the fluctuation contribution to the one-electron 
DOS can be carried out within the diagrammatic technique [|93|]94l ■ 



Here we refer the energy E to the Fermi level, where we assume E — 0. 
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Fig. 15. The one-electron Green function with the first order fluctuation correction. 

Let us start from the discussion of a clean superconductor. As is well 
known the one-electron DOS is determined by the imaginary part of the 
retarded Green function integrated over momentum. This definition permits 
us to express the appropriate fluctuation correction in terms of the fluctuation 
propagator: 

8^ c \E,e) = --Im / ^-^8G R (p,E) = --lmR R {E) (238) 

where R R (E) is the retarded analytical continuation of the expression corre- 
sponding to the diagram of Fig. [l5|: 

R(s n )=Tj2 J ^L( q ,n k ) J G 2 (p, £n )G(q-p,a- £n )-(239) 

The result of the integration of the last expression depends strongly of 
the electron spectrum dimensionality: for the two important cases of isotropic 
3-D and 2D electron spectra one finds [M 



Gz (3 , c )Re^= x (240) 



Sv$(E,e) ( 47r) 3 / 2 
, (3) (0) - 7C(3) V-^"~^ T -i_ 2i£ , + xgTc 

l r G l -iE + t g I' 2 [r G l L - 2iE + ^T c ] 1/2 



-^W~ ~ -TW) GH2 ' C) [E* + xIT c t-1] (241) 



E , E + V E 2 + xlT c r G l 

: In ■ 



E 2 + K 2 2 T C T G 1 "VTcTGL 



where xd = nyTrD/7Q($). 

In a dirty superconductor the calculations may be carried out in a similar 
way with the only difference that the impurity renormalization of the Cooper 
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Fig. 16. The theoretical curve of the energy dependence for the normalised correc- 
tion to the single-particle density of states vs energy for a clean two-dimensional 
superconductor above T c . 



vertices has to be taken into account |)3| . The value of the fluctuation dip at 
the Fermi level can be written in the form: 



^ (d) (0) _ / v/G^e- 3 / 2 , D = 3 
1/(0) ~ \Gi^ d) e-\ D = 2 



(242) 



At large energies E ^> t g \ the DOS recovers its normal value, according to 



the same laws (242) but with the substitution e — > E/T c . It is interesting 
that the critical exponents of the fluctuation correction of the DOS change 
when moving from a dirty to a clean superconductor |94[ |: the analysis of 
(p40|)-(p4l|) gives 

S J^M „ { V^^ 1/2 > D = 3 f243 ) 
i/(0) \G» (a , c) e-\ £> = 2 ' ^ J 

Another important respect in which the character of the DOS renormal- 
ization differs strongly for the clean and dirty cases is the energy scale at 
which this renormalization occurs. In the dirty case this energy turns out to 
be H E ( d) ~T-T c ~TqI, while in the clean case £$ c) ~ y/T c (T-T c ) @. 
To understand this important difference one has to study the character of the 
electron motion in both cases |)4| . The relevant energy scale in the dirty case 
is the inverse of the time necessary for the electron to diffuse over a distance 
equal to the coherence length £(T). This energy scale coincides with the in- 
verse relaxation time: tj = T>£~ 2 (T) ~ ~ T — T c . In the clean case, the 
ballistic motion of the electrons gives rise to a different characteristic energy 
scale t- 1 ~ vp^HT) ~ (T c r G 2) 1/2 ~ ^T C {T -T c ). 
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One can check that the integration of ( |24CD -(241) over all positive energies 
gives zero: 

5u{E)dE = (244) 



This "sum rule" is a consequence of a conservation law: the number of quasi- 
particles is determined by the number of cells in the crystal and cannot be 
changed by the interaction. So the only effect which can be produced by the 
inter-electron interaction is a redistribution of the energy levels near the Fermi 



energy The sum rule (244) plays an important role in the understanding of 
the manifestation of the fluctuation DOS renormalization in the observable 
phenomena. As we will see in the next Section the singularity in the tunnel- 
ing current (at zero voltage), due to the density of states renormalization, 
turns out to be much weaker than that in the DOS itself (lne instead of e _1 
or e~ 2 , see ( |242| )-( |24^ )). A similar smearing of the DOS singularity occurs 
in the opening of the pseudo-gap in the c-axis optical conductivity, in the 
NMR relaxation rate etc. These features are due to the fact that we must 
always form the convolution of the DOS with some slowly varying function: 
for example, a difference of Fermi functions in the case of the tunnel current. 
The sum rule then leads to an almost perfect cancellation of the main singu- 
larity at low energies. The main non-zero contribution then comes from the 
high energy region where the DOS correction has its 'tail'. Another impor- 



tant consequence of the conservation law (244) is the considerable increase 
of the characteristic energy scale of the fluctuation pseudo-gap opening with 
respect to Eq: this is eVo = irT for tunneling and ui ~ r _1 for the c-axis 
optical conductivity. 



The effect of fluctuations on the tunnel current [179]. It is quite 
evident that the renormalization of the density of states near the Fermi level, 
even of only one of the electrodes, will lead to the appearance of anomalies 
in the voltage-current characteristics of a tunnel junction. The quasiparticle 
current flowing through it may be written as a convolution of the densities of 
states with the difference of the electron Fermi distributions in each electrode 
(L and R): 

1 

Igp ~ eiWOW(O) X (245) 
( tanh - tanh —J v L {E)u R {E + eV)dE, 

where R n is the Ohmic resistance per unit area and ^l(O), vr(0) are the 
densities of states at the Fermi levels in each of electrodes in the absence of 
interaction. One can see that for low temperatures and voltages the expression 
in parenthesis is a sharp function of energy near the Fermi level. Nevertheless, 
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depending on the properties of the DOS functions, the convolution (245) may 



exhibit different properties. If the energy scale of the DOS correction is much 



larger than T , the expression in parenthesis in (245) acts as a delta- function 
and the zero-bias anomaly in the tunnel conductivity strictly reproduces the 
anomaly of the density of states around the Fermi level: 

sc(y) _ Mey) 

G„(0) ~ u(0) ' 

where G(V) is the differential tunnel conductance and G„(0) is the back- 
ground value of the Ohmic conductance supposed to be bias independent, 
SG(V) = G(V) — G n (Q). T his s ituation, for instance, occurs in a junction with 



one amorphous electrode 183 , where the dynamically screened Coulomb in- 
teraction is strongly retarded, which leads to a considerable suppression of 
the density of states in the vicinity of the Fermi level, within r _1 ^> T. 

It is worth stressing that the proportionality between the tunneling cur- 
rent and the electron DOS of the electrodes is widely accepted as an axiom, 
but generally speaking this is not always so. As one can see from the previous 
subsection, the opposite situation occurs in the case of the DOS renormal- 
ization due to the electron-electron interaction in the Cooper channel: in this 
case the DOS correction varies strongly already in the scale of Eq ~ E^ 0I <C T 



and the convolution in fl245| ) with the DOS (241) has to be carried out without 
the simplifying approximations assumed to obtain ( |246| ) . We will show that 
the fluctuation induced pseudo-gap like structure in the tunnel conductance 



differs drastically from the anomaly of the density of states (241), both in its 
temperature singularity near T c and in the energy range of its manifestation. 

Let us first discuss the effect of the fluctuation suppression of the density 
of states on the properties of a tunnel junction between a normal metal and 
a superconductor above T c . The effect under discussion turns out to be most 
pronounced in the case of thin superconducting films (d <C and layered 

superconductors like HTS cuprates. In order to derive the explicit expression 
for the fluctuation contribution to the differential conductance of a tunnel 



junction with one thin film electrode close to its T c we differentiate (245) 



with respect to volt age, and substitute the DOS correction given by (241) 



This results in (see [179 



"W*>- 1 r . Sv M {Et e)~ (247) 



G„(0) 2T J_ x cosh 2 

„ Gi (2) In ( - \ ) Re< ( - ~ — 

It is important to emphasize several nontrivial features of the result ob- 
tained. First, the sharp decrease (e~ 2 W) of the density of electron states in 
the immediate vicinity of the Fermi level generated by fluctuations surpris- 
ingly results in a much more moderate growth of the tunnel resistance at zero 
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Fig. 17. The theoretical prediction for the fluctuation-induced zero-bias anomaly in 
tunnel-junction resistance as a function of voltage for reduced temperatures e = 0.05 
(top curve), e — 0.08 (middle curve) and e = 0.12 (bottom curve). The insert shows 
the experimentally observed differential resistance as a function of voltage in an 
Al-I-Sn junction just above the transition temperature 



voltage (lnl/e). Second, in spite of the manifestation of the DOS renormal- 



ization at the characteristic scales E, 



id) 



T-Tr 



E, 



(cl) 



V characteristic is much 



the energy scale of the anomaly developed in the / - 
larger: eV = nT > E (see Fig. |l7| ). 

In the inset of Fig. [17| the result of measurements of the differential re- 
sistance of the tunnel junction Al — I — Sn at temperatures slightly above 
the critical temperature of Sn electrode is presented. This experiment was 
done [184] with the purpose of checking the theory proposed [179]. The non- 
linear differential resistance was precisely measured at low voltages which 
permitted the observation of the fine structure of the zero-bias anomaly. The 
reader can compare the shape of the measured fluctuation part of the dif- 
ferential resistance (the inset in Fig. [l7|) with the theoretical prediction. It 
is worth mentioning that the experimentally measured positions of the min- 
ima are eV ~ ±3T C , while the theoretical prediction following from (247) is 
eV = ±7rT c . Recently similar results on an aluminium film with two regions of 
different superconducting transition temperatures were reported [ 185 1 . The 
observations of the pseudogap anomalies in tunneling experiments at tem- 
peratures above T c obtained by a variety of experimental techniques were 
reported in p^p^ , p8t^89Ml9C [ . . 

We will now consider the case of a symmetric junction between two su- 
perconducting electrodes at temperatures above T c . In this case, evidently, 
the correction fl247j ) has to be multiplied by a factor of "two" because of 
the possibility of fluctuation pairing in both electrodes. Furthermore, in view 
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of the extraordinarily weak (~ In 1 /e) temperature dependence of the first 
order correction, different types of high order corrections may manifest them- 
selves on the energy scale eV ~ T - T c or ^T C {T -T c ). Among the m are 
the familiar AL and MT corrections which take place in the first order of Gi 
but in the second order of the barrier transparency. Another type of higher 
order correction appears in the first order of barrier transparency but in the 



second of fluctuation strength (~ Gi ) [179|. Such corrections are generated 



by the interaction of fluctuations through the barrier and they can be evalu- 
ated directly from (245) applied to a symmetric junction . The second order 
correction in Gi can be written as [179]: 



(0,e) 



dE 



cosh 2 (-Mi) 



Gi\i) 



(248) 



This nonlinear fluctuation correction turns out to be small by Gi 2 but its 
strong singularity in temperature and opposite sign with respect to 5G$ 
make it interesting. Apparently it leads to the appearance of a sharp maxi- 
mum at zero voltage in G(V) with a characteristic width eV ~ T — T c in the 
immediate vicinity of T c (one can call this peak as the hyperfine structure). 
This result was confirmed in |191] but to our knowledge such corrections were 
never observed in tunneling experiments. 

\/Gi, 



One can see that SG j, and 6G K f,' become of the same order at e 



(2) 



T fl aiiuuLT /Z 

i.e. the critical region where nonlinear fluctuations effects become important 
in the problem under consideration starts much before the thermodynamical 
criterion e cr ~ Gi. In the next Section we will discuss this early manifestation 
of nonlinear fluctuation effects in transport phenomena. 



6.6 Nonlinear fluctuation effects [14 



As we have already seen in the temperature region Gi <C e <C 1 the thermo- 
dynamic fluctuations of the order parameter \P can be considered to be Gaus- 
sian. Nevertheless the example of the previous Section demonstrates that in 
transport phenomena nonlinear effects, related with the interaction of fluctu- 
ations (higher order c orre ctions) can manifest themselves much earlier. It has 
been found in paper [|l78| , that nonlinear fluctuation phenomena restrict the 
Gaussian region in the fluctuation conductivity of a superconducting film to 



a new temperature scale: y/Gi^d) «f<l (see also ]1 1 0| , |l 1 79|]l 1 5| , |l 92 



In this Section we obtain expressions for the conductivity in the temperature 
region Gi^d) < e < •JGi^d)i where both the perturbation theory works well 
and the nonlinear fluctuation effects are important. 

Let us start from the correlator ( |173| ) which can be expressed by means 
of the Gi(2d) number: 

T 1 327T 3 T 2 

= -7TW = ffl a "STI' (219) 
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The lone-wave- length fluctuations with k 2 < k^ in = %Te/irV can be consid- 
ered as a local condensate. They lead to the formation of the pseudogap 



A 



pg 



d 2 k 

fc2<fe2 iin (2-7T) 2 



1 1/2 



TJGi 



(2d)- 



(250) 



in the single-particle spectrum of excitations. 



Not very close to the transition (e > yJGi^d)) only excitations with ener- 
gies E > A pg are important. The pseudogap does not play any role for them. 
Thus, in this region of temperatures it is sufficient to consider fluctuations 
in the linear approximation only (see [|]J^,Q]). However, in the temperature 
region e < ^JGi^d) the nonlinear fluctuation contribution of the excitations 
with energies E < A pg becomes essential. 

To take into account the spatial dependence of the order parameter we 



will use the results obtained in 193 ] . It was shown there that the spatial vari- 



ations of Apg act on single-particle excitations in the same way as magnetic 
impurities do (the analogy between the effect of fluctuations and magnetic 
impurities was observed in many papers, see for example, [ 194 1 ) . In this case, 
the total pairbreaking rate r can be written as a sum of the pairbreaking 
rate due to the magnetic impurities and the fluctuation term. Thus, the self- 



consistent equation for r can be written in the following form [ 193 1 : 



r 



d 2 k 



(2tt) 2 E- 



\Vk 2 



r 



In the region E < T, T > Te we obtain from Eqs.(251), (B4g) 



r ~T(Gi 



(2d) 



>V2 



A r . 



(251) 



(252) 



which coincides with the results obtained in |11C,101|. 

Let us note, that the pair-breaking rate r was found to be of the order of 
the pseudogap A pg . Thus, a wide maximum appears in the density of states 
at E ~ Apg. As we already saw (201), in purely 2D case the Maki-Thompson 



correction to the conductivity saturates for Te < T (where r = 8Tj v /tt) 
and takes the form fl4|| : 



5a 



MT 



T n . . ttT 

— GiroA\ m . 

r {2d > 8Te 



(253) 



As it can be seen from Eqs.(252) such a saturation takes place when e < 



\JGi(2d)- Similar results have been obtained in |11C,111,115|, with slightly 
different numerical coefficients^]. However, its exact value is not very impor- 
tant since in the region Te < T the Maki-Thompson correction is less singular 



28 Note, that the numerical coefficient in Eq.(253) depends on the definition of 
Gi(2d) an d how the summation of higher order diagrams is made. 
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than the Aslamazov-Larkin one and can be neglected. The latter does not 
saturate when T tends to T c but becomes more and more singular. 

In the presence of the pseudogap if there is no equilibrium, the fluctuat- 
ing Cooper pair lifetime increases with respect to the GL one: Tfi = arch 
(a > 1). Recall, that analogous changes in the coefficient a in the TDGL equa- 
tions appear below the transition temperature (see e.g. p^,^2,195,19€,197|). 



The growth of the coefficient a and, consequently, the increase of the fluc- 
tuation lifetime is because the quasiparticles require more time to attain 
thermal equilibrium (the corresponding time we denote as r e ). A rough es- 
timate gives a ~ A pg T e . In the case of weak energy relaxation, r e has to be 
determined from the diffusion equation taking account of the pseudogap (see 



[ 196 197, 198| ). Note, that in this complicated case the coefficient a becomes a 
non-local operator. Rough estimates give the following value for the thermal 
equilibrium transition time r e ~ (^^min) ~ C^ e ) ■ Taking into account 



Eq.(250) we obtain from ( |ll| ) for the paraconductivity contribution in the 



discussed limit of the weak energy relaxation [|L_ 

p.3/2 

^ ~ qaa (25 4) 

Let us discuss now the role of the energy relaxation processes, character- 
ized by a quasiparticle lifetime t £ . Nonelastic electron scattering off phonons 
and other possible collective excitations can decrease r e significantly. These 
processes together with additional pairbreaking processes (due to magnetic 
impurities or a magnetic field) lead to a decrease of the nonlinear effects. In 
view of these processes, one can write the following interpolation formula for 
the non-linear fluctuation conductivity |l4f| : 

5a - {2d) (255) 



hr 1 + 



Note that Eqs.( 254 -255) are valid only if the parameters r and r e are such 
that the correction to conductivity da is larger than the usual Aslamazov- 
Larkin correction Eq.@. If r > T, Tt £ < y?Gi {2d ) or if T 2 t £ /T < Gi {2d) , 



than nonlinear effects are negligible and the usual result (|99|) is valid for all 
e > Gi(2d)- 

We see that the paraconductivity can exceed the value of the normal 
conductivity a„ in the region Gi^d) < ( < Gi%d)' us reca l^ that in this 
region corrections to all the thermodynamic coefficients are still small and 
the linear theory is well applicable. 



6.7 The effect of fluctuation on the optical conductivity [199] 

The optical conductivity of a layered superconductor can be expressed by 

:«) 

x/3 



the same analytically continued electromagnetic response operator Q^g{u>) 
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(see Exp. (180)) but in contrast to the d.c. conductivity case, calculated with- 
out the assumption to — ► 0. Let us recall that the paraconductivity tensor 
in an a.c. field was already studied in Section 4 in the framework of the 
TDGL equation [[u] and the most interesting asymptotics for our discussion 
(125)-( p"26| ), valid for lo <C T in the 2D regime, were calculated there. The 
microscopic calculation of the AL diagram |)8| shows that in the vicinity of 
T c and for lo <C T the leading singular contribution to the response oper- 
ator Q^g ^ arises from the fluctuation propagators rather than from the 
B a blocks, which confirms the TDGL results. Nevertheless the DOS and MT 
corrections can be calculated only by the microscopic method, as was done 
in |§fi9§.. 

Let us note that the external frequency lo u enters in the expression for 
the DOS contribution to Q a p{u) only by means of the Green's function 
G(p, lj ti+u ) and it is not involved in q integration. So, near T c , even in the 
case of an arbitrary external frequency, we can restrict consideration to the 
static limit, taking into account only the propagator frequency fl^ = 0, and 



to get [19£] 



Recr 



DOS 

a/3 



M = - n K ( W > T ' T ) A <*P In 

Z7TS 



Ve + r + Ve. 



where the anisotropy tensor A a p was introduced in (189). Let us stress that, 
in contrast to the AL frequency dependent contribution, this result has been 
found with only the assumption e <^ 1, so it is valid for any frequency, 



and impurity concentration. The function k(oj,T, t) was calculated in [ 199 1 
exactly but we present here only its asymptotics for the clean and dirty cases: 



k d (w,T«r- 1 ) 



rem 

/T\2 



T< r- 1 < to 



k c i (u,T > r n i 



28C(3) 



(Tr) 2 , lo^t- 1 <T 
(5) 2 , r-i«to«T 
-4(i) 3 ,r-i«T« W 



The general expression for the MT contribution is too cumbersome, so we 
restrict ourselves here to the important 2D overdamped regime (r « e < 7,,): 



MT(an)(2D) ( , A _ 



LO > T, 
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Let us discuss the results obtained. Because of the large number of pa- 
rameters entering the expressions we restrict our consideration to the most 
interesting c-axis component of the fluctuation conductivity tensor in the 2D 
region (above the Lawrence-Doniach crossover temperature). 

The AL contribution describes the fluctuation condensate response to the 
applied electromagnetic field. The current associated with it can be treated as 
the precursor phenomenon of the screening currents in the superconducting 
phase. As was demonstrated above the characteristic "binding energy " of 
fluctuation Cooper pair is of the order of T — T c , so it is not surprising 
that the AL contribution decreases when the electromagnetic field frequency 
exceeds this value. Indeed u AL ~ T — T c is the only relevant scale for a AL : its 
frequency dependence does not contain T, t v and r. The independence from 
the latter is due to the fact that elastic impurities do not present obstacles 
for the motion of Cooper pairs. The interaction of the electromagnetic wave 
with the fluctuation Cooper pairs resembles, in some way, the anomalous 
skin-effect where the reflection is determined by the interaction with the free 
electron system. 

The anomalous MT contribution also is due to fluctuation Cooper pairs, 
but this time they are formed by electrons moving along self-intersecting tra- 
jectories. Being the contribution related with the Cooper pair electric charge 
transfer it does not depend on the elastic scattering time but it turns out to 
be extremely sensitive to the phase-breaking mechanisms. So two character- 
istic scales turn out to be relevant in its frequency dependence: T — T c and 
r" 1 . In the case of HTS, where t^ 1 has been estimated as at least 0.1T C , for 
temperatures up to 5 4- 10 K above T c the MT contribution is overdamped, 
it is determined by the value of t v and is almost temperature independent. 

The DOS contribution to Reer(o;) is quite different from those above. 
In the wide range of frequencies u> <C r _1 the lack of electron states at 
the Fermi level leads to the opposite sign effect in comparison with the AL 
and MT contributions: Re<T DOS (u>) turns out to be negative and this means 
an increase of the surface impedance, or, in other words, decrease of the 
reflectance. Nevertheless, the applied electromagnetic field affects the electron 
distribution and at very high frequencies lo ~ r _1 the DOS contribution 
changes its sign. It is interesting that the DOS contribution, as a one-electron 
effect, depends on the impurity scattering in a similar manner to the normal 
Drudc conductivity. The decrease of Reer (w) starts at frequencies u> ~ 
min{T, r -1 } which for HTS are much higher than T — T c and t~ . 

The w-dependence of Recr'°' with the most natural choice of parameters 
(T c r < T c e < t" 1 < min{T, r^ 1 }) is presented in Fig. ||. 

Let us discuss it referring to a strongly anisotropic layered superconduc- 
tor. The positive AL and MT contributions to cr'"' , being suppressed by the 
square of the interlayer transparency, are small in magnitude and they vary in 
the low frequency region uj ~ min{T — T Cl t^ 1 }. The DOS contribution is pro- 
portional to the first order of transparency and remains in this region almost 



Fluctuation Phenomena in Superconductors 111 



WW] 




Fig. 18. The theoretical dependence [199] of the real part of the conductivity, 
normalized by the Drude normal conductivity, on lj/T, [o"'(w)] = Re [cr(a>)] /u 11 . 
The dashed line refers to the afc-plane component of the conductivity tensor whose 
Drude normal conductivity is = N(0)e 2 rvp. The solid line refers to the c-axis 
component whose Drude normal conductivity is <r" = <r?«7 2 s 2 /v%. In this plot we 
have put Tt = 0.3, E F /T = 50, r = 0.01, e = 0.04, 7% = 4. 



invariable. With a further increase of frequency min{T — T c , t^ 1 } < ui the AL 
and MT contributions decay; Recrj^ remains negative up to w ~ min{T, t -1 }, 
then it changes its sign at u> ~ t , reaches maximum and rapidly decreases. 
The following high frequency behavior is governed by the Drude law. So 
one can see that the characteristic pseudo-gap-like behavior in the frequency 
dependence of the c-axis optical conductivity takes place: a transparency 
window appears in the range w£ [T- T c , r^ 1 ]. 

In the case of the afc-plane optical conductivity the two first positive 
contributions are not suppressed by the interlayer transparency, and exceed 
considerably the negative DOS contribution in a wide range of frequencies. 
Any pseudo-gap like behavior is therefore unlikely in crj.°*(w): the reflectivity 
will be of the metallic kind. 



6.8 Thermoelectric power above the superconducting transition 



203|j201 | 



Thermoelectric effects are difficult both to calculate and to measure if com- 
pared with electrical transport properties. At the heart of the problem lies 
the fact that the thermoelectric coefficients in metals are the small resultant 
of two opposing currents which almost completely cancel. In calculating the 
thermoelectric power one finds that the electrons above the Fermi level carry 
a heat current that is nearly the negative of that carried by the electrons 
below Ep. In the model of a monovalent metal in which band structure and 
scattering probabilities are symmetric about Ep, this cancellation would be 
exact; in a real metal a small asymmetry survives. Because of their compen- 
sated nature, thermoelectric effects are very sensitive to the characteristics 
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of the electronic spectrum, presence of impurities and peculiarities of scatter- 
ing mechanisms. The inclusion of many-body effects, such as electron-phonon 
renormalization, multi-phonon scattering, drag effect, adds even more com- 
plexity to the problem of calculating the thermoelectric power. Among such 
effects, there is also the influence of thermodynamical fluctuations on the 
thermoelectric transport in a superconductor above the critical temperature. 
This problem has been attracting the a ttent ion of theoreticians for more than 



twenty years, since the paper of Maki [202] appeared, where the logarithmi- 



cally divergent AL contribution was predicted for the two-dimensional case. 
So the AL term turns out to be less singular compared with the corresponding 
correction to conductivity. 

In every case where the main AL and MT fluctuation corrections are sup- 
pressed for some reason, the contribution connected with fluctuation renor- 
malization of the one-electron density of states (DOS) can become impor- 
tant. The analogous situation also occurs in the case of the thermoelectric 



coefficient [ 20C , 201 1 . Although the DOS term has the same temperature de- 



pendence as the AL contribution [202 203], it turns out to be the leading 



fluctuation contribution in both the clean and dirty cases, due to its specific 
dependence on the electron mean free path. 

We introduce the thermoelectric coefficient •& in the framework of linear 
response theory as 

1 lm[QW R (w)] 
v = — lim 

T oj^o uj 

where Q^ eh ' R {ui) is the Fourier representation of the retarded correlation 
function of electric J e and heat J h current operators in Hciscnberg represen- 
tation: 

Q(eh)R( X _ X ') = _ t ')(([j h (X), J e {X')})). 

Here X = (r, t) and ((• • •)) represents both thermodynamical averaging and 
averaging over random impurity positions. The correlation function Q( eh ^ R in 
the diagrammatic technique is represented by a bubble with two exact elec- 
tron Green's functions and two external field vertices, the first, ev, associated 
with the electric current operator and the second, \{s n v, associated 



with the heat current operator (e„ is fermionic Matsubara frequency) [ 1 72 ] . 
The first order fluctuation corrections to Q^ eh '{uj v ) are represented by the 
same diagrams as for conductivity (see Fig.||). 

The first diagram describes the AL contribution to thermoelectric coeffi- 



cient and was calculated in [202 203] with the electron-hole asymmetry factor 
taken into account in the fluctuation propagator. Diagrams 2-4 represent the 
Maki-Thompson contribution, neither anomalous nor regular parts of these 



diagrams contribute to $ in any order of electron- hole asymmetry 172, 203 1. 
The contribution from diagrams 5-10 describes the correction to $ due to 
fluctuation renormalization of the one-electron density of states. Evaluating 
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it in the same way as ( |193[ ) but with one heat current vertex one obtains a 
vanishing result if electron- hole asymmetry is not taken into account. The 
first possible source of this factor is contained in the fluctuation propagator; 
it was used in ]203[ for the AL diagram but for the DOS contribution this cor- 
rection results in non-singular contributions to only and can be neglected. 
Another source of electron-hole asymmetry is connected with expansion of 
energy-dependent functions in powers of £/Ep near the Fermi level: 



"(Ov 2 (0="(0)v 2 (0)+£ 



d(K£)v 2 (g)) 



(256) 



Only the second term in Eq. (256) contributes to the thermoelectric coef- 
ficient. Performing the integration over £ , summations over fermionic fre- 
quencies and analytical continuation of the result obtained we find that the 
contribution to the thermoelectric coefficient associated with the DOS renor- 
malization takes the form 

~d{vv 2 )~ 



,QDOS 



where 

k*(Tt) 



1 



Air 2 u(0)v% 



n* (Tt) ln[ 7T —= \ 



8Tt 



Tt 

7C(3) J T 

(Tt) 



4ttTt , 



(257) 



(258) 



9, 4Tr Tt > 1 
^ Tt < 1 



Summing Eq. ( |257| ) with the AL contribution [ 203 one can find the total 
correction to the thermoelectric coefficient in the case of a 2D superconduct- 
ing film of thickness d : 



d DOS + 6 



AL 



= -0.17- 



1 1 



In 



EpT pf<1 \T — T t 



T- 



k,*(T c t) 



5.3 In ^ 



Assuming \xi{u)d/T c ) ~ 2 one finds that the DOS contribution dominates 
the AL one for any value of impurity concentration: n* has a minimum at 
Tt f» 0.3 and even at this point the DOS term is twice as large. In both 
limiting cases Tt <C 1 and Tt ^> 1 this difference strongly increases. 

In practice, although the Seebeck coefficient S = /a is probably the 
easiest to measure among the thermal transport coefficients, the comparison 
between experiment and theory is complicated by the fact that S cannot be 
calculated directly; it is rather a composite quantity of the electrical conduc- 
tivity and thermoelectric coefficient. As both $ and a have corrections due to 
superconducting fluctuations, the total correction to the Seebeck coefficient 
is given by 



(259) 
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We see that the fluctuations result in a decrease of the absolute value of the 
overall Seebeck coefficient as the temperature approaches T c . 

The situation is complicated additionally in HTS materials, where the 
temperature behavior of the background value of the thermoelectric power 
remains unknown. This does not permit to extract precisely from the ex- 
perimental data the fluctuation part Ad to compare it with the theoretical 
prediction. Nevertheless the very sharp maximum in the Seebeck coefficient 

] seems to be unrelated 
to the fluctuation effects. This conclusion is supported by recent analysis 
of the temperature dep endence of the thermoelectric coefficient close to the 
transition in Refs. |2(J7| . 



experimentally observed in a few papers |204|205| : p06| 



6.9 The effect of fluctuations on NMR characteristics [176] 



Preliminaries. In this Section we discuss the contribution of superconduct- 
ing fluctuations to the spin susceptibility and the NMR relaxation rate. For 
both these effects the interplay of different fluctuation contributions is un- 
usual with respect to the case of the conductivity. Like in the case of the 
optical conductivity, the fluctuation contributions to the spin susceptibility 
and the NMR relaxation rate can manifest themselves as the opening of a 
pseudogap already in the normal phase, a phenomenon which is characteristic 
to HTS compounds. 

We begin with the dynamic spin susceptibility x± (k, w) = x± (k, ito v — ► 
i0 + ) where 



uJ 



1/T 



X±(k,iv u ) = j^ dre^ T ((f T (5 + (k,T)^(-k,0)))). (260) 

Here S± arc the spin raising and lowering operators, T T is the time ordering 
operator, and the brackets denote thermal and impurity averaging in the 
usual way. The uniform, static spin susceptibility is given by Xs — X± s ~~ * 
0, uj = 0) while the dynamic NMR relaxation rate is given by 

=±= = lim - / -^3 Im X i fl) (k, w) (261) 
T\T ^owj (2tt) 

where A is a positive constant involving the gyromagnetic ratio. 

For non- interacting electrons ^(k, w^) is determined by the usual loop 
diagram. Simple calculations lead to the well known results for T <C Ep\ 
Xs = v (Pauli susceptibility) and (1/TiT) = Ani> 2 (Korringa relaxation). 
We will present the fluctuation contributions in a dimensionless form by 
normalizing to the above results. 

To leading order in Gi the fluctuation contributions to x± can be discussed 
with the help of the same diagrams drawn for the conductivity in Fig. ^|. It is 
important to note that the role of the external vertices (electron interaction 
with the external field) is now played by the S±(k, t) operators. This means 
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that the two fcrmion lines attached to the external vertex must have opposite 
spin labels (up and down). Consequently, the Aslamazov-Larkin diagram for 
X± does not exist since one cannot consistently assign a spin label to the 
central fermion for spin-singlet pairing. The next set of diagrams to consider 
is the Maki-Thompson contribution. While the MT diagrams for \± appear 
to be identical to those for the conductivity, there is an important difference 
in topology which arises from their spin structure. It is easy to see, by drawing 
the fluctuation propagator explicitly as a ladder of attractive interaction lines, 
that the MT diagram is a non-planar graph with a single fermion loop. In 
contrast the MT graph for the conductivity is planar and has two fermion 
loops. The number of loops, in accordance with the rules of diagrammatic 
technique j9^] , affects the sign of the contribution. 

The diagrams 5 and 6 represent the effect of fluctuations on the single- 
particle self energy, leading to a decrease in the DOS. The DOS diagrams 7 
and 8 include impurity vertex corrections (note that these have only a single 
impurity scattering line as additional impurity scattering in the form of a 
ladder has a vanishing effect). Finally 9 and 10 are the DOS diagrams with 
the Cooperon impurity corrections. 



Spin Susceptibility p08 , 1 76 1 . We note that, when the external frequency 
and momentum can be set to zero at the outset, as is the case for Xs, there 
is no anomalous MT piece (which as we shall see below is the most singu- 
lar contribution to 1/Ti). The MT diagram 2 then yields a result which is 
identical to the sum of the DOS diagrams 5 and 6. 

In the clean limit (T c r 3> 1) the fluctuation contribution is given by 
xf = Xs2 + Xs5+Xs6' au °ther diagrams turn out to be negligible. In the dirty 
case (T c t -c I), the DOS diagrams 5 and 6, together with the regular part of 
the MT diagram (2), yield the same result as in the clean limit (of the order 
0(T c /Ep)). One can see, that this contribution is negligible in comparison 
with the expected dominant one for the dirty case of the order 0(1/ Eft). 
A thorough study of all diagrams shows that the important graphs in the 
dirty case are those with the Cooperon impurity corrections MT 3 and 4, 
and the DOS ones 9 and 10. This is the unique example known to us where 
the Cooperons, which play a central role in the weak localization theory, 
give the leading order result in the study of superconducting fluctuations. 
Diagrams 3 and 4 give one half of the final result given below; diagrams 9 
and 10 provide the other half. The total fluctuation susceptibility is xf = 
Xs3 + Xs4 + Xsg + Xsio- Interesting, that in both the clean and dirty cases 
Xs/Xs^ can be expressed by the same formula if one expresses the coefficient 
in terms of the GL number Gi( 2 ) fll76| ): 

2Gi {2) In ( 2 . ) . (262) 



X S (0) (> \V~t + ^+~r 

It is tempting to explain the negative sign of the fluctuation contribution 
to the spin susceptibility in Eq. (262) as arising from a suppression of the DOS 
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at the Fermi level. But one must keep in mind that only the contribution of 
diagrams 5 and 6 can strictly be interpreted in this manner; the MT graphs 
and the coherent impurity scattering described by the Cooperons do not 
permit such a simple interpretation. 



Relaxation Rate [ p09| , pTo|j204| , |r76| . The calculation of the fluctuation 



contribution to 1/Ti requires rather more care than Xs because of the sub- 
tleties of analytic continuation. Let us define the local susceptibility 

%)= [{dk) x +-Qi,w v ). 

In order to write down the fluctuation contribution to 1/Ti for the case 
of an arbitrary impurity concentration including the ultra-clean case let us 
start from the anomalous MT contribution and evaluate it using the standard 
contour integration techniques 

Urn -\mK^ R (ui) = -%- /(dq)T(q,0)/C(q), (263) 

UJ^O id 8 J 

AC(q) =2t [ 



dz 



-oo cosh 2 (z/4Tr) (yp g 2 _ ( z _ j) 2 _ 



(V*V - (z + i) 2 - l) 



(264) 



We have used the impurity vertices in the general form ( |233| ) . The first simple 
limiting case for ( |264| ) is Iq <§; 1, when the square roots in the denominator 
can be expanded and /C(q) = 27r/Dq 2 . As we already know from Section 
8.4 this corresponds to the usual local approximation and covers the domain 
T c t «1/ y/e. Introducing the pair breaking rate r y v as an infrared cut off one 
can find: 

s{l/Ti) MT {an) (3) x 

(1/TO = — Gl ^— V Hthv) - (265) 
The other limiting case is the "ultra-clean limit" when the characteristic q- 



values satisfy Iq 3> 1. This is obtained when Tr » \j\fk ^> 1. From (264) 
we then find /C(q) = 4\n(lq)/vq, which leads to 



5{l / Tl) MTjan) ^ ^ j_ 

(m)° yi4C(3) Gj(2 -Ve 



Gi [%cl) — ln(Trv^). (266) 



We note that in all cases the anomalous MT contribution leads to an en- 
hancement of the NMR relaxation rate over the normal state Korringa value. 
In particular, the superconducting fluctuations above T c have the opposite 
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sign to the effect for T«T C (where 1/Ti drops exponentially with T). One 
might argue that the enhancement of 1/Ti is a precursor to the coherence 
peak just below T c . Although the physics of the Hebel-Slichter peak (pile-up 
of the DOS just above gap edge and coherence factors) appears to be quite 
different from that embodied in the MT process, we note that both effects 
are suppressed by strong inelastic scattering. 

We now discuss the DOS and the regular MT contributions which are 
important when strong dephasing suppresses the anomalous MT contribution 
discussed above. The local susceptibility arising from diagrams 5 and 6 can 
be easily evaluated. The other remaining contribution is from the regular 
part of the MT diagram. It can be shown that this regular contribution is 
exactly one half of the total DOS contribution from diagrams 5 and 6. All 
other diagrams either vanish (as is the case for graphs 7 and 8) or contribute 
at higher order in I /Eft (this applies to the graphs with the Cooperon 
corrections). The final results can be presented in a unique way for the clean 
(but not ultra-clean) and dirty cases by means of the Gi( 2 ) number : 



The negative DOS contribution to the NMR relaxation rate is evident 
from the Korringa formula and it sign seems very natural while the sign of 
the positive Maki-Thompson contribution can generate a questions about its 
physical origin. Let us consider a self-intersecting trajectory and the motion 
of the electron along it with fixed spin orientation (let us say "spin up"). If, 
after passing a full turn, the electron interacts with the nucleus and changes 
its spin state and momentum to the opposite value it can pass again along 
the previous trajectory moving in the opposite direction . Interaction of the 
electron with itself on the previous stage of the motion is possible due to the 
retarded character of the Cooper interaction and such a pairing process, in 
contrast to the AL one, turns out to be an effective mechanism for relax- 
ation near T c . This purely quantum process opens a new mechanism of spin 
relaxation, and so contributes positively to the relaxation rate 1/Ti. 

In the case of the nuclear magnetic relaxation rate calculations, the elec- 
tron interaction causing nuclear spin flip is considered. If one would try to 
imagine an AL process of this type he would be in trouble, because the 
electron-nuclei scattering with spin-flip evidently transforms the initial sin- 
glet state of the fluctuation Cooper pair in a triplet-one, which is forbidden 
in the scheme discussed. So the formally discovered absence of the AL con- 
tribution to the relaxation rate is evident enough. 

It is worth mentioning that the cancellation of the MT and DOS con- 
tributions to conductivity found in Section 8.4 is crucial for the fluctuation 
contributions to the NMR relaxation rate. In fact, the MT and DOS con- 
tributions here have the same structure as in the conductivity while the AL 
contribution is absent. So the full fluctuation correction to the NMR relax- 
ation rate in clean superconductor simply disappears. 




(267) 
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Discussion. The main results of this Section, valid for e -C 1, can be sum- 
marized as follows: 

(1) Fluctuations lead to a suppression of the spin susceptibility Xs , due 
to the combined effect of the reduction of the single particle density of states 
arising from the self energy contributions, and of the regular part of the MT 
process. 

(2) "Cooperon" impurity interference terms, involving impurity ladders 
in the particle-particle channel, are crucial for the Xs suppression in the dirty 
limit. 

(3) The processes which dominate the results in (1) and (2) above have 
usually been ignored in fluctuation calculations (conductivity, 1/Ti, etc.). 
The spin susceptibility is unusual in that the AL and the anomalous MT 
terms, which usually dominate, are absent. 

(4) For weak pair-breaking (1/t^ <C T c ), an enhance ment of \jT\T , 
coming from the positive anomalous MT term, takes place [ 210| , 176| . 

(5) Strong dephasing suppresses the anomalous MT contribution, and 
1/Ti is then dominated by the less singular DOS and the regular MT terms. 
Being negative, these contributions lead to a suppression of spectral weight 
and a decrease in 1/Ti. 

An intensive controversy took place in recent years in relation to the 
magnetic field dependence of the fluctuation contribution to 1/Ti. The sit- 
uation here resembles much the situation with the magnetoconductivity: 
a positive MT contribution is suppressed by the magnetic field while the 
magnetic field dependent part of the DOS contribution increases with the 
growth of the field. But in contrast to the magnetoconductivity, which can 
be measured extremely precisely, the NMR relaxation rate measurements 
are much more sophisticated. The result of this delicate competition, de- 
pending on many parameters (r, 7 ¥ >,t,), was found in HTS materials to be 
qualitatively different in experiments of various groups. The absence of a 
strong posit ive AL contribution, possible d-pairing, killing the MT contri- 
bution |21C|, small magnitude of the sum of MT and DOS effects even in 
the case of s-pairing, lack of t he precise values of r, 7m, t, leading to contra- 



dictive theoretical predictions[176,113,18C,212], the dispersion in the qual- 



i ty of samples an d expe rimental methods were the reason of this discussion 

p 3[|2 1 ||fl[|||Jy| Me ] . 



7 Conclusions 

Several comments should be made in conclusion. As was mentioned in the 
Introduction the first "fluctuation boom" took place at the end of 60's - 
beginning of 70's, just after the discovery of the fluctuation smearing of the 
superconducting transition and formulation of the microscopic theory of fluc- 
tuations. The discovery of HTS reanimated this interest and, in order to ac- 
count for the specifics of these layered structures with high critical tempera- 
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tures, low charge carrier concentration and other particularities, considerable 
progress in studies of fluctuation phenomena was achieved (see for instance 
the conference proceedings ]217 , |218 and the extensive review article J76|). As 
it is recognized now the optimally or overdoped phases of HTS compounds 
present an example of a "bad" Fermi liquid. The accounting for superconduct- 
ing fluctuations is identical to including of the electron-electron interaction 
beyond the Fermi-liquid approximation. As a result a lot of anomalies of the 
normal state properties of such HTS compounds can be explained. The situ- 
ation was found to be much more sophisticated in underdoped phases where 
the quasiparticle approach, which from we have started this Chapter, fails. 

In the fluctuation theory discussed above, as in modern statistical physics 
in general, two methods have been used mainly: they are the diagrammatic 
technique and the method of functional (continual) integration over the or- 
der parameter. Each of them as we have seen, has its own advantages and 
disadvantages and in different parts of this review we used the former or the 
latter. 

The years of the fluctuation boom coincided with the maximum devel- 
opment of the diagrammatic methods of many body theory in Condensed 
Matter Theory. This methods turns out to be extremely powerful: any phys- 
ical problem, after its clear formulation and writing down the Hamiltonian, 
can be reduced to the summation of some classes of diagrams. The diagram- 
matic technique is especially comfortable for problems containing some small 
parameter. In the theory of superconducting fluctuations such a small pa- 
rameter exists: as we have seen, it is the Ginzburg-Levanyuk number G«m) 
which is expressed as some powers of the ratio max{T c , r~ 1 }/Ep. This is 
why superconducting fluctuations led to the appearance of the small correc- 
tions to different physical values in a wide range of temperatures, and due 
to this smallness these corrections can be evaluated quantitatively. On the 
other hand their specific dependence on nearness to the critical temperature 
T — T c permits to separate them in experiment from other effects. 

In those cases when fluctuations are small it is possible to restrict their 
summation to the ladder approximation only. The diagrammatic technique 
permits in a unique way to describe the quantum and classical fluctuations, 
the thermodynamical and transport effects 

In the description of thermodynamic fluctuations the method of func- 
tional integration turned out to be simpler. The ladder approximation in 
the diagrammatic approach is equivalent to the Gaussian approximation in 
functional integration. The method of functional integration turns out to 
be more effective too in the case of strong fluctuations, for instance, in the 
immediate vicinity of the phase transition. The final equations of the renor- 
malization group carried out by means of functional integrations turn out to 
be equivalent to the result of the summation of the parquet diagrams series. 
Nevertheless the former derivation is much more simple. 
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There is one another reason why we have tried to use both methods and 
even to carry out some results in both ways. In its explosive development of 
the last decades physics became an "oral science" . In the process of such 
direct communication near a blackboard it is difficult to write and to read 
some cumbersome formulas. The language of diagrams is much more com- 
prehensive: by drawing them the speaker demonstrates that this one is small 
and that one has to be taken into account for this and that reasons clear 
for the experienced listener. The success of the diagrammatic technique in 
some sense is similar to the success of geometry in Ancient Greece, where the 
science was "oral" too. 

This advantage of the diagrammatic technique transforms into its dis- 
advantage when there is no direct communication between the speaker and 
listener. It is difficult to learn the diagrammatic technique by a textbook on 
your own, when no one helps you to find the necessary insight on a complex 
graph. May be because of similar reasons geometry disappeared in Middle 
Ages when direct communications between scientists was minimal while the 
"written" algebra had continued to develop. Operating with Osvald Spen- 
gler "prosymbols" we can say that the diagrammatic technique belongs more 
to the Ancient Greece culture style with its "finite body" prosymbol, while 
functional integration, side by side with the Vikings travels to unknown lands 
and Leibnitz analysis of infinitesimals, is an evident modern contribution to 
West-European culture with its "infinite space" prosymbol. 

That is why, suspecting that the modern physics in the near future can 
fall down to a "New Middle Ages" period, we have carried out some results 
by means of functional integration instead of the diagrammatic technique. 
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